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Topological insulators are new states of quantum matter which can not be adiabatically connected to 
conventional insulators and semiconductors. They are characterized by a full insulating gap in the bulk 
and gapless edge or surface states which are protected by time-reversal symmetry. These topologi- 
cal materials have been theoretically predicted and experimentally observed in a variety of systems, 
including HgTe quantum wells, BiSb alloys, and Bi2Te3 and Bi2Se3 crystals. We review theoreti- 
cal models, materials properties and experimental results on two-dimensional and three-dimensional 
topological insulators, and discuss both the topological band theory and the topological field theory. 
Topological superconductors have a full pairing gap in the bulk and gapless surface states consisting 
of Majorana fermions. We review the theory of topological superconductors in close analogy to the 
theory of topological insulators. 
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I. INTRODUCTION 

Ever since the Greeks invented the concept of the atom, fun- 
damental science has focused on finding ever smaller build- 
ing blocks of matter. In the 19th century, the discovery of 
elements defined the golden age of chemistry. Throughout 
most of the 20th century, fundamental science was dominated 



2 



by the search for elementary particles. In condensed matter 
physics, there are no new building blocks of matter to be dis- 
covered: one is dealing with the same atoms and electrons 
as those discovered centuries ago. Rather, one is interested 
in how these basic building blocks are put together to form 
new states of matter. Electrons and atoms in the quantum 
world can form many different states of matter: for example, 
they can form crystalline solids, magnets and superconduc- 
tors. The greatest triumph of condensed matter physics in the 
last century is the classification of these quantu m states by 
the p rinciple of spontaneous symmetry breaking ( Anderson! 
19971) . For example, a crystalline solid breaks translation sym- 
metry, even though the interaction among its atomic build- 
ing blocks is translationally invariant. A magnet breaks rota- 
tion symmetry, even though the fundamental interactions are 
isotropic. A superconductor breaks the more subtle gauge 
symmetry, leading to novel phenomena such as flux quantiza- 
tion and Josephson effects. The pattern of symmetry breaking 
leads to a unique order parameter, which assumes a nonvan- 
ishing expectation value only in the ordered state, and a gen- 
eral effective field theory can be formulated based on the order 
parameter. The ef fective field theory, generally called Landau- 
Ginzburg theory dLandau and Lifshitzl Il980l) . is determined 



by general properties such as dimensionality and symmetry 
of the order parameter, and gives a universal description of 
quantum states of matter. 

In 1980, a new quantum stat e was discovered which d oes 



not fit into this simple paradigm dvon Klitzing et al. 11980b . In 
the quantum Hall (QH) state, the bulk of the two-dimensional 
(2D) sample is insulating, and the electric current is carried 
only along the edge of the sample. The flow of this unidi- 
rectional current avoids dissipation and gives rise to a quan- 
tized Hall effect. The QH state provided the first example of 
a quantum state which is topologically distinct from all states 
of matter known before. The precise quantization of the Hall 
conductance is explained by the fact that it is a topological in- 
variant, which can only take i nteger values in units of e 2 / h, in- 
depen dent of material details ( Laughlinl 1981 ; Thouless et al. 



1982b . Mathematicians have introduced the concept of topo- 
logical invariance to classify different geometrical objects into 
broad classes. For example, 2D surfaces are classified by the 
number of holes in them, or genus. The surface of a perfect 
sphere is topologically equivalent to the surface of an ellip- 
soid, since these two surfaces can be smoothly deformed into 
each other without creating any holes. Similarly, a coffee cup 
is topologically equivalent to a donut, since both of them con- 
tain a single hole. In mathematics, topological classification 
discards small details and focuses on the fundamental distinc- 
tion of shapes. In physics, precisely quantized physical quan- 
tities such as the Hall conductance also have a topological ori- 
gin, and remain unchanged by small changes in the sample. 

It is obvious that the link between physics and topology 
should be more general than the specific case of QH states. 
The key concept is that of a "smooth deformation". In math- 
ematics, one considers smooth deformations of shapes with- 
out the violent action of creating a hole in the deformation 



process. The operation of smooth deformation groups shapes 
into topological equivalence classes. In physics, one can con- 
sider general Hamiltonians of many-particle systems with an 
energy gap separating the ground state from the excited states. 
In this case, one can define a smooth deformation as a change 
in the Hamiltonian which does not close the bulk gap. This 
topological concept can be applied to both insulators and su- 
perconductors with a full energy gap, which are the focus of 
this review article. It cannot be applied to gapless states such 
as metals, doped semiconductors, or nodal superconductors. 
According to this general definition, if we put in contact two 
quantum states belonging to the same topological class, the in- 
terface between them does not need to support gapless states. 
On the other hand, if we put in contact two quantum states be- 
longing to different topological classes, or put a topologically 
nontrivial state in contact with the vacuum, the interface must 
support gapless states. 

From these simple arguments, we immediately see that the 
abstract concept of topological classification can be applied to 
condensed matter system with an energy gap, where the no - 
tion of a smooth deformation can be defined (IZhangl 12008b . 
Further progress can be made through the concepts of topolog- 
ical order parameter and topological field theory (TFT), which 
are powerful tools describing topological states of quantum 
matter. Mathematicians have expressed the intuitive concept 
of genus in terms of an integral, called t opological invari- 
ant, over the local curvature of the surface ( Nakahara , 1990b . 
Whereas the integrand depends on details of the surface ge- 
ometry, the value of the integral is independent of such details 
and depends only on the global topology. In physics, topo- 
logically quantized physical quantities can be similarly ex- 
presse d as invariant integrals over the freque ncy-momentum 



space ( ThoulessL 11998c iThouless et all 11982b . Such quanti 



ties can serve as a topological order parameter which uniquely 
determines the nature of the quantum state. Furthermore, the 
long-wavelength and low-energy physics can be completely 
described by a TFT, leading to powerful predictions of experi- 
mentally measurable topological effects ( Zhangj 1992 ). Topo- 



logical order parameters and TFTs for topological quantum 
states play the role of conventional symmetry-breaking order 
parameters and effective field theories for broken-symmetry 
states. 

The QH states belong to a topological class which ex- 
plicitly breaks time-reversal (TR) symmetry, for example, 
by the presence of a magnetic field. In recent years, 
a new topological class of materials has been theoreti- 



cally predicted and experimentally observed (Bernevig et al. 


20061 Chen et al. . 120091 Fu et al. 


. 2007: Hsiehef al. 


. 2008 


Konig et al. . 2007t Xia et al. . 2009 : Zhang et al. . 


2009b 



These new quantum states belong to a class which is invari- 
ant under TR, and where spin-orbit coupling (SOC) plays 
an essential role. Some important concepts were devel- 



oped in earlier works ( 


Bernevig and Zhangl 2006: Haldane 


1988: KaneE 


ndW 


ele. 


2005t Murakami et al. . 20031 12004 


Sinova et al. , 


2004 


; Zhang and Hul 2001 ), culminating in the 



construction of the topological band theory (TBT) and the 
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TFT of 2D and 3 D topological in s ulators ( Fu and Kanel l2007t 



Fuet all 12007b iKane and Melel 120051: iMoore and BalentsL 
2007t lOi et all 120081: iRovl l2009h . All TR invariant insula- 



tors in nature (without ground state degeneracy) fall into two 
distinct classes, classified by a Z2 topological order parame- 
ter. The topologically nontrivial state has a full insulating gap 
in the bulk, but has gapless edge or surface states consisting 
of an odd number of Dirac fermions. The topological prop- 
erty manifests itself more dramatically when TR symmetry 
is preserved in the bulk but broken on the surface, in which 
case the material is fully insulating both inside the bulk and 
on the surface. In this case, Maxwell's laws of electrodynam- 
ics are dramatically altered by a topological term with a pre- 
cisely quantized coefficient, similar to the case of the QH ef- 
fect. The 2D topological insulator, synonymously called the 
quantum spin Hall (QSH) insulator, w as first theoretically pre- 
dicted i n 2006 (Bernevig et all 20061) and exp erimentally ob- 
served dKonig et ali l2007t iRoth et all. l20Q9f) in HgTe/CdTe 
quantum wells (QW). A topologically trivial insulator state 
is realized when the thickness of the QW is less than a crit- 
ical value, and the topologically nontrivial state is obtained 
when that thickness exceeds the critical value. In the topolog- 
ically nontrivial state, there is a pair of edge states with op- 
posite spins pr opagating in opposit e directions. Four-terminal 
measurements (iKonig et all 120071) show that the longitudinal 
conductance in the QSH regime is quantized to 2e 2 /h, inde- 
pendently of the width of the samp l e. Su bsequent nonlocal 
transport measurements (Ro th et al. . 120091) confirm the edge 
state transport as predicted by theory. The first discovery of 
the QSH topological insulator in HgTe was ranked by Sci- 
ence Magazine as one of the top ten breakthroughs among 
all sciences in year 2007, and the subje ct quickly b ecame 
mainstream in condensed matter physics (IDavi |2008). The 
3D topological insulator was predicted i n the Bii 



r Sb r allo y 



within a certain range of compositions x (IFu and KaneL 120071) . 
and angle-resolved photoemission spectroscopy (ARPES) 
measurements soon observed a n odd number of topolog- 



ically nontrivial surface states dHsieh et all 120081) . Sim 



pier versions of the 3D topological insulator were the oret- 
ically p redicted in Bi 2 Te-j, Sb2Te,c> ( Zhang et al. 120091) and 



Bi 2 Se 3 (IXia et a/.ll2009tlZhang et al. 



2009) compounds with 



a large bulk gap and a gapless surface state consisting of a 
single Dirac cone. ARPES experiments indeed observed the 



linear dispersion relati on of these surface states (C hen et al. 



20091: IXiaef aZll2009l) . These pioneering theoretical and ex 
perimental works opened up the exciting field of topolog- 
ical i nsulators, and the field is now expanding at a rapid 



pace dHasan and Kanel l2010tlKanell2008TlKonig et aZ.|] 2008 



Moorei l2009t lOi and ZhaTg. boiollzhangl . l2008b . Beyond 
the topological materials mentioned above, more than fifty 
new co mpounds have been predicted to be topological in 



Yan et al. 



predicte d to b e topologi cal ln- 
51 iFranzl 120101 iLin et all 1201 (A 



sulators (Chadov et al. , 2010t Franz , 2010t Lin et aZlEoiC: 
2010l) , and two of them have been experim entally 
;centlv dChen et al.[ l2010t ISato et all \20ldti . This 



observed recently 
collective body of work establishes beyond any reasonable 
doubt the ubiquitous existence in nature of this new topo- 
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FIG. 1 Analogy between QH and QSH effects: (a) A spinless ID 
system has both forward and backward movers. These two basic de- 
grees of freedom are spatially separated in a QH bar, as expressed by 
the symbolic equation "2 = 1 + 1". The upper edge supports only a 
forward mover and the lower edge supports only a backward mover. 
The states are robust and go around an impurity without scattering, 
(b) A spinful ID system has four basic degrees of freedom, which 
are spatially separated in a QSH bar. The upper edge supports a for- 
ward mover with spin up and a backward mover with spin down, and 
conversely for the lower edge. That spatial separati on is expressed 
by th e symbolic equation "4 = 2 + 2". Adapted from[Qi and Zhana, 
120101 . 



logical state of quantum matter. It is remarkable that such 
topological effects can be realized in common materials, pre- 
viously used for infrared detection or thermoelectric applica- 
tions, without requiring extreme conditions such as high mag- 
netic fields or low temperatures. The discovery of topological 
insulators has undoubtedly had a dramatic impact on the field 
of condensed matter physics. 

After briefly reviewing the history of the theoretical pre- 
diction and the experimental observation of the topological 
materials in nature, we now turn to the history of the con- 
ceptual developments, and retrace the intertwined paths taken 
by th eorists. An imp ortant step was taken in 1988 by Hal- 
dane dHaldanelll988l). w ho borrowed the concept of the parity 
anomaly ( Redlichl 1984 ; Semenoff , 1984 ) in quantum electro- 
dynamics to construct a theoretical model of the QH state on 
the 2D honeycomb lattice. This model does not require an 
external magnetic field nor the associated orbital quantization 
and Landau levels (LLs). However, it is in the same topo- 
logical class as the ordinary QH states, and requires both two 
dimensionality and the breaking of the TR symmetry. There 
was a misconception at the time that topological quantum 
states could only exist under these conditions. Another im- 
portant step was the construction in 1989 of a TFT of the QH 
effect based on the Chern-Simons (CS) term dZhangl 1 1 992b . 
This theory captures the most important topological aspects 
of the QH effect in a single and unified effective field the- 
ory. At this point, the path towards generalizing the QH 
states became clear: since the CS term can exist in all even 
spatial dimensions, the topological physics of the QH states 
can be generalized to such dimensions. However, it was un- 
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clear at the time what kind of microscopic interactions could 
be resp onsible for these topo logical states. In 2001, Zhang 
and Hu dZhang and Hull200ll) explicitly constructed a micro- 
scopic model for the generalization of the QH state in 4D. 
A crucial ingredient of this model is its invariance under TR 
symmetry, in sharp contrast to the QH state in 2D which ex- 
plicitly breaks TR symmetry. This fact can also be seen di- 
rectly from the CS effective action in 4 + 1 spacetime dimen- 
sions, which is invariant under TR symmetry. With this gener- 
alization of the QH state, two basic obstacles, the breaking of 
TR symmetry and the restriction to 2D, were removed. Partly 
because of the mathematical complexity involved in this work, 
it was not appreciated by the general community at the time — 
but is clear now — that this state is the root state from which 
all T R invariant topo logical insulators in 3D and 2D are de- 
rived (lOi et q/.l l2008). TR invariant topological insulators can 
be classified in the form of a family tree, where the 4D state 
is the "grandfather" state and begets exactly two generations 
of descendants, the 3D and 2D topol ogical insulators, by the 



20081 iRyu e/^l2010tlSchnvder et al 



200 



procedure of dimensional reducti on (Kitae vj 120091: lOi et al 



Motivated by the construction of a TR invariant topolog- 
ical state, theorists started to look for a physical realization 
of this new topological class , and discovered the intrinsic spi n 



Hall effect (Murakami et al, 2003 



Murakami etal. (IMurakami et al II2003I) state their motivation 



2004; Sinova etal. 2004). 



clearly in the introduction: "Recently, the QH effect has been 
generalized to four spatial dimensions [...]. The QH response 
in that system is physically realized through the SOC in a TR 
symmetric system". Soon after, it was realized in 2004 that the 
two key ideas, TR symmetry and SOC, can also be applied to 
ins ulators as well, leading to the concept of spin Hall insula- 
tor ( IMurakami et q/.l 120041) . The spin Hall effect in insulators 
is dissipationless, similarly to the QH effect. The concept of 
spin Hall insulator motivated Kane and Mele in 2005 to inv es- 
tigate the QSH effect in graphene (iKane and Melel 120051) . a 
material first discovered experimentally that same year. Work- 
ing independently, Bernevig and Zhang studied the QSH ef- 
fect in strained semiconductors, where SOC generates LLs 
witho ut the breaking of TR symmetry ( Bernevig and Zhana, 
l2006h . Unfortunately the energy gap in graph ene caused b 



the intrinsic SOC is insignificantly small (IMin et all 1200' 



Yao et all 120071) . Even though neither models have been 



experimentally realized, they played important roles for the 
concep tual developments. In 2006, Bernevig, Hughes and 
Zhang ( Bernevig et al. , 2006b successfully predicted the first 
topological insulator to be realized in HgTe/CdTe QWs. 

The QSH state in 2D can be roughly understood as two 
copies of the QH state, where states with opposite spin 
counter-propagate at the edge. A natural question arises as 
to whether the edge state s of the QSH state are stable. In 
a deeply insightful paper (IKane and Melel 120051) . Kane and 
Mele showed in 2005 that the stability depends on the num- 
ber of pairs of edge states. An odd number of pairs is sta- 
ble, whereas an even number of pairs is not. This obser- 
vation led Kane and Mele to propose a I2 classification of 



TR invariant 2D insulators. In addition, they devised a pre- 
cise algorithm for the computation of a Z2 topological invari- 
ant within TBT TBT was soon extende d to 3D by Fu , Kane 



and Mele, Moore and Balents and Roy dFu and Kanel 12007 ; 
Fu et all 120071: iMoore and Balentsl l2007t iRovl l2009b . where 
sixteen topologically distinct states are possible. Most of these 
states can be viewed as stacked 2D QSH insulator planes, but 
one of them, the strong topological insulator, is genuinely 3D. 
The topological classification according to TBT is only valid 
for noninteracting systems, and it was not clear at the time 
whether these states are stable under more general topological 
deformations including interactions. Qi, Hu ghes and Zhan g 
introduced the TFT of topological insulators dOi et a/.l 120081) . 
and demonstrated that these states are indeed generally sta- 
ble in the presence of interactions. Furthermore, a topolog- 
ically invariant topological order parameter can be defined 
within the TFT as a experimentally measurable, quantized 
topological magnetoelectric effect. The standard Maxwell's 
equations are modified by the topological terms, leading to 
the axion electrodynamics of the topological insulators. This 
work also showed that the 2D and 3D topological insulators 
are descenda nts of the 4D topolog ical insulator state discov- 
ered in 2001 ( IZhang and Hul 120011) . and motivated this series 
of recent developments. At this point, the two different paths 
based on the TBT and TFT converged, and an unified theoret- 
ical framework emerged. 

There are a number of excellent reviews on this sub- 



iect dHasan and Kanell201oHKonig et aZ.ll2008l:lMoorel.l2009t 



Oi and Zhangl 120101) . This article attempts to give a simple 
pedagogical introduction to the subject and reviews the cur- 
rent status of the field. In Sec. HI] and Sec. Hill we review the 
standard models, materials and experiments for the 2D and the 
3D topological insulators. These two sections can be under- 
stood without any prior knowledge of topology. In Sec. |IV] 
we review the general theory of topological insulators, pre- 
senting both the TFT and the TBT. In Sec. [V] we discuss an 
important generalization of topological insulators-topological 
superconductors. 



II. TWO-DIMENSIONAL TOPOLOGICAL INSULATORS 

The QSH state, or the 2D topological insulator was first dis- 
covered in t he HgTe/CdTe quantu m wells. Bernevig, Hughes 
and Zhang ( iBernevig et all 120061) initiated the search for the 
QSH state in semiconductors with an "inverted" electronic 
gap, and predicted a quantum phase transition in HgTe/CdTe 
quantum wells as a function of the thickness c£qw of the quan- 
tum well. The quantum well system is predicted to be a con- 
ventional insulator for cIqw < d c , and a QSH insulator with 
a single pair of helical edge states for cIqw > 4> where d c is 
a critical thickness. The first experimental confirmation of the 
existence of the QSH state i n HgTe/CdTe quantu m wells was 



carried out by Konig et al. dKonig et all 120071) . This work 



reports the observation of a nominally insulating state which 
conducts only through ID edge channels, and is strongly in- 
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fluenced by a TR symme try-breaking magn etic field. Further 
transport measurements ( Roth etai , 2009t) reported unique 
nonlocal conduction properties due to the helical edge states. 

The QSH insulator state is invariant under TR, has a charge 
excitation gap in the 2D bulk, but has topologically pro- 
tected ID gapless edge states that lie inside the bulk insu- 
lating gap. The edge states have a distinct helical prop- 
erty: two states with op posite spin polarization counter- 



propagate at a given edge (iKane and Meld 12005b IWu et al. 
2006t IXu and Moorel 120061) . For this reason, they are also 



called helical edge states, i.e. the sp in is correlated with the 
direction of motion dWu et all 120061) . The edge states come 



in Kramers doublets, and TR symmetry ensures the crossing 
of their energy levels at special points in the Brillouin zone. 
Because of this level crossing, the spectrum of a QSH insu- 
lator cannot be adiabatically deformed into that of a topo- 
logically trivial insulator without helical edge states. There- 
fore, in this precise sense, the QSH insulator represents a new 
topologically distinct state of matter. In the special case that 
SOC preserves a U (l) s subgroup of the full SU(2) spin rota- 
tion group, the topological properties of t he QSH state can b e 
characterized by the spin Chern number dSheng et al. . 2006). 
More generally, the topological properties of the QSH state 
are mathematically charact erized by a Z2 topological invari- 
ant dKane and Melel 120051) . States with an even number of 
Kramers pairs of edge states at a given edge are topologi- 
cally trivial, while those with an odd number are topologically 
nontrivial. The Z2 topological quantum number can also be 
defined for generally interacting systems and experimentally 
measured in term s of the fraction al charge and quantized cur- 
rent on the edge ( Oi et all 2008 ) , and spin-charg e separation 



in the bulk dOi and Zhang[|2008j 



Ran et all 120081) 



In this section, we shall focus on the basic theory of the 
QSH state in the HgTe/CdTe system because of its simplicity 
and experimental relevance, and provide an explicit and ped- 
agogical discussion of the helical edge states and their trans- 
port properties. There are several other theoreti cal propos- 
als fo r the QSH state, including bilayer bismuth dMurakami , 
20061) . and t he "broken-gap" type-II AlSb/InAs/GaSb quan- 



tum wells dLiu et all 120081) . Initial experiments in the 



AlSb/ InAs/GaSb system already show encouraging signa 
tures dKnez et aliuOld) . The QSH system has also been pro 



posed for the transition metal oxide Na2lr03 dShitade et al. 



20091) . The concept of fractio nal QSH state was p r opose d 
at the same time as the QSH ( Bernevig and Zhansa, 2006 ). 
and ha s been investigated theoretically in more details re- 
cently dLevin and Sternll2009tlYoung et a/.ll2008l) . 



A. Effective model of the two-dimensional time-reversal 
invariant topological insulator in HgTe/CdTe quantum 
wells 

In this section we review the basic electronic structure of 
bulk HgTe and CdTe, and present a simple model first in 
traduced by Bernevig, Hughes and Zhang dBernevig et al. 
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FIG. 2 (a) Bulk band structure of HgTe and CdTe; (b) schematic pic- 
ture of quantum w ell geometry and lowe st subbands for two different 
thicknesses. From lBernevig et qZll2006l . 



2006h (BHZ) to describe the physics of those subbands of 
HgTe/CdTe quantum wells that are relevant for the QSH ef- 
fect. HgTe and CdTe crystallize in the zincblende lattice struc- 
ture. This structure has the same geometry as the diamond 
lattice, i.e. two interpenetrating face-centered-cubic lattices 
shifted along the body diagonal, but with a different atom 
on each sublattice. The presence of two different atoms per 
lattice site breaks inversion symmetry, and thus reduces the 
point group symmetry from Oh (cubic) to Td (tetrahedral). 
However, even though inversion symmetry is explicitly bro- 
ken, this only has a small effect on the physics of the QSH 
effect. To simplify the discussion, we shall first ignore this 
bulk inversion asymmetry (BIA). 

For both HgTe and CdTe, the important bands near 
the Fermi level are close to the T point in the Brillouin 
zone [Fig. |2|a)]. They are a s-type band (Tq), and a p-type 
band split by SOC into a J = 3/2 band (T 8 ) and a J = 1/2 
band (TV). CdTe has a band ordering similar to GaAs with a s- 
type (Te) conduction band, and p-type valence bands (Tg, Tj) 
which are separated from the conduction band by a large en- 
ergy gap (~ 1.6 eV). Because of the large SOC present in 
the heavy element Hg, the usual band ordering is inverted: 
the negative energy gap of —300 meV indicates that the Tg 
band, which usually forms the valence band, is above the Tq 
band. The light-hole Tg band becomes the conduction band, 
the heavy-hole band becomes the first valence band, and the 
s-type band (Tq) is pushed below the Fermi level to lie be- 
tween the heavy-hole band and the spin-orbit split-off band 
(TV) [Fig. |2j a)]. Due to the degeneracy between heavy-hole 
and light-hole bands at the T point, HgTe is a zero-gap semi- 
conductor. 

When HgTe-based quantum well structures are grown, the 
peculiar properties of the well material can be utilized to tune 
the electronic structure. For wide QW layers, quantum con- 
finement is weak and the band structure remains "inverted". 
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FIG. 3 Energy levels of the QW as a function of QW width. 
From lKonig et ali l2008t 



However, the confinement energy increases when the well 
width is reduced. Thus, the energy levels will be shifted and, 
eventually, the energy bands will be aligned in a "normal" 
way, if the QW thickness g?qw falls below a critical thickness 
d c . We can understand this heuristically as follows: for thin 
QWs the heterostructure should behave similarly to CdTe and 
have a normal band ordering, i.e. the bands with primarily Tq 
symmetry are the conduction subbands and the Tg bands con- 
tribute to the valence subbands. On the other hand, as c2qw 
is increased, we expect the material to behave more like HgTe 
which has inverted bands. As <iq>w increases, we expect to 
reach a critical thickness where the Tg and Tq subbands cross 
and become inverted, with the Tg bands becoming conduc- 
tion subban ds and the Ts bands becoming valence su bbands 
[Fig. Hb)] dBernevig et ail l2006t lNovik efa/.L I2005I) . The 
shift of energy levels with c?qw is depicted in Fig. [3] The QW 
states derived from the heavy-hole Tg band are denoted by 
H n , where the subscript n = 1,2,3,... describes well states 
with increasing number of nodes in the z direction. Similarly, 
the QW states derived from the electron Tq band are denoted 
by E n . The inversion between E\ and H\ bands occurs at a 
critical thickness gJqw = d c ~ 6.3 nm [Fig. [3]]. In the follow- 
ing, we develop a simple model and discuss why we expect 
QWs with cZqw > d c to form TR invariant 2D topological 
insulators with protected edge states. 

Under our assumption of inversion symmetry, the rele- 
vant subbands, E\ and Hi, must be doubly degenerate since 
TR symmetry is present. We express states in the basis 
{\E 1 +),\H 1 +),\E 1 -),\H 1 -}}, where {E^) and ^i) 
are two sets of Kramers partners. The states \Ei±) and \H±±) 
have opposite parity, hence a Hamiltonian matrix element that 
connects them must be odd under parity. Thus, to lowest or- 
der in k, QEi+), and QEi-), |iTi-)) will each be 
coupled generically via a term linear in k. The \H\+) heavy- 
hole state is formed from the spin-orbit coupled p-orbitals 
\Px+iPy, T), while the \ H\— ) heavy-hole state is formed from 
the spin-orbit coupled p-orbitals | — (p x — ip y ),l). Therefore, 
to preserve rotation symmetry around the growth axis z, the 
matrix elements must by proportional to k± = k x ± ik y . The 



only terms allowed in the diagonal elements are terms that 
have even powers of k including fc-independent terms. The 
subbands must come in degenerate pairs at each k, so there 
can be no matrix elements between the + state and the — state 
of the same band. Finally, if there are nonzero matrix elements 
between | E\ +) , \H\ — ) or \E\ — ) , \H\ +), this would induce a 
higher-order process coupling the ± states of the same band 
and splitting the degeneracy. Therefore, these matrix elements 
are forbidden as well. These simple arguments led to the fol- 
lowing model, 



/ h(k) 
^ h*{-k) 

/ l (k)= £ (k)I 2x2 + <J a (kK, 

where l2 X 2 is the 2x2 identity matrix, and 

e(k) = C-D(kl + k 2 y ), 
d a (k) = (Ak x ,-Ak y ,M(k)) 
M(k) = M - B(kl + k y ), 



(1) 
(2) 



(3) 



where A, B,C, D, M are material parameters that depend on 
the QW geometry, and we choose the zero of energy to be the 
valence band edge of HgTe at k = [Fig. 12). 

The bulk energy spectrum of the BHZ model is given by 



E. 



i(k) ± V 'd a da 



e(k) ± J A\kl + k 2 y ) + M 2 (fc). 



(4) 

(5) 



For B = 0, the model reduces to two copies of the massive 
Dirac Hamiltonian in (2 + 1)D. The mass M corresponds to 
the energy difference between the E\ and Hi levels at the T 
point. The mass M changes sign at the critical thickness d c , 
where E\ and H\ become degenerate. At the critical point, 
the system is described by two copies of the massless Dirac 
Hamiltonian, one for each spin, and at a single valley k 
0. This situation is similar to graphene (iCastro Neto et al. 



20091) . which is also described by the massless Dirac Hamil- 
tonian in (2 + 1)D. However, the crucial difference lies in the 
fact that graphene has four Dirac cones, consisting of two val- 
leys and two spins, whereas we have two Dirac cones, one 
for each spin, and at a single valley. For c?qw > d c , the Ei 
level falls below the H\ level at the T point, and the mass M 
becomes negative. A pure massive Dirac model does not dif- 
ferentiate between a positive or negative mass M . Since we 
are dealing with a nonrelativistic system, the B term is gener- 
ally allowed. In order to make the distinction clear, we call M 
the Dirac mass, and B the Newtonian mass, since it describes 
the usual nonrelativistic mass term with quadratic dispersion 
relation. We shall show later that the relative sign between the 
Dirac mass M and the Newtonian mass B is crucial to deter- 
mine whether the model describes a topological insulator state 
with protected edge states or not. 

HgTe has a crystal structure of the zincblende type which 
lacks inversion symmetry, leading to a BIA term in the Hamil- 
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d(A) 


A (eV A) 


B (eV-A 2 ) 


D (eV) 


M (eV) 


A z (eV) 


55 


3.87 


-48.0 


-30.6 


0.009 


0.0018 


61 


3.78 


-55.3 


-37.8 


-0.00015 


0.0017 


70 


3.65 


-68.6 


-51.2 


-0.010 


0.0016 



TABLE I Material parameters for HgTe/CdTe quantum wells with 
different well thicknesses d. 



tonian, given to leading order by dKonig et ali 120081) 



Hbia 



/ -A 2 \ 

A z 

A z 

\ -A 2 / 



(6) 



This term plays an important role in determining the spin 
orientation of the helical edge state. The topological phase 
transiti on in the presence of BIA has been inves tigated re- 



cently (IKonig et all l2008t Murakami et all 120071) . In addi 



tion, in an asymmetric QW structural inversion symmetry can 
be broken by a build-in electric field, leadi ng to a SOC term of 



Rashba type in the effective Hamiltonian (IRothe et aZl l2010t 



Strom et a/l kOlOl) . For simplicity, we will focus on symmet- 
ric QW witout SIA. In Table H] we give the parameters of the 
BHZ model for various values of <iq>w- 



For the purposes of studying the topological properties of 
this system, as well as the edge states, it is sometimes con- 
venient to work with a lattice regularization of the continuum 
model (Q]i which gives the energy spectrum over the entire 
Brillouin zone, i.e. a tight-binding representation. Since all 
the interesting physics at low energy occurs near the T point, 
the behavior of the dispersion at energies much larger than the 
bulk gap at the T point is not important. Thus, we can choose 
a regularization to simplify our calculations. This simplified 
lattice model consists of replacing ([3]i by 

e(k) = C — 2Da~ 2 (2 — cos k x a — cos k y a), 
d a (k) = [AaT 1 sin k x a, — Aa^ 1 sin k y a, M(k)) , 
M(k) = M — 2Ba~ 2 (2 - cos k x a - cos k y a) . (7) 

It is clear that near the T point, the lattice Hamiltonian reduces 
to the continuum BHZ model in Eq. ((TJ. For simplicity, below 
we work in units where the lattice constant a = 1. 



B. Explicit solution of the helical edge states 

The existence of topologically protected edge states is an 
important property of the QSH insulator. The edge states 
can be obtained by solving the BHZ model (f2]i with an open 
boundary condition. Consider the model Hamiltonian (O de- 
fined on the half-space x > in the xy plane. We can divide 



the model Hamiltonian into two parts, 

H = H a + H i , 

/ M(k x ) Ak x 

Ho = l{k x ) + 



(8) 



Fi = -Dk: 







-M(k 










V 








/ 


-Bkl 








Bkl 










V 












M(k x ) 
Ak x 






-Bkl 






-M(k x ) J 
\ 





Bk 2 y ) 



(9) 



with e(k x ) = C - Dk 2 x and M{k x ) = M - Bk 2 x . All k x - 
dependent terms are included in Hq . For such a semi-infinite 
system, k x needs to be replaced by the operator — id x . On 
the other hand, translation symmetry along the y direction is 
preserved, so that k y is a good quantum number. For k y = 0, 
we have Hi = and the wave equation is given by 

H (k x -> -id x )y(x) = E*(x). (10) 

Since Hq is block-diagonal, the eigenstates have the form 



*t(a=) 





(11) 



where is a two-component zero vector, and &f(x) is related 
to ^i(x) by TR. For the edge states, the wave function if>o(x) 
is localized at the edge and satisfies the wave equation 



e(-id x ) + 



M(-id x ) -iA x d x 
-iA x d x -M(-id x ) 



Ms) = EMV$) 



which has been solved analytically for open boundary 



conditions using different methods ( Konig e t al. , 12008 ; 
Under et all 120091: iLu et al.[ l2010t Izhou et all 120081) . In or- 
der to show the existence of the edge states and to find the re- 
gion where the edge states exist, we briefly review the deriva- 
tion of the explicit form of the edge states by neglecting e for 
simplicity dKonig et q/.l 120081) . 

Neglecting e, the wave equation (fl~2l has particle-hole sym- 
metry. Therefore, we expect that a special edge state with 
E = can exist. With the wave function ansatz i/jq = ipe Xx , 
the above equation can be simplified to 

(M + BX 2 ) t v 4> = AX<j>, 



(13) 



therefore the two-component wave function <j> should be an 
eigenstate of the Pauli matrix r y . Defining a two-component 
spinor <fi± by T y 4>± = ±(f>±, Eq. ( TT3l is simplified to a 
quadratic equation for A. If A is a solution for </> + , then — A is 
a solution for 0_ . Consequently, the general solution is given 
by 

be X2X U+ + (ce~ XlX + de- X2X ^ 



ipo(x) = 
where Ai,2 satisfy 

Al,2 : 



1 

2B 



A±yfA T 



4MB) 



(15) 



The coefficients a, b, c, d can be determined by imposing the 
open boundary condition ip(Q) = 0. Together with the nor- 
malizability of the wave function in the region x > 0, the 
open boundary condition leads to an existence condition for 
the edge states: 3?Ai )2 < (c = d = 0) or 3?Ai )2 > 
(a = b = 0), where !ft stands for the real part. As seen from 
Eq. ( TT5b . these conditions can only be satisfied in the inverted 
regime when M/B > 0. Furthermore, one can show that 
when A/B < 0, we have 3iAi )2 < 0, while when A/B > 0, 
we have 5ftAi 2 > 0. Therefore, the wave function for the edge 
states at the T point is given by 



(a) 



0.05 



-\2X\ 



"* x )tf> + , A/B<0; 



!»_, A/B > 0. 



(16) 



The sign of A/B determines the spin polarization of the edge 
states, which is key to determine the helicity of the Dirac 
Hamiltonian for the topological edge states. Another im- 
portant quantity characterizing the edge states is their decay 
length, which is defined as l c = max {ISftAi^l" 1 }- 

The effective edge model can be obtained by projecting the 
bulk Hamiltonian onto the edge states and fljj, defined in 
Eq. (fTTT) . This procedure leads to a 2 x 2 effective Hamil- 
tonian defined by H<f gc (k y ) = (* a | (h + H^j |^>. To 
leading order in k y , we arrive at the effective Hamiltonian for 
the helical edge states: 



H, 



edge 



= Ak y o z . 



(17) 



For HgTe QWs, we have A ~ 3.6 eV-A dKonig ef 120081) . 
and the Dirac velocity of the edge states is given by v = 
A/h ~ 5.5 x 10 5 m/s. 

The analytical calculation above can be confirmed by exact 
numerical diagonalization of the Hamiltonian (ffl on a strip 
of finite width, which can also include the contribution of 
the e(k) term [Fig. 0). The finite decay length of the heli- 



llHou et al. 


. 2009; Strorr 


and Johannesson 


1 120091 iTeo and KaneL 


20091: Zhou et al. 



C. Physical properties of the helical edge states 

1 . Topological protection of the helical edge states 

From the explicit analytical solution of the BHZ model, 
there is a pair of helical edge states exponentially localized 
at the edge, and described by the effective helical edge the- 
ory (Ffy. In t his co ntext, the concept of "helical" edge 



state ( Wu et al. , 2006I) refers to the fact that states with op- 
posite spin counter-propagate at a given edge, as we see from 
the edge state dispersion relation shown in Fig. Efb), or the 
real space picture shown in Fig.[TJb). This is in sharp contrast 
to the "chirai" edge states in the QH state, where the edge 
states propagate in one direction only, as shown in Fig. [Ha). 




1.02 -0.01 0.01 0.02 

MA" 1 ) 




-0 



- ( U 5 02 -0.01 
k(A" 



0.01 0.02 



FIG. 4 Energy spectrum of the effective Hamiltonian (O in a cylin- 
der geometry. In a thin QW, (a) there is a gap between conduction 
band and valence band. In a thick QW, (b) there are gapless edge 
states on the left and right edge (red and blue lines, respectively). 
The dashed line stands for a typic al value of the chemi cal potential 
within the bulk gap. Adapted from lOi and Zhan3.l20ld . 



In the QH effect, the chirai edge states can not be backseat - 
tered for sample widths larger than the decay length of the 
edge states. In the QSH effect, one may naturally ask whether 
backscattering of the helical edge states is possible. It turns 
out that TR symmetry prevents the helical edge states from 
backscattering. The absence of backscattering relies on the 
destructive interference between all possible backscattering 
paths taken by the edge electrons. 

Before giving a semiclassical argument why this is so, 
we first consider an analogy from daily experience. Most 
eyeglasses and camera lenses have an antireflective coating 
[Fig.|2a)], where light reflected from the top and bottom sur- 
faces interfere destructively, leading to no net reflection and 
thus perfect transmission. However, this effect is not robust, 
as it depends on a precise matching between the wavelength 
of light and the thickness of the coating. Now we turn to the 
helical edge states. If a nonmagnetic impurity is present near 
the edge, it can in principle cause backscattering of the heli- 
cal edge states due to SOC. However, just as for the reflec- 
tion of photons by a surface, an electron can be reflected by 
a nonmagnetic impurity, and different reflection paths inter- 
fere quantum-mechanically. A forward-moving electron with 
spin up on the QSH edge can make either a clockwise or a 
counterclockwise turn around the impurity [Fig. [2b)]. Since 
only spin down electrons can propagate backwards, the elec- 
tron spin has to rotate adiabatically, either by an angle of tt 
or — 7r, i.e. into the opposite direction. Consequently, the two 
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FIG. 5 (a) On a lens with antireflective coating, light reflected by 
top (blue line) and bottom (red line) surfaces interferes destructively, 
leading to suppressed reflection, (b) Two possible paths taken by an 
electron on a QSH edge when scattered by a nonmagnetic impurity. 
The electron spin rotates by 180° clockwise along the blue curve, and 
counterclockwise along the red curve. A geometrical phase factor 
associated with this rotation of the spin leads to destructive interfer- 
ence between the two paths. In other words, electron backscattering 
on the QSH edge is suppressed in a way similar to how the reflec- 
tion of photons is suppres sed by an antireflective coating. Adapted 
from lOi and Zhand.l20ld 



paths differ by a full tt — {—tt) = 2tt rotation of the electron 
spin. However, the wave function of a spin- 1/2 particle picks 
up a negative sign under a full 2tt rotation. Therefore, two 
backscattering paths related by TR always interfere destruc- 
tively, leading to perfect transmission. If the impurity carries 
a magnetic moment, TR symmetry is explicitly broken, and 
the two reflected waves no longer interfere destructively. In 
this way, the robustness of the QSH edge state is protected by 
TR symmetry. 

The physical picture described above ap plies only to the 
case of a single pair o f QSH edge states ( Kane and Mele , 
20051: IWu et all. l2006t IXu and Moorel |2006j). If there are 
two forward-movers and two backward-movers on a given 
edge, an electron can be scattered from a forward-moving to 
a backward-moving channel without reversing its spin. This 
spoils the perfect destructive interference described above, 
and leads to dissipation. Consequently, for the QSH state to 
be robust, edge states must consist of an odd number of for- 
ward (backward) movers. This even-odd effect is the key rea- 
son why the QSH ins ulator is characterized by a Z2 topolog 



ical quantum number dKane and M ele. 2005c IWu et al 
Xuand Moorel 120061). 



2006; 



The general properties of TR symmetry are important for 
understanding the properties of the edge theory. The anti- 
unitary TR operator T takes different forms depending on 
whether the degrees of freedom have integer or half-odd- 
integer spin. For half-odd-integer spin, we have T 2 = — 1 
which implies, by Kramers' theorem, that any single-particle 
eigenstate of the Hamiltonian must have a degenerate part- 
ner. From Fig. Hfb), we see that the two dispersion branches 
at one given edge cross each other at the TR invariant k = 
point. At this point, these two degenerate states exactly satisfy 



Kramers' theorem. If we add TR invariant perturbations to the 
Hamiltonian, we can move the degenerate point up and down 
in energy, but cannot remove the degeneracy. In this precise 
sense, the helical edge states are topologically protected by 
TR symmetry. 

If TR symmetry is not present, a simple "mass" term can 
be added to the Hamiltonian so that the spectrum becomes 
gapped: 



dk 
27 



+ h.c.) , 



where h.c. denotes Hermitian conjugation, and ipt±,ipk± are 
creation/annihilation operators for an edge electron of mo- 
mentum k, with ± denoting the electron spin. The action of 
TR symmetry on the electron operators is given by 



(18) 



which implies 



Consequently, H mass is a TR symmetry breaking perturba- 
tion. More generally, if we define the "chirality" operator 



C = N+ - = 



dk , 



any operator that changes C by 2(2n — 1), 11 E Z is 
odd under TR. In other words, TR symmetry only allows 
2n-particle backscattering, described by operators such as 
ipk+^l'+^p-^p'- (f° r n — !)■ Therefore, the most relevant 
perturbation V'L-V'fc'- i s forbidden by TR symmetry, which is 
essential for the topological stability of the edge states. This 
edge state effective theory is nonchiral, and is qualitatively 
different from the usual spinless or spinful Luttinger liquid 
theories. It can be considered as a new class of I D critical 
theories, dubbed a "helical liquid" dWu et al. . 120061) . Specifi- 
cally, in the noninteracting case no TR invariant perturbation 
is available to induce backscattering, so that the edge state is 
robust. 

Consider now the case of two flavors of helical edge states 
on the boundary, i.e. a ID system consisting of two left- 
movers and two right-movers with Hamiltonian 



\^l s+ vk^ ks+ - ipi s _vkip ks 



A mass term such as m J || (ijj{ 1+ 'ipk2- - V'L-V'm- + h- c 
(with m real) can open a gap in the system while preserving 
time-reversal symmetry. In other words, two copies of the 
helical liquid form a a topologically trivial theory. More 
generally, an edge system with TR symmetry is a nontrivial 
helical liquid when there is an odd number of left- (right-) 
movers, and trivial when there is an even number of them. 
Thus the topology of QSH systems are characterized by a Z2 
topological quantum number. 
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2. Interactions and quenched disorder 



We now review the effect of interacti ons and q i 
disorder on t he Q SH edge liquid dWu et al.l 
Xu and Moorel 120061) . Only two TR invariant nonchiral 



ue nched 
20061: 



interactions can be added to Eq.[l7] the forward and Umklapp 
scatterings 



Hf = g dx ip+ip + ip' < _ip- 



(19) 



9u 



dx e 



-iik F x i t 



tp + (x)ip + (x + a) 



Xip^(x + a)-0_(x) + h.c, 



where the two-particle operators ifj^ifj\tpip are point-split 
with the lattice constant a which plays the role of a short- 
distance cutoff. The chiral interaction terms only renormalize 
the Fermi velocity v, and are thus ignored. It is well known 
that the forward scattering term gives a nontrivial Luttinger 
parameter K = (v — g)/(v + g), but keeps the system gap- 
less. Only the Umklapp term has the potential to open up a gap 
at the commensurate filling k,F = tt/2. The bosonized form 
of the Hamiltonian reads 



H = dx 



g u cos - 



2(7ra) 2 



(21) 



where v = \Jv 2 — g 2 is the renormalized velocity, and we 
define nonchiral bosons <fi = 4>r + 4>l and 8 = 4>r — (f>L, 
respectively, where (f>R and </>l are chiral bosons describ- 
ing the spin up (down) right-mover and the spin down (up) 
left-mover, respectively. <fi contains both spin and charge 
degrees of freedom, and is equivalent to the combination 
<t> c — Q s in the spinful Luttinger liqui d, with (j> c and 8 ? the 
charge and spin bosons, respectively (IGiamarchi , 120031) . It 
is also a compact variable with period y/n. A renormaliza- 
tion group analysis shows that the Umklapp term is relevant 
for K < 1/2 with a pinned value of <fi. Consequently, a gap 
A ~ a^ 1 {g u ) 2 ~ iK opens and spin transport is blocked. The 
mass order parameters N XtV the bosonized form of which is 
N x 



2na 



sinv47r</>, N y 



2-Ka 



cos V 4-7r</>, are odd un- 



K < 3/8 dGiamarchi and Schulzl Il988t IWu et all l2006t 
Xu and Moorel 120061) . At T = 0. N X} y(x) exhibits glassy 



behavior, i.e. disordered in the spatial direction but static in 
the time direction. Spin transport is thus blocked and TR 
symmetry is again spontaneously broken at T = 0. At low 
but finite T, the system remains gapped with TR symmetry 
restored. 

In the above, we have seen that the helical liquid can in 
principle be destroyed. However, for a reasonably weak inter- 
acting system, i.e. K « 1, the one-component helical liquid 
remains gapless. In an Ising ordered phase, the low-energy ex- 
citations on the ed ge are Ising dom ain walls which carry frac- 



(20) tional e/2 charge ( Oi et al. , 2008 ). The properties of multi 



component h elical liquids in the pre sence of disorder has also 
been studied (IXu and Moore , 2006). 

A magnetic impurity on the edge of a QSH insulator is ex- 
pected to act as a local mass term for the edge theory, and 
thus is expected to lead to a suppression of the edge conduc- 
tance. While this is certainly true for a static magnetic im- 
purity, a quantum magneti c impurity, i.e. a Kondo impurity , 
leads to subtler behavior (IMacieiko et al. , l2009t IWu et al 



2006). In the presence of a quantum magnetic impurity, due 



to the combined effects of interactions and SOC one must 
also ge nerally consider local two -particle backscattering pro- 
cesses (IMeidan and Ored, 12005b similar to Eq. d20T >. but oc- 
curring only at the position of the impurity. At high tem- 
peratures, both weak Kondo and weak two-particle backscat- 
tering are expected to give rise to a logar ithmic tempera 



ture d ependence as in the usual Kondo effect (M acieiko et al. 



2009), and their effect is not easily distinguishable. How- 



ever, at low temperatures the physics depends drastically on 
the strength of Coulomb interactions on the edge, parameter- 
ized by the Luttinger parameter K. For weak Coulomb in- 
teractions K > 1/4, the edge conductance is restored to the 
unitarity limit 2e 2 /h with unusual power laws characteristic 
of a "local helical liquid" dMaciejko et al. . 120091: IWu et al. . 
2006h . For strong Coulomb interactions K < 1/4, the con- 
ductance vanishes at T = 0, but is restored at low T by 
a fra ctionalized tunneling current of charge e/2 quasiparti- 
cles ( IMacieiko et al.l 120091) . The tunneling of a charge e/2 
quasiparticle is described by an instanton process which is the 
time counterpart to the stati c e/2 charge on a spatial magnetic 



der TR. For g u < 0, 4> is pinned at either or ^/tt/2, and the 

N y order is Ising-like. At T = 0, the system is in a Ising or- domain wall along the edge dOi et al.ll2008l) . In addition to the 
dered phase, and TR symmetry is spontaneously broken. On 
the other hand, when < T <^ A, N y is disordered, the gap 
remains, and TR symmetry is restored by thermal fluctuations. 
A similar reasoning applies to the case g u > where N x is 
the order parameter. 

There is also the possibility of two-particle backscattering 
due to quenched disorder, described by the term 



single-channel Kondo effect just described, the possibility of 
an even more exotic two-channel Kondo effect on the edge of 
the QSH insulator has also been studied recently dLaw et al. 
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H dis = I dx cos 1 



T&^((f>(x,T) +a(x)), (22) 



where the scattering strength g u (x) and phase a(x) 
are Gaussian random variables. The standard replica 
analysis shows that disorder becomes relevant at 



3. Helical edge states and the holographic principle 

There is an alternative way to understand the qualitative 
difference between an even and odd num ber of edge state s 



in terms of a "fermion doubling" theorem dWu et al.l 120061) . 



This theorem states that there is always an even number of 
Kramers pairs at the Fermi energy for a TR invariant, but 
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FIG. 6 (a) Energy dispersion of a ID TR invariant system. The 
Kramers degeneracy is required at k = and k — ir, so that the 
energy spectrum always crosses 4n times the Fermi level 6f- (b) 
Energy dispersion of the helical edge states on one boundary of the 
QSH system (solid lines). At k = the edge states are Kramers 
partners, while at k — n they merge into the bulk and pair with 
the edge states of the other boundary (dash lines). In both (a) and 
(b), r ed and blue lines re present the two partners of a Kramers pair. 
From lKonig et al. 1 120081 



number of Kramers pairs, the coupling between them can only 
annihilate an even number of Kramers pairs if TR is preserved. 
As a result, at least one pair of gapless edge states can survive. 

This fermion doubling theorem can be generalized to 3D 
in a straightforward way. In the 2D QSH state, the simplest 
helical edge state consists of a single massless Dirac fermion 
in (1 + 1)D. The simplest 3D topological insulator contains 
a surface state consisting of a single massless Dirac fermion 
in (2 + 1)D. A single massless Dirac fermion would also 
violate the fermion doubling theorem and cannot exist in a 
purely 2D system with TR symmetry. However, it can ex- 
ist holographically, as the boundary of a 3D topological in- 
sulator. More generically, there is a one-to-one correspon- 
dence between topological in sulators and robust gap l ess the- 
ories in one lower dimens ion! Freedman et al. , 2010t Kitaev . 
20091: iTeo and Kanell20ldl) 



otherwise arbitrary ID band structure. A single pair of he- 
lical states can occur only "holograhically", i.e. when the 
ID system is the boundary of a 2D system. This fermion 
doubling theorem is a TR invariant generalization of the 
Niels en-Ninomiya no-go theorem f or chiral fermions on a lat- 
tice dNielsen and Ninomivai 119811) . For spinless fermions, 
there is always an equal number of left-movers and right- 
movers at the Fermi level, which leads to the fermion dou- 
bling problem in odd spatial dimensions. A geometrical way 
to understand this result is that for periodic functions (i.e. en- 
ergy spectra of a lattice model), "what goes up must eventually 
come down". Similarly, for a TR symmetric system with half- 
odd-integer spins, Kramers' theorem requires that each eigen- 
state of the Hamiltonian is accompanied by its TR conjugate 
or Kramers partner, so that the number of low-energy channels 
is doubled. A Kramers pair of states at k = must recombine 
into pairs when k goes from to it and 2ir, which requires the 
bands to cross the Fermi level 4n times [Fig.|6ta)]. However, 
there is an exception to this theorem, which is analogous to the 
reason why a chiral liquid can exist in the QH effect. A helical 
liquid with an odd number of fermion branches can occur if 
it is holographic, i.e. if it appears at the boundary (edge) of a 
2D system. In this case, the edge states are Kramers partners 
at k = 0, but merge into the bulk at some finite k c , such that 
they do not have to be combined at k = it. More accurately, 
the edge states on both left and right boundaries becomes bulk 
states for k > k c and form a Kramers pair [Fig.|6jb)]. This 
is exactly the behavior discussed in Sec. III.Bl in the context of 
the analytical solution the edge state wave functions. 

The fermion doubling theorem also provides a physical un- 
derstanding of the topological stability of the helical liquid. 
Any local perturbation on the boundary of a 2D QSH sys- 
tem is equivalent to the action of coupling a "dirty surface 
layer" to the unperturbed helical edge states. Whatever per- 
turbation is considered, the "dirty surface layer" is always ID, 
such that there is always an even number of Kramers pairs of 
low-energy channels. Since the helical liquid has only an odd 



4. Transport theory of the helical edge states 

In conventional diffusive electronics, bulk transport sat- 
isfies Ohm's law. Resistance is proportional to the length 
and inversely proportional to the cross-sectional area, im- 
plying the existence of a local resistivity or conductivity 
tensor. However, in systems such as the QH and QSH 
states, the existence of edge states necessarily leads to non- 
local transport which invalidates the concept of local resis- 
tivity. Such nonlocal transport has been experimentally ob- 
served in t he QH regime in the presence o f a large mag- 
netic field (IBeenakker and van Houtenl Il99ll) . and the non- 
local transport is well described by a quantum transport 
theory based on the Landauer-Biittiker formalism (IBiittiker , 
19881) . A similar transport theory has been developed for 
the helical edge states of the QSH states, and the nonlo- 
cal transpo rt experiments ar e in excellent agreement with 



the theory (IRoth et all 120091) . These measurements are now 



widely acknowledged as constituting definitive experimen- 
tal evid ence for the e xistence of edge states in the QSH 
regime dButtikeriE009l) . 



Within the general Landauer-Biittiker formalism (IBiittiker 



1986), the current- voltage relationship is expressed as 



h 



(23) 



where Ii is the current flowing out of the ith electrode into the 
sample region, Vt is the voltage on the ith electrode, and Tji is 
the transmission probability from the ith to the jth electrode. 
The total current is conserved in the sense that J^, Ii = 0. 
A voltage lead j is defined by the condition that it draws no 
net current, i.e. Ij = 0. The physical currents remain un- 
changed if the voltages on all electrodes are shifted by a con- 
stant amount \i, implying that = ^ n a TR 
invariant system, the transmission coefficients satisfy the con- 
dition Tij = T^. 

For a general 2D sample, the number of transmission chan- 
nels scales with the width of the sample, so that the transmis- 
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sion matrix Tjj is complicated and nonuniversal. However, 
a tremendous simplification arises if the quantum transport is 
entirely dominated by the edge states. In the QH regime, chi- 
ral edge states are responsible for the transport. For a standard 
Hall bar with N current and voltage leads attached, the trans- 
mission matrix elements for the v = 1 QH state are given 
by T(QH)i + i,j = 1, for i = 1,...,N, and all other ma- 
trix elements vanish identically. Here we periodically iden- 
tify the i = N + 1 electrode with i = 1. Chiral edge states 
are protected from backscattering, therefore, the zth electrode 
transmits perfectly to the neighboring (t + l)th electrode on 
one side only. In the example of current leads on the elec- 
trodes 1 and 4, and voltage leads on the electrodes 2, 3, 5 and 
6, (see the inset of Fig. [12] for the labeling), one finds that 
h = -h = hi, V 2 - V 3 = and V x - V 4 = ^I 14 , giving 
a four-terminal resistance of R 



14,23 



= and a two-terminal 

resistance of i?i4.i4 = \. 

The helical edge states can be viewed as two copies of chi- 
ral edge states related by TR symmetry. Therefore, the trans- 
mission matrix is given by T(QSH) = T(QH) + Tt(QH), 
implying that the only nonvanishing matrix elements are given 
by 



T(QSH)i 



T(QSH) 



1. 



(24) 



Considering again the example of current leads on the elec- 
trodes 1 and 4, and voltage leads on the electrodes 2, 3, 5 
and 6, one finds that I\ = —I4 = In, V% — V3 = j^Iu 
and V\ — Va = fr^i4, giving a four-terminal resistance of 
-Ri4,23 = an d a two-terminal resistance of i?i4,i4 = ^5. 
Four terminal resistance with different configurations of volt- 
age and current probes can be predicted in the same way, 
which are all rational fractions ofh/e 2 . The experimental data 
[Fig. Q2) ne atly confirms all th ese highly nontrivial theoretical 
predictions (IRoth et ^Zl l2009h . For two micro Hall bar struc- 
tures that differ only in the dimensions of the area between 
the voltage contacts 3 and 4, the expected resistance values 
-Ri4,23 = 2^ ^14-14 = ^ aie indeed observed for gate 
voltages for which the samples are in the QSH regime. 

As mentioned earlier, one might sense a paradox between 
the dissipationless nature of the QSH edge states and the finite 
four-terminal longitudinal resistance -R14.23, which vanishes 
in the QH state. We can generally assume that the microscopic 
Hamiltonian governing the voltage leads is invariant under 
TR symmetry. Therefore, one would naturally ask how such 
leads could cause the dissipation of the helical edge states, 
which are protected form backscattering by TR symmetry? In 
nature, TR symmetry can be broken in two ways, either at 
the level of the microscopic Hamiltonian, or at the level of 
the macroscopic irreversibility in systems whose microscopic 
Hamiltonian respects TR symmetry. When the helical edge 
states propagate without dissipation inside the QSH insulator 
between the electrodes, neither forms of TR symmetry break- 
ing are present. As a result, the two counter-propagating chan- 
nels can be maintained at two different quasi-chemical poten- 
tials, leading to a net current flow. However, once they enter 



the voltage leads, they interact with a reservoir containing a 
large number of low-energy degrees of freedom, and TR sym- 
metry is effectively broken by the macroscopic irreversibility. 
As a result, the two counter-propagating channels equilibrate 
at the same chemical potential, determined by the voltage of 
the lead. Dissipation occurs with the equilibration process. 
The transport equation d23l breaks the macroscopic TR sym- 
metry, even though the microscopic TR symmetry is ensured 
by the relationship TJy = Tji. In contrast to the case of the 
QH state, the absence of dissipation in the QSH helical edge 
states is protected by Kramers' theorem, which relies on the 
quantum phase coherence of wave functions. Thus, dissipa- 
tion can occur once phase coherence is destroyed in the metal- 
lic leads. On the contrary, the robustness of QH chiral edge 
states does not require phase coherence. A more rigorous and 
microscopic analysis of the different role played b y a metallic 



lead i n QH and QSH states has been performed (R oth et al. 



20091) . the result of which agrees with the simple transport 
equations d23l and d24l) . These two equations correctly de- 
scribe the dissipationless quantum transport inside the QSH 
insulator, and the dissipation inside the electrodes. As shown 
in Sec. III. F. 41 these equations can be put to more stringent 
experimental tests. 

The unique helical edge states of the QSH state can be 
used to construct devices with interesting transport proper- 
ties ( Akhmero v et al. 120091: iKharitonovl (20101 IZhang et all 
20091) . Besides the edge state transport, the Q SH state also 



leads to interesting bulk transport properties (iNovik et al. 
2010h . 



D. Topological excitations 

In the previous sections, we discussed the transport prop- 
erties of the helical edge states in the QSH state. Unlike the 
case of the QH state, these transport properties are not ex- 
pected to be precisely quantized, since they are not directly 
related to the Z2 topological invariant which characterizes the 
topological state. In this section, we show that it is possi- 
ble to measure the Z2 topological quantum number directly 
in experiments. We shall discuss two examples. The first is 
the fractional charge an d quantized curr ent experiments at the 



edge of a QSH system (lOi et all 120081) . Second, we discuss 



the spin -charge separation effect occurring in the bulk of the 
sample dOi and Zhang[|2008l:lRan et adl2008l) . 



1 . Fractional charge on the edge 

The first theoretical proposal we discuss is that of a local- 
ized fractional charge at the edge of a QSH sample when 
a magnetic domain wall is present. The concept of frac- 
tional charge in a condensed matter system induced at a mass 
domain wall go e s bac k to the Su-Schrieffer-Heeger (SSH) 



model dSu et al. Il979h . For spinless fermions, a mass do- 



main wall induces a localized state with one-half of the elec- 



I 1 



tron charge. However, for a real material such as poly- 
acetylene, two spin orientations are present for each elec- 
tron, and because of this doubling, a domain wall in poly- 
acetylene only carries integer charge. The beautiful pro- 
posal of SSH, and its counterpart in fie ld theory, the Jackiw- 
Rebbi model dJackiw and Rebbi. 1 1976b . have never been ex- 
perimentally realized. As mentioned earlier, conventional ID 
electronic systems have four basic degrees of freedom, i.e. 
forward- and backward-movers with two spins. However, a 
helical liquid at a given edge of the QSH insulator has only 
two: a spin up (down) forward-mover and a spin down (up) 
backward-mover. Therefore, the helical liquid has half the de- 
grees of freedom of a conventional ID system, and thus avoids 
the doubling problem. Because of this fundamental topologi- 
cal property of the helical liquid, a domain wall carries charge 
e/2. In addition, if the magnetization is rotated periodically, a 
quantized charge current will flow. This pr ovides a d i rect re - 
alization of the Thouless topological pump ( Thouless , ll983b . 

We begin with the edge Hamiltonian given in Eq. ( fTTI i. 
These helical fermion states only have two degrees of free- 
dom; the spin polarization is correlated with the direc- 
tion of motion. A mass term, being proportional to the 
Pauli matrices a 1 ' 2 ' 3 , can only be introduced in the Hamil- 
tonian by coupling to a TR symmetry breaking external 
field such as a mag netic field, aligned magnetic impuri- 
ties (IGao et q/.l 120091). or interactio n-driven ferromagnetic or- 
der on the edge ( Kharitonov 2010b . To leading order in per- 
turbation theory, a magnetic field generates the mass terms 



H M = dx^ m a (x,t)a a ^ 

^ a=l,2,3 

= / 'crf^U-BiOMK*, (25) 
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FIG. 7 (a) Schematic picture of the half-charge on a domain wall. 
The blue arrows show a magnetic domain wall configuration and the 
purple line shows the mass kink. The red curve shows the charge den- 
sity distribution, (b) Schematic picture of the pumping induced by 
the rotation of magnetic field. The blue circle with arrow shows th e 
rotation of the magnetic field vector. Adapted from lOi et fl/.Ll2008l . 



wall of 6 [Fig. Ha)] dJackiw and Rebbill 1976b . Similarly, the 
charge pumped by a purely time-dependent 9(t) field in a time 
interval [t u t 2 ] is AQ pump |£ = [9(t 2 ) - 0(h)] /2tt. When 9 
is rotated from to 2ir adiabatically, a quantized charge e is 
pumped through the ID system [Fig.|7tb)]. 

From the linear relation m a = t a iBi, the angle 9 can be 
determined for a given magnetic field B. Independent from 
the details of t a i, opposite magnetic fields B and — B always 
correspond to opposite mass, so that 9(B) = 6(— B) + n. 
Thus the charge localized on an anti-phase magnetic domain 
wall of magnetization field is always e/2 mod e, which is a 
direct manifestation of the 7i 2 topological quantum number of 
the QSH state. Such a half charge is dete ctable in a spec ially 
designed single-electron transistor device dOi et qZ.[|2008b . 



where = (^ + , ) and the model-dependent coefficient 
matrix t a i is determined by the coupling of the edge states to 
the magn etic field. According to the w ork of Goldstone and 
Wilczek (IGoldstone and WilczekL 1 1 98 11) . at zero temperature 
the ground-state charge density jo = p and current j\ = j in 
a background field m a (x, t) is given by 



1 



1 



2-7T yjm a m a 



i^e^m a d v m p , a, = 1,2, 



with /j,, v = 0,1 corresponding to the time and space 
components, respectively, and 7713 does not enter the long- 
wavelength charge-response equation. If we parameterize the 
mass terms in terms of an angular variable 9, i.e. mi = 
mcos9, m 2 = ms'm9, the response equation is simplified 
to 



^-d x 9{x,t), j = ~9t6(x,t). 



(26) 



Such a response is topological in the sense that the net charge 
Q in a region [x-y, x 2 ] at time t depends only on the boundary 
values of 9(x, t) i.e. Q — [9(x 2 , t) — 9(x\, t)} /2tt. In par- 
ticular, a half-charge ±e/2 is carried by an anti-phase domain 



2. Spin-charge separation in the bulk 

In addition to the fractional charge on the edge, there have 
been t heoretical proposals for a bulk spin-c harge separation 
effect dOi and Zhana, 12008c iRan et al.l 12008 ). These ideas ar e 



similar to the Z 2 spin pump proposed in dFu and KaneL 120061) . 
We first present an argument which is physically intuitive, but 
only valid when there is at least a U(l) s spin rotation symme- 
try, e.g. when S z is conserved. In this case, the QSH effect is 
simply defined as two copies of the QH effect, with opposite 
Hall conductances of ±e 2 /h for opposite spin orientations. 
Without loss of generality, we first consider a disk geometry 
with an electromagnetic gauge flux of <fo = <j>i = hc/2e, or 
simply 7r in units of h = c = e = 1, through a hole at the 
center [Fig. [8). The gauge flux acts on both spin orientations, 
and the it flux preserves TR symmetry. We consider adia- 
batic processes <fo(t) and 4>\.{t), where <pf(t) = 4>±{t) = 
at t = 0, and </>-]- (t) = 4>i{t) = ±7r at t = 1. Since the flux 
of 7r is equivalent to the flux of — n, there are four different 
adiabatic processes all reaching the same final flux configura- 
tion. In process (a), (pf(t) = —<f>\,{t) and 4>\{t = 1) = 7T. In 
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process (b), fa(t) = —fa(t) and fa(t = 1) = — tt. In pro- 
cess (c), fa (t) — fa(t) and fa(t = 1) = tt. In process (d), 
fa(t) = 4>i(t) and fa(t = 1) = — tt. These four processes 
are illustrated in Fig. [8] Processes (a) and (b) preserve TR 
symmetry at all intermediate stages, while processes (c) and 
(d) only preserve TR symmetry at the final stage. 

We consider a Gaussian loop surrounding the flux. As the 
flux fa(t) is turned on adiabatically, Faraday's law of induc- 
tion states that a tangential electric field E-j- is induced along 
the Gaussian loop. The quantized Hall conductance implies a 

2 

radial current j-j. = ^-z x E-j-, resulting in a net charge flow 
AQ^ through the Gaussian loop: 



dt 



o 
he 



dt- 



dn.j t = -- 



e 2 he 
he 2e 



dt / dl ■ E t 



dt 



(27) 



An identical argument applied to the spin down component 
shows that AQi = — e/2. Therefore, this adiabatic process 
creates the holon state with AQ = AQf + AQi = — e and 

AS Z = AQ t - AQ t = 0. 



t(®?5 





A Q=-e 



(d) 

-it 

AS,=-1/2 





01 = 0! 
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FIG. 8 Four different adiabatic processes from fa = fa = to 
= fa = ±7r. The red (blue) curve stands for the flux fa(i)(t), 
respectively. The symbol "0" ("(g)") represents increasing (decreas- 
ing) fluxes, and the arrows show the current into and out of the Gaus- 
sian loop, induced by the changing flux. Charge is pumped in the 
processes with fa(t) = —<j)i(t), while spi n is pumped in those with 
fa(t) = fa{t). From lOi and Zhangl, 12008 . 



and the charg e quantum numbers are sharpl y defined quan- 
tum numbers (IKivelson and Schrieffen,ll982l) . The insulating 
state has a bulk gap A, and an associated coherence length 
£ ~ A/ A where A is the Dirac parameter in Eq.Q] As long 
as the radius of the Gaussian loop tq far exceeds the coher- 
ence length, i.e., tq ^> £, the spin and the charge quantum 
numbers are sharply defined with exponential accuracy. 

When the spin rotation symmetry is broken but TR symme- 
try is still pre sent, the concept of sp in-charge separation is still 



well defined dOi and Zhangj 12008). A spinon state can be de 



fined as a Kramers doublet without any charge, and a holon or 
a chargeon is a Kramers singlet carrying ch arge ±e. By com- 
bining the spin and charge flux threading dEssin and Moore , 
2007b . it can be shown generally that these spin-charge 
sepa rated quantum numbers are localized near a <fi = tt 
flux dOi and Zhangj 120081: IRan et aZ.ll2008l) . 



E. Quantum anomalous Hall insulator 

Although TR invariance is essential in the QSH insulator, 
there is a TR symmetry breaking state of matter which is 
closely related to the QSH insulator: the quantum anomalous 
Hall (QAH) insulator. The QAH insulator is a band insula- 
tor with quantized Hall conductance but with out orbital mag- 
netic field. Nearly two decades ago, Haldane (IHaldanelll988h 
proposed a model on a honeycomb lattice where the QH is 
realized without any external magnetic field, or the breaking 
of translational symmetry. However, the microscopic mecha- 
nism of circulating current loops within one unit ce ll has not 



been realized in any materials. Qi, Wu and Zhang dOi et al 



20061) proposed a simple model based on the concept of the 
QAH insulator with ferromagnetic moments interacting with 
band electrons via the SOC. This simple m odel can be real 



ized in real mat erials. Two recent proposals dLiu et ail 12008 



Yu et ali [2oTob make use of the properties of TR invariant 
topological insulators to realize the QAH state by magnetic 
doping. This is not accidental, but shows the deep relation- 
ship between these two states of matter. Thus we give a brief 
review of the QAH state in this subsection. 

As a starting point, consider the upper 2x2 block of the 
QSH Hamiltonian ©: 



/i(k)=e(k)I 2x2 +d Q (k)cr a . 



(28) 



If we consider only these two b ands, thi s mode l describes a 



Applying similar arguments to process (b) gives AQ t = TR symmetry breaking system dOi et ali 120061) . As long as 
AQj, = e/2, which leads to a chargeon state with AQ = e there is a gap between t he two bands, the Hall conductance of 
and AS Z = 0. Processes (c) and (d) give AQ t = - AQ± = the system is quantized dThouless e t a/.[ |l982h . The quantized 
e/2 and AQ^ = — AQi = —e/2 respectively, which yield the 
spinon states with AQ = and AS Z = ±1/2. The Hamilto- 
nians H(t) in the presence of the gauge flux are the same at 
t = and t = 1, but differ in the intermediate stages of the 
four adiabatic processes. Assuming that the ground state is 
unique at t = 0, we obtain four final states at t = 1, which are 
the holon, chargeon and the two spinon states. Both the spin 



Hall conductance is determined by the first Chern number of 
the Berry phase gauge field in the Brillouin zone, which, for 
the generic two-band model (l28l . reduces to the following for- 
mula: 



<7 H = 



R 1L 

h Att 



dk„ 



dky d 



ad 

v 



(29) 
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which is e 2 /h times the winding number of the unit vector 
d(k) = d(k)/|d(k)| around the unit sphere. The d(k) vec- 
tor defined in Eq. (01 has a skyrmion structure for M/B > 
with winding number 1, while the winding number is for 
M/B < 0. Just as in an ordinary QH insulator, the system 
with nontrivial Hall conductance e 2 /h has one chiral edge 
state propagating on the edge. For the QSH system described 
by Eq. (|2]), the lower 2x2 block has the opposite Hall conduc- 
tance, so that the total Hall conductance is zero, as guaranteed 
by TR symmetry. The chiral edge state of the QAH and its TR 
partner form the helical edge states of the QSH insulator. 
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FIG. 9 Evolution of band structure and edge states upon increasing 
the spin splitting. For (a) Ge < and Gh > 0, the spin down 
states \E1, — ) andn \H1, — } in the same block of the Hamiltonian 
(O first touch each other, and then enter the normal regime. For (c) 
Ge > and Gh > 0, gap closing occurs between \E1, +) and 
\H1,—}, which belong to different blocks of the Hamiltonian, and 
thus will cross each other without opening a gap, (b) Behav ior of the 
edge states during the level crossing. From lLiu et a/.|,|2008l 

When TR symmetry is broken, the two spin blocks are no 
longer related, and their charge Hall conductances no longer 
cancel exactly. For example, we can consider a different mass 
M for the two blocks, which breaks TR symmetry. If one 
block is in the trivial insulator phase (M/ B < 0) and the other 
block is in the QAH phase (M/ B > 0), the whole system 
becomes a QAH state with Hall conductance ±e 2 /h. Physi- 
cally, this can be realized by exchange coupling with magnetic 
impurities. In a system doped with magnetic impurities, the 
spin splitting term induced by the magnetization is generically 
written as 



( G E 



V o 




Gh 








-Ge 
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-G H J 



(30) 



where Ge and Gh describe the splitting of El and HI bands 
respectively, which are generically different. Adding H s to 



the Hamiltonian (O, we see that the mass term M for the up- 
per block is replaced by M+ (Ge — Gh) /2, while that for the 
lower block is replaced by M — (Ge — Gh) /2. Therefore, the 
two blocks do acquire a different mass, which makes it possi- 
ble to reach the QAH phase. After considering the effect of the 
identity term (Ge + Gji)/% the condition for the QAH phase 
is given by G E G H < 0. When G E G H > and G E ^ G H , 
the two blocks still acquire a different mass, but the system be- 
comes metallic before the two blocks develop an opposite Hall 
conductance. Physically, we can also understand the physics 
from the edge state picture [Fig.|9jb)]. On the boundary of a 
QSH insulator there are counter-propagating edge states car- 
rying opposite spin. When the spin splitting term increases, 
one of the two blocks, say the spin down block, experiences a 
topological phase transition at M = (Ge — Gh)/2- The spin 
down edge states penetrate deeper into the bulk due to the de- 
creasing gap and eventually disappear, leaving only the spin 
up state bound more strongly to the edge. Thus, the system 
has only spin up edge states and evolves from the QSH state 
to the QAH state [Fig. [9{b)] . Although the discussion above 
is based on the specific model (ffj), the mechanism to gener- 
ate a QAH insulator from a QSH insulator is generic. A QSH 
insulator can always evolve into a QAH insulator once a TR 
symmetry breaking perturbation is introduced. 

Fortunately, in Mn-doped HgTe QWs the condition 
GeGh < is indeed satisfied, so that the QAH phase ex- 
ists in this system as long as the Mn spins are polarized. The 
microscopic reason for the opposite sign of Ge and Gh is 
the opposite sign of the s-d and p-d exchange couplings in 
this system ( ILiu et a/.l 120081) . Interestingly, in anot her family 



of QS H insulators, Bi2Se3 and Bi2Te3 thin films dLiu et al. 



201 Oh . the condition GeGh < is also satisfied when mag- 
netic impurities such as Cr or Fe are introduced into the sys- 
tem, but for a different physical reason. In HgTe QWs, the 
two bands in the upper block of the Hamiltonian (0 have 
the same direction of spin, but couple with the impurity spin 
with an opposite sign of exchange coupling because one band 
originates from s-orbitals while the other originates from p- 
orbitals. In Bi2Se3 and Bi2Te3, both bands originate from 
p-orbitals, which have the same sign of exchange coupling 
with the imp urity spin , but th e sign of spin in the upper block 

Consequently, the condition 
More details on the properties 



20101) . 



is opposite ( lYu et al. , 
GeGh < is still satisfied, 
of the Bi2Se3 and Bi2Te3 family of materials can be found 
in the next section, since as bulk materials they are both 3D 
topological insulators. 



F. Experimental results 

1 . Quantum well growth and the band inversion transition 

As shown above, the transition from a normal to an in- 
verted band structure coincides with the phase transition from 
a trivial insulator to the QSH insulator. In order to cover 
both the normal and the inverted band structure regime, HgTe 
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QW samples with a QW width i n the range from 4.5 nm 



to 12.0 nm were grown dKonigl l2007t IKonig et all 12008 



100 



Konig et all I2007I) by molecular beam epitaxy (MBE). Sam- 
ples with mobilities of several 10 5 cm 2 /(V-s), even for low 
densities n < 5 x 10 11 cm -2 , were available for transport 
measurements. In such samples, the mean free path is of the 
order of several microns. For the investigation of the QSH 
effect, devices in a Hall bar geometry [Fig. [T2l inset] of vari- 
ous dimensions were fabricated from QW structures with well 
widths of 4.5 nm, 5.5 nm, 6.4 nm, 6.5 nm, 7.2 nm, 7.3 nm, 
8.0 nm and 12.0 nm. 

For the investigation of the QSH effect, samples with a low 
intrinsic density n(V g = 0) < 5 x 10 11 cm -2 were studied. 
When a negative gate voltage V g is applied to the top gate 
electrode of the device, the usual decrease in electron density 
is observed. In Fig. fTOl a). measurements of the Hall resistance 
R xy are presented for a Hall bar with length L = 600 /im and 
width W — 200 /im. The decrease of the carrier density is 
reflected in an increase of the Hall coefficient when the gate 
voltage is lowered from V to —1 V. In this voltage range, 
the density decreases linearly from 3.5 x 10 11 cm -2 to 0.5 x 
10 11 cm~ 2 [Fig. [lOt b)]. For even lower gate voltages, the 
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FIG. 10 (a) Hall resistance R xy for various gate voltages, indicating 
the transition from n- to p-conductance. (b) Gate-volta ge dependent 
carrie r density deduced from Hall measurements. From lKonig et all 
120081 . 

sample becomes insulating, because the Fermi energy Ep is 
shifted into the bulk gap. When a large negative voltage V g < 
—2 V is applied, the sample becomes conducting again. It can 
be inferred from the change in sign of the Hall coefficient that 
the device is p-conducting. Thus, Ep has been shifted into the 
valence band, passing through the entire bulk gap. 

The peculiar band structure of HgTe QWs gives rise to 
a unique LL dispersion. For a normal band structure, i.e., 
dqw < d c , all LLs are shifted to higher energies for increas- 
ing magnetic fields [Fig. [TTT all. This is the usual behavior 
and can be commonly observed in most materials. When the 
band structure of the HgTe QW is inverted for <1qw > d c , 




FIG. 11 Landau level dispersion for quantum well thicknesses 
of (a) 4.0 nm, (b) 15.0 nm. The qualitative behavior is indica- 
tive f or samples with (a ) normal and (b) inverted band structure. 
From lKonig et aZ.1120081 . 



however, a significant change is observed for the LL disper- 
sion [Fig. []T|(b)]. Due to the inversion of electron-like and 
hole-like bands, states near the bottom of the conduction band 
have predominantly p character. Consequently, the energy 
of the lowest LL decreases with increasing magnetic field. 
On the other hand, states near the top of the valence band 
have predominantly s character, and the highest LL shifts to 
higher energies with increasing magnetic field. This leads to 
a crossing of these two peculiar LLs for a special value of 
the magnetic field. This behavior has been observed earlier 
by the Wiirzburg group and c an now be demonstr ated analyt- 
ically within the BHZ model (IKonig et all 120081) . The exact 
magnetic field S cross at which the crossing occurs depends on 
^qw- The existence of the LL crossing is a clear signature 



of an inverted band structure, which corresponds to a negative 
energy gap with M/B < in the BHZ model. The cross- 
ing of the LLs from the conduction and valence bands can 
be observed in experiments [Fig. [T2la)l. For gate voltages 



V g > -1.0 V and V g < -2.0 V, E F is clearly in the conduc- 
tion band and valence band, respectively. When Ep is shifted 
towards the bottom of the conduction band, i.e. V g < — 1.0 V, 
a transition from a QH state with filling factor v = 1, i.e. 
Rxy = h/e 2 = 25.8 kf2, to an insulating state is observed. 
Such behavior is expected independently of the details of the 
band structure, when the lowest LL of the conduction band 
crosses Ep for a finite magnetic field. When Ep is located 
within the gap, a nontrivial behavior can be observed for de- 
vices with an inverted band structure. Since the lowest LL of 
the conduction band lowers its energy with increasing mag- 
netic field, it will cross Ep for a certain magnetic field. Subse- 
quently, one occupied LL is below Ep, giving rise to the usual 
transport signatures of the quantum Hall regime, i.e. R xy is 
quantized at h/e 2 and R xx vanishes. When the magnetic field 
is increased, the LLs from the valence and conduction band 
cross. Upon crossing, their "character" is exchanged, i.e. the 
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FIG. 12 (a) Hall resistance R xy of a (L x W) = (600 x 200) /im 2 
QW structure with 6.5 nm well width for different carrier concen- 
trations obtained for different gate voltages V g in the range from 
— 1 V to —2 V. For decreasing V g , the n-type carrier concentration 
decreases and a transition to p-type conduction is observed, passing 
through an insulating regime between —1.4 V and —1.9 V at zero 
field, (b) Landau level fan chart of a 6.5 nm quantum well obtained 
from an eight-band k • p calculation. Black dashed lines indicate 
the position of the Fermi energy, Ef, for gate voltages —1.0 V and 
—2.0 V. Red and green dashed lines indicate the position of Ep for 
the red and green Hall resistance traces in (a). The crossing points 
of Ef with the re spective Landau lev els are marked by arrows of the 
same color. From Konig et ai, 20071 . 



level from the valence band turns into a conduction band LL 
and vice versa. The lowest LL of the conduction band now 
rises in energy for larger magnetic fields. Consequently, it 
will cross the Ep for a certain magnetic field. Since Ep will 
be located within the fundamental gap again afterwards, the 
sample will become insulating again. Such a reentrant n-type 
QH state is shown in Fig. fTZl a) for V g = —1.4 V (green trace). 
For lower gate voltages, a corresponding behavior is observed 
for a p-type QH state (e.g. red trace for V g = —1.8 V). As 
Fig.[T2|b) shows, the experimental results are in good agree- 
ment with the theoretically calculated LL dispersion. The 
crossing point of the LLs in magnetic field, £> cross , can be 
determined accurately by tuning Ep through the energy gap. 
Thus, the width of the Q W layer can be verified experimen- 



tally ( IKonig etall 120071) 



for the existence of the QSH state. In contrast, trivial insu- 
lating behavior is obtained for devices with dqw < d c . For 
a normal band structure, the energy gap between the lowest 
LLs of the conduction and valence bands increases in mag- 
netic field [Fig. [TTJa)]. Thus, a sample remains insulating in 
magnetic field, if Ep is located in the gap at zero field. The 
details of the physics of this reentrant QH state can be under- 
stood witWnJheBHZmodel with an added orbital magnetic 
field (IKonig et all 120081) . This nontrivial LL c rossing could 



also be detected optically ( Schmidt et al. , 20091) . 



2. Longitudinal conductance in the quantum spin Hall state 

Initial evidence for the QSH state was revealed when Hall 
bars of dimensions (L x W) = (20.0 x 13.3) ^m 2 with dif- 
ferent thickness dqw are studied. For thin QW devices with 
(Iqw < d c and a normal band structure, the sample shows 
trivial insulating behavior [Fig. \l3l . A resistance of several 
megaohms is measured when the Fermi level lies within the 
bulk insulating gap. This value can be attributed to the noise 
level of the measurement setup, and the intrinsic conductance 
is practically zero. For a thicker device with g?qw > d c 
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The observation of a reentrant QH state is a clear indication 
of the nontrivial insulating behavior, which is a prerequisite 



FIG. 13 Longitudinal resistance of a 4.5 nm QW [dashed (black)] 
and a 8.0 nm QW [so lid (red)] as a function of gate voltage. 
From lKonig etall I2008L 

and an inverted band structure, however, the resistance does 
not exceed 100 kO. This behavior is reproduced for vari- 
ous Hall bars with a QW width in the range from 4.5 nm to 
12.0 nm. While devices with a normal band structure, i.e. 
dqw < d c ~ 6.3 nm, show trivial insulating behavior, a finite 
conductance in the insulating regime is observed for samples 
with an inverted band structure. 

The obtained finite resistance R pa 100 kfl is significantly 
higher than the four-terminal resistance h/(2e 2 ) ps 12.9 kfl 
one anticipates for the geometry used in the experiments. 
The enhanced resistance in these samples with a length of 
L = 20 fim can be understood as a consequence of inelas- 
tic scattering. While, as discussed above, the helical edge 
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states are robust against single-particle elastic backscatter- 
ing, inelastic mechanisms can cause backscattering. For n- 
doped HgTe quantum wells, the typical mobility of the or- 
der of 10 5 cm 2 / (V-s) implies an elast ic mean free path of the 



order of 1 fm\ (IDaumer et ail 120031) . Lower mobilities can 



be anticipated for the QSH regime. The inelastic mean free 
path, which determines the length scale of undisturbed trans- 
port by the QSH edge states, can be estimated to be several 
times larger due to the suppression of phonons and the re- 
duced electron-electron scattering at low temperatures. Thus, 
the inelastic scattering length is of the order of a few microns. 

For the observation of the QSH conductance, the sample 
dimensions were reduced below the estimated inelastic mean 
free path. When Hall bars with a length L = 1 /im are stud- 
ied, a four-terminal resistance close to h/(2e 2 ) is observed. 
The threshold voltage Vth 1S defined such that the QSH regime 
is in the vicinity of V g = Vth- The slight deviation of R 



3. Magnetoconductance in the quantum spin Hall state 

Another indication that the observed nontrivial insulating 
state is caused by the QSH effect is obtained by measure- 
ments in a magnetic field. The following experimental results 
were obtained on a Hall bar with dimensions (L x W) = 
(20.0 x 1 3.3) /xm 2 in a vector magnet syst em at a temperature 
of 1.4 K dKonigl, 120071 iKonig et al[\200lh . When a magnetic 
field is applied perpendicular to the QW layer, the QSH con- 
ductance decreases significantly already for small fields. A 
cusp-like magnetoconductance peak is observed with a full 
width at half-maximum £>fwhm °f 28 ml Additional mea- 
surements show that the width of the magnetoconductance 
peak decreases with decreasing temperature. For example, 
£>fwhm = 10 mT is observed at 30 mK. For various devices 
of different sizes, a qualitatively similar behavior in magnetic 
field is observed. 




FIG. 14 Longitudinal resistance as a function of gate voltage for two 
devices with L = 1 ^im. The width W is 1 /im [solid (black) and 
dotted (blue)] and 0.5 ^im [dashed (red)]. The solid and dashed traces 
were obtained at a tempe rature of 1.8 K, and the dotted one at 4.2 K. 
From lKonig et q/.Ll2008t 



from the quantized value h/(2e 2 ) can be attributed to some 
residual scattering. This is an indication that the length of the 
edge states still exceeds the inelastic mean free path. The re- 
sults presented in Fig. [14] provide evidence that transport in 
the QSH regime indeed occurs due to edge states. The two 
devices with W = 1.0 fjm and W — 0.5 /jm were fab- 
ricated from the same QW structure. The resistance of the 
two devices differ significantly in the n-conducting regime, 
where transport is determined by bulk properties. In the QSH 
regime, however, both devices exhibit the same resistance, 
even though the width of the devices differs by a factor of 
two. This fact clearly shows that the conductance is due to the 
edge states, which are independent of the sample width. 




0,10 



FIG. 15 Four-terminal magnetoconductance G14.23 in the QSH 
regime as a function of tilt angle between the QW plane and the ap- 
plied magnetic field for ad = 7.3 nm QW structure with dimensions 
(L x W) — (20 x 13.3) /1m 2 measured in a ve ctor field cryostat at 
a temperature of 1.4 K. From lKonig et a/.Ll2008l . 



When the magnetic field is tilted towards the plane of the 
QW, the magnetoconductance peak around B = widens 
steadily [Fig. [T5l . For a tilt angle a = 90°, i.e. when 
the magnetic field is in the QW plane, only a very small 
decrease in the conductance is observed. The decrease of 
the conductance for an in-plane field can be described by 
Bfwhm ~ 0.7 T for any in-plane orientation. From the re- 
sults shown in Fig. [15] it is evident that a perpendicular field 
has a much larger influence on the QSH state than an in-plane 
field. The magnetoresist ance in the QSH regime has been in- 
vestigated theo r eticall y ( Chu et al\ 2009t Konig et al. , 2008 ; 
Macieiko et al. 120091: iTkachov and Hankiewiczl l2010t) . The 
large anisotropy can be understood by a slightly modified ver- 
sion of the BHZ model with the inclusion of BIA terms and 



aniso tropy in the g-factor (IKonig et a/.L l2008: Ma ciejko et al. 
20091) . The cusp behavior in the magnetoconductance is pos 



sibly due to the presence of strong disorder; numerical simu- 
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FIG. 16 Experimental measurements of the four- and two-terminal 
resistance: (a) i?i4,23 (red line) and 7?i4 ; i4 (green line) and (b) 
^13,56 (red line) and -R13.13 (green line). The dotted blue lines indi- 
cate the exp ected resistance va lue from the theory of the helical edge 
states. From lRofh et al. 



lations ( iMacieiko et al. . 120091) are in good agreement with the 
experimental results. 
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FIG. 17 Nonlocal four-terminal resistance and two- terminal resis- 
tance measured on an H-bar device: i?i4,23 (red line) and 7?i4,i4 
(green line). The dott ed blue line repre sents the theoretically ex- 
pected resistance. From lRofh et q/.L |2009| . 



4. Nonlocal conductance 

Further confidence in the helical edge state transport can 
be gain ed by performing m ore extended multi-terminal exper- 



iments ( Roth et al . 2009). The longitudinal resistance of a 



device was measured by passing a current through contacts 
1 and 4 [Fig. [161 and by detecting the voltage between con- 
tacts 2 and 3 (i?i4,23)- For this case, a result similar to the re- 
sults found previously, i.e. a resistance hj (2e 2 ) when the bulk 
of the device is gated into the insulating regime [Fig. [ToTa)l. 
However, the longitudinal resistance is significantly differ- 
ent in a slightly modified configuration, where the current is 
passed through contacts 1 and 3 and the voltage is measured 
between contacts 4 and 5 (^13,45) [Fig. [ToTb)l. The result 
is i?i3,45 « 8.6 kfi, which is markedly different from what 
one would expect for either QH transport, or purely diffusive 
transport, where this configuration would be equivalent to the 
previous one. However, the application of the transport equa- 
tions d23l and d24i i indeed predicts that the observed behavior 
is what one expects for helical edge channels. One easily finds 
that this resistance value can be expressed as an integer frac- 
tion of the inverse conductance quanta e 2 /h: Ri3,45 = h/ie 2 . 
This result shows that the current through the device is influ- 
enced by the number of ohmic contacts in the current path. As 
discussed earlier, these ohmic contacts lead to the equilibra- 
tion inside the contact of the chemical potentials of the two 
counter-propagating helical edge channels. 

Another measurement that directly confirms the nonlocal 
character of the helical edge channel transport in the QSH 
regime is shown in Fig. [T7] This figure shows data obtained 
from a device in the shape of the letter "H". In this four- 
terminal device the current is passed through contacts 1 and 
4 and the voltage is measured between contacts 2 and 3. In 
the metallic n-type regime (low gate voltage) the voltage sig- 



nal tends to zero. In the QSH regime, however, the nonlo- 
cal resistance signal increases to « 6.5 kfi, which again fits 
perfectly to the result of Laudauer-Biittiker considerations: 
-^14,23 = h/4e 2 sa 6.45 kf2. Classically, one would expect 
only a minimal signal in this configuration (from Poisson's 
equation, assuming diffusive transport, one estimates a signal 
of about 40 fi), and certainly not one that increases so strongly 
when the bulk of the sample is depleted. The signal measured 
here is fully nonlocal, and can be taken (as was done twenty 
years ago for the QH regime) as definite evidence of the exis- 
tence of edge channel transport in the QSH regime. 



III. THREE-DIMENSIONAL TOPOLOGICAL INSULATORS 

The model Hamiltonian for the 2D topological insulator in 
HgTe QWs also gives a basic template for generalization to 
3D, leading to a simple model Hami ltonian for a class o f mate- 
rials: Bi2Se3, Bi2Te3, and Sb2Te3 (IZhang et al.l\2QQ% . Sim- 
ilar to their 2D counterpart the HgTe QWs, these materials 
can be described by a simple but realistic model, where SOC 
drives a band inversion transition at the V point. In the topo- 
logically nontrivial phase, the bulk states are fully gapped, 
but there is a topologically protected surface state consisting 
of a single massless Dirac fermion. The 2D massless Dirac 
fermion is "helical", in the sense that the electron spin points 
perpendicularly to the momentum, forming a left-handed he- 
lical texture in momentum space. Similarly to the ID helical 
edge states, a single massless Dirac fermion state is "holo- 
graphic", in the sense that it cannot occur in a purely 2D sys- 
tem with TR symmetry, but can exist as the boundary of a 3D 
insulator. TR invariant single-particle perturbations cannot in- 
troduce a gap for the surface state. A gap can open for the 
surface state when a TR breaking perturbation is introduced 
on the surface. Moreover, the system becomes full insulating, 
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both in the bulk and on the surface. In this case, the topolog- 
ical properties of the fully gapped insulator are characterized 
by a novel topological magnetoelectric effect. 

Soon after the theoretical predict ion of the 3D topo logi- 
cal insu lator in th e Bi 2 Te^. Sb^Te^ ()Zhang etal. , 2009 ) and 



Bi 2 Se 3 dXia et al.l 120091: IZhang et al.l 120091) class of materi 



als, angle-resolved photoemission spectroscopy (AR PES) ob 



serve d the surface s t ates w i th a single Dirac c one dChen et al 



20091: iHsieh et al.l 120091 IXia et all l2009h . Furthermore, 



spin-resolved ARPES measurements indeed observed the 
left-han ded helical spin texture of the massless Dirac 
fermion (Hsieh et ^Zl l2009f) . These pioneering theoretical and 
experimental works inspired much of the subsequent develop- 
ments which we review in this section. 

We take advantage of the model simplicity of the Bi2Se3, 
Bi2Te3, Sb2Te3 class of 3D topological insulators and give 
a pedagogical introduction based on this particular material 
system. In the next section, we shall introduce the general 
theory of the topological insulators. The electronic structure 
of the Bi2Se3, Bi2Te3, Sb2Te3 class of topological insulators 
is simple enough to be captured by a simple model Hamil- 
tonian. However, more powerful methods are needed to de- 
termine the topological properties of materials with a more 
complex electronic struc ture. In this regard, the TBT has 



played an importa nt role (IFu et al.l 120071 : iMoore and Balents 



20071: iRovl l2009h. In par ticular, a method due to Fu and 



Kane (IFu and KaneL 120071) gives a simple algorithm to de- 
termine the topological properties of an arbitrarily complex 
electronic structure with inversion symmetry. This method 
predicts that Bi^Sbi-^ is a topological insulator for a certain 



range of composition x. ARPES measurements dHsieh et al 



20081) have indeed observed topologically nontrivial surface 



states in this system, giving the first example of a 3D topo- 
logical insulator. The topological properties of this material 
have bee n further investigated both theoretically and experi 



menta lly (INishide et aZl l2010l: iTeo et al.l l2008t IZhang et al.l 



2009). however, the surface states in Bi^Sbi-^ are rather 



complicated, and cannot be described by simple model Hamil- 
tonians. For this reason, we focus on the Bi2Se3, Bi2Te3, 
Sb2Te3 class of topological insulators in this section. 



A. Effective model of the three-dimensional topological 
insulator 

In this review w e focus on an effect ive model for 3D topo- 
logical insulators (IZhang et 'all I2009T) which, simply by ad- 
justing parameters, is valid for studying the properties of 
Bi2Se3, Bi2Te3, and Sb2Te3. Bi2Se3, Bi2Te3, and Sb2Te3 
share the same rhombohedral crystal structure with space 
group D^ d {Rim) and five atoms per unit cell. For exam- 



ple, the crystal structure of Bi2Se3 is shown in Fig. [T8f a), and 
consists of a layered structure where individual layers form a 
triangular lattice. The important symmetry axes are a trigonal 
axis (three-fold rotation symmetry) defined as the z axis, a bi- 
nary axis (two-fold rotation symmetry) defined as the x axis, 



Quintuple 
layer 





FIG. 18 (a) Crystal structure of Bi2Se3 with three primitive lattice 
vectors denoted by ti,2,3- A quintuple layer with Sel-Bil-Se2-Bil'- 
Sel' is indicated by the red box. (b) Top view along the z direc- 
tion. Triangular lattice in one quintuple layer has three inequiva- 
lent positions, denoted by A, B and C. (c) Side view of the quin- 
tuple layer structure. Along the z direction, Se and Bi atomic lay- 
ers are stacked in the sequence C(Sel')-A(Sel)-B(Bil)-C(Se2)- 

A(Bil')-B(Sel')-C(Sel)- • • . The Sel (Bil) layer is related to the 
Sel' (Bil') layer by inversio n,where Se2 atoms play the role of in- 
version center. Adapted from Z hang et al. 



and a bisectrix axis (in the reflection plane) defined as the y 
axis. The material consists of five-atom layers stacked along 
the z direction, and known as quintuple layers. Each quintuple 
layer consists of five atoms per unit cell with two equivalent Se 
atoms denoted by Sel and Sel' in Fig.[T87b). two equivalent 
Bi atoms denoted by Bil and Bil' in Fig. [T87 b). and a third Se 
atom denoted by Se2 in Fig.[T"87b). The coupling between two 
atomic layers within a quintuple layer is strong, while that be- 
tween quintuple layers is much weaker, and predominantly of 
the van der Waals type. The primitive lattice vectors ti 2 j3 and 
rhombohedral primitive unit cells are shown in Fig.[T8]'a). The 
Se2 site plays the role of an inversion center. Under inversion, 
Bil is mapped to Bil' and Sel is mapped to Sel'. 

To get a better understanding of the band structure and or- 
bitals involved, we start from the atomic energy levels and 
then consider the effects of crystal field splitting and SOC on 
the energy eigenvalues at the T point in momentum space. 
This is summarized schematically in three stages (I), (II) and 
(III) [Fig. [T9la)l. Since the states near the Fermi level are 
primarily from p-orbitals, we will neglect the s-orbitals and 
start from the atomic p-orbitals of Bi (electronic configuration 
6s 2 6p 3 ) and Se (4s 2 4p 4 ). In stage (I), we consider chemical 
bonding between Bi and Se atoms within a quintuple layer, 
which corresponds to the largest energy scale in this problem. 
First, we can recombine the orbitals in a single unit cell ac- 
cording to their parity. This results in three states (two odd, 
one even) from each Se p-orbital and two states (one odd, 
one even) from each Bi p-orbital. The formation of chemical 
bonds hybridizes the states on the Bi and Se atoms, and pushes 
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down all the Se states and lifts all the Bi states. In Fig. [T9l a), 
these five hybridized states are labeled as \Pl x z yz "), |P2^ yz ) 
and | P0~ y z ), where the superscripts ± stand for the parity of 
the corresponding states. In stage (II), we consider the effect 
of crystal field splitting between different p-orbitals. Accord- 
ing to the point group symmetry, the p z orbital is split from 
the p x and p y orbitals while the latter two remain degenerate. 
After this splitting, the energy levels closest to the Fermi en- 
ergy turn out to be the p z levels |P1+) and \P2~). In the last 
stage (III), we take into account the effect of SOC. The atomic 
SOC Hamiltonian is given by Hso = AL • S, with L, S the 
orbital and spin angular momentum, respectively, and A the 
strength of SOC. The SOC Hamiltonian mixes spin and or- 
bital angular momenta while preserving the total angular mo- 
mentum. This leads to a level repulsion between |Pl+,f) 
and |Pl+ +i ,4-)> and between similar combinations. Conse- 
quently, the energy of the |Pl+,f (J,)) state is pushed down 
by the effect of SOC, and the energy of the |P2~ , f (4)) state 
is pushed up. If SOC is larger than a critical value A > A c , 
the order of these two energy levels is reversed. To illustrate 
this inversion process explici tly, the energy level s |P1 +) and 
|P2j) have been calculated (IZhang et al.l 120091) for a model 
Hamiltonian of Bi2Se3 with artificially rescaled atomic SOC 
parameters A(Bi) = xAo(Bi), A(Se) = a;Ao(Se), as shown in 
Fig.QStb). Here A (Bi) = 1.25 eV and A (Se) = 0.22 eV 
are the actual values of the SOC streng th for Bi and Se atoms, 
respectively (IWittel and Mannel 1 19741) . From Fig. [T9l b). one 
can clearly see that a level crossing occurs between |Plj~) 
and |P2~) when the SOC strength is about 60% of its actual 
value. Since these two levels have opposite parity, the inver- 
sion between them drives the system into a t opological insu 
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FIG. 19 (a) Schematic picture of the evolution from the atomic 
Px, y ,z orbitals of Bi and Se into the conduction and valence bands 
of Bi2Se3 at the V point. The three different stages (I), (II) and 
(III) represent the effect of turning on chemical bonding, crystal field 
splitting, and SOC, respectively (see text). The blue dashed line rep- 
resents the Fermi energy, (b) The energy levels \Plf) and \P27) 
of Bi2Se3 at the T point versus an artificially rescaled atomic SOC 
A(Bi) = zAo(Bi) = 1.25a; [eV], A(Se) = zA (Se) = 0.22s [eV] 
(see text). A level crossin g occurs between these two states at 
~ 0.6. Adapted from Z hang et a/.Ll2009l . 



lator p hase, similar to the case of HgTe QWs dBernevig et al 
120061) . Therefore, the mechanism for the occurrence of a 3D 
topological insulating phase in this system is closely analo- 
gous to the mechanism f or the 2D QSH effect (2D topolog- 



ical insulator) in HgTe ( Bernevig et al. . 20061) . More pre 



ants (IWinklerl 120031) at a finite wavevector k. The important 
symmetries of the system are TR symmetry T, inversion sym- 
metry /, and three-fold rotation symmetry C 3 aroung the z 
axis. In the basis {|P1+ t) , |P2" t) , |Pl+4> , 1^2" 4-)}, 
the representation of these symmetry operations is given 
by T = ia v JC (g> I 2x2 , I = I2X2 ® t 3 and C 3 = 
exp (if cr z ® I2X2), where l nxn is the n x n identity ma- 
trix, 1C is the complex conjugation operator, and (j x ' y,z and 
T x ' y ' z denote the Pauli matrices in the spin and orbital space, 
respectively. By requiring these three symmetries and keeping 
only terms up to quadratic order in k, we obtain the following 
generic form of the effective Hamiltonian: 



cisely, to determine whether or not an inversion-symmetric 
crystal is a topological insulator, we must have full knowl- 
edge of the states at all of t he eight TR invariant momenta 
(TRIM) dFu and Kanel 120071) . The system is a (strong) topo- 
logical insulator if and only if the band inversion between 
states with opposite parity occurs at odd number of TRIM. 
The parity of the Bloch states at all TRIM have been studied 
by ab initio methods for the four materia ls Bi2Se3, Bi 2 Te 3 , 

' l2009h . Comparing the M-B 1 k 2 z -B 2 k 
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Sb 2 Se 3 , and Sb 2 Te 3 (IZhang et al. 



with£; ± = k x ±ik y , e (k) = C+Dik 2 z +D 2 k\ andTW(k) 



Bloch states with and without SOC, one conclude that Sb 2 Se 3 
is a trivial insulator, while the other three are topological insu- 
lators. For the three topological insulators, the band inversion 
only occurs at the Y point. 

Since the topological nature is determined by the physics 
near the Y point, it is possible to write down a simple effective 
Hamiltonian to characterize the low-energy, long-wavelength 
properties of the system. Starting from the four low-lying 
states |Pl+,t(4)) and |P2~,t (•!•)) at the Y point, such 
a Hamiltonian can be constructed by the theory of invari- 



The parameters in the effective model can 
be determined by fitting the energy spectru m of the effective 
Hamiltonian to that of ab initio calculations ( ILiu et qfl l2010l 



Zhang et al. The fitting leads to the parameters 

displayed in Tableim dLiu et al[\l0m . 

Except for the identity term eo(k), the Hamiltonian OTI ) is 
similar to the 3D Dirac model with uniaxial anisotropy along 
the z direction, but with the crucial difference that the mass 
term is k-dependent. From the fact that M, B>i,B 2 > 
we can see that the order of the bands |Tl+,f (4)) and 
|T2j,f (4)) is inverted around k = compared with large 
k, which correctly characterizes the topologically nontrivial 
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TABLE II The parameters in the model Hamilto- 
nian d31t obtai ned from fitting to ab initio calculation. 
Adapted from dLiuef a/.Ll2010h . 





Bi 2 Se s 


Bi 2 Te 3 Sb 2 Te 3 


A^eV-A) 


2.26 


0.30 


0.84 


A 2 (eV-k) 


3.33 


2.87 


3.40 


C(eV) 


-0.0083 


-0.18 


0.001 


Di(eV-A 2 ) 


5.74 


6.55 


-12.39 


D 2 {eV-k 2 ) 


30.4 


49.68 


-10.78 


M(eV) 


0.28 


0.30 


0.22 


Bi(eV-A 2 ) 


6.86 


2.79 


19.64 


B 2 {eV-A 2 ) 


44.5 


57.38 


48.51 



nature of the system. In addition, the Dirac mass M, i.e. the 
bulk insulating gap, is ~ 0.3 eV, which allows the possibil- 
ity of having a room-temperature topological insulator. Such 
an effective model can be used for further theoretical study of 
the Bi2Se3 system, as long as low-energy properties are con- 
cerned. 

Corrections to the effective Hamiltonian d3TT i that are of 
higher order in k can also be considered. To cubic (k 3 ) 
order, some new terms can break the continuous rotation 
symmetry around the z axis to a discrete three-fold rotation 
symmetry C3. Correspondingly, the Fermi surface of the 
surface state acquires a hexagonal shape jFul 120091) . which 
leads to important consequences for experiments on topo- 
logic al insulators such as surface state quasipa r ticle interfer- 



logical insulato rs such as surface s tate quasipart icle lnterter 
ence ( Alpichshev et al.i 2010t Lee et all 2009 ; Zhang et al. 
20091: Izhou et al. ' ' 



2009) 



A modified version of the effective 
model OTb taking into account corrections up to k 3 has been 
obtained for the three topological insulator s Bi2Se,3, I^Tes 



and Sb2Te3 based on ab initio calcu lations dLiu et ail 2010h . 
In this same work ( Liu et al. , 2010l) . an eight-band model is 
also proposed for a more quantitative description of this fam- 
ily of topological insulators. 



B. Surface states with a single Dirac cone 

The existence of topological surface states is one of the 
most important properties of topological insulators. The sur- 
face states can be directly extracted from ab initio calcula- 
tions by constructing maximally localized Wannier functions 
and calculating the loca l density of states on an open bound- 



ary (IZhang et all 120091) . The result for the Bi2Se3 family of 
materials is shown in Fig. l20t a)-(d). where one can clearly 
see the single Dirac-cone surface state for the three topo- 
logically nontrivial materials. However, to obtain a better 
understanding of the physical origin of topological surface 
states, it is helpful to sho w how the surface states emerge from 
the effec t ive model d3TT> dLinder et all 120091 iLiu et all l2010t 



Lu et al. , 2010t Zhang et al. , 2009 ). The surface states can be 
obtained in a similar way as the edge states of the BHZ model 
(Sec.HLUi. 



Consider the model Hamiltonian (l3TT l on the half-space z > 
0. In the same way as in the 2D case, we can divide the model 
Hamiltonian into two parts, 



H = Hn 



H = i{k z ) + 



Hi 



D 2 k\ 



( M{k z ) 
Aik z 


v 
/ -B 2 k\ 



V A 2 k + 



Aik z 
-M(k z ) 





B 2 k\ 
A 2 k+ 







M(k z . 
-Aik ; 



-B 2 k\ 




(32) 

\ 



-A x k z 
-M(k z ) J 

A 2 k_ \ 








, (33) 



with e(k z ) = C + D x k 2 z and M(k z ) = M - B x k 2 z . H 
in Eq. (0 and Eq. d32l are identical, with the parameters 
A,B,C,D,M in Eq. ® replaced by A 1 ,B 1 ,C,D 1 ,M in 
Eq. d32t . Therefore, the surface state at k x = k y = is deter- 
mined by the same equation as that for the QSH edge states. A 
surface state solution exists for M/B\ > 0. In the same way 
as in the 2D case, the surface state has a helicity determined 
by the sign of A\/B\. (Here and below we always consider 
the case with B X B 2 > 0, A X A 2 > 0.) 

In analogy to the 2D QSH case, the surface effective model 
can be obtained by projecting the bulk Hamiltonian onto 
the surface states. To the leading order in k x ,k y , the ef- 
fectiv e surface Hamiltonian H SUT { has th e following matrix 
form dLiu et al.[ l2010t IZhang et all l2009k : 



H snli (k x ,k y ) = C + A 2 (a x k y - a v k x ) 



(34) 



Higher order terms such as k 3 terms break the axial symme- 
try around the z axis down to a three-fold ro tation symmetry . 



which has been studied in the literature (|Fu, 2009; Li u et al. 
2010h . For A 2 = 4.1eV-A, he velocity of the surface states is 
given by v = A 2 /h~ 6.2x 10 5 m/s, which agrees reasonably 
with ab initio results [Fig. [20) v ~ 5.0 x 10 5 m/s. 

To understand the physical properties of the surface states, 
we need to analyze the form of the spin operators in this 
system. By using the wave function from ab initio calcu- 
lations and projecting the spin operators onto the subspace 
spanned by the four basis states, we obtain the spin operators 
for our model Hamiltonian, with matrix elements between sur- 
face states given by (^/ a \S x \^p) = S x0 a^, (^ a \S y \^p) = 
Syoaf and (^ a \S z \^fp) = S z0 <jf, with S x(ViZ)0 some pos- 
itive constants. Therefore, we see that the Pauli er matrix 
in the model Hamiltonian ( f34b is proportional to the physi- 
cal spin. As discussed above, the spin direction is determined 
by the sign of the parameter A1/B1, which depends on mate- 
rial properties such as the atomic SOC. In the Bi2Se3 family 
of materials, the upper Dirac cone has a left-handed helicity 
when looking from above the surface [Fig.l20l'e).(f)l. 

From the discussion above, we see that the surface state 
is described by a 2D massless Dirac Hamiltonian ( f34b . An- 
other well-known system with a similar property i s graphene, 
a single sheet of graphite ( Castro Neto et al. , 2009b . However, 
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(f) 




FIG. 20 (a)-(d) Energy and momentum dependence of the local den- 
sity of states for the Bi2Se3 family of materials on the [111] surface. 
A warmer color represents a higher local density of states. Red re- 
gions indicate bulk energy bands and blue regions indicate a bulk 
energy gap. The surface states can be clearly seen around V point 
as red lines dispersing inside the bulk gap. (e) Spin polarization of 
the surface states on the top surface, where the z direction is the sur- 
fac e normal, pointi ng outwards. Adapted from IZhang et all |2009| 
and lLiu et al. 



pomtn 



there is a key difference between the surface state theory for 
3D topological insulators and graphene or any 2D Dirac sys- 
tem, which is the number of Dirac cones. Graphene has four 
Dirac cones at low energies, due to spin and valley degener- 
acy. The valley degeneracy occurs because the Dirac cones 
are not in the vicinity of k = but rather near the two Bril- 
louin zone corners K and K. This is generic for a purely 2D 
system: only an even number of Dirac cones can exist in a TR 
invariant system. In other words, a single 2D Dirac cone with- 
out TR symmetry breaking can only exist on the surface of a 
topological insulator, which is also an alternative way to un- 
derstand its topological robustness. As long as TR symmetry 
is preserved, the surface state cannot be gapped out because 
no purely 2D system can provide a single Dirac cone. Such a 
surface state is a "holographic metal" which is 2D but deter- 
mined by the 3D bulk topological property. 

In this section we discussed the surface states of an insula- 
tor surrounded by vacuum. This formalism can be straightfor- 



wardly gejiendjzedto tf^int^^ two i nsu- 

lators dFradkin et all 1 1986t [Volkov and Pankratovi 1 1985b . In 
these pioneering works, the interface states between PbTe and 
SnTe were investigated. The interface states consist of four 
Dirac cones. Therefore, they are topologically trivial and not 
generally stable under TR invariant perturbations. The surface 
states of topological insulators are als o similar to th e domain 
wall fermions of lattice gauge theory dKaplanll 19921) . In fact, 
domain wall fermions are precisely introduced to avoid the 
fermion doubling problem on the lattice, which is similar to 
the concept of a single Dirac cone on the surface of a topolog- 
ical insulator. 

The helical spin texture described by the single Dirac cone 
equation fl3"4t leads to a general relation between charge cur- 
rent density j (x) and spin density S(x) on the surface of the 
topological insulator ( Raghu etall 2010b : 



j(x) = w[i/) t (x)crV'(x) x z] = t>S(x) x z. (35) 

In partic ular, the plasmon mode on the surface generally c ar- 
ries spin ( Burkov and Hawthorn , 20ld : Raghu et al. , 20101) . 



C. Crossover from three dimensions to two dimensions 

From the discussion above, one can see that the models 
describing 2D and 2D topological insulators are quite simi- 
lar. Both systems are described by lattice Dirac-type Hamil- 
tonians. In particular, when inversion symmetry is present, 
the topologically nontrivial phase in both models is charac- 
terized by a band inversion between two states of opposite 
parity. Therefore, it is natural to study the relation between 
these two topological states of matter. One natural question 
is whether a thin film of 3D topological insulator, viewed as 
a 2D system, is a trivial insulator or a QSH insulator. Be- 
sides theoretical interest, this problem is also relevant to ex- 
periments, especially in the Bi2Se3 family of materials. In- 
deed, these materials are l ayered and can be easily grown as 



thin f ilms either by MBE dLi et all [2oTol 



20091) . catalyst-free vapor-so lid growth dKong et q/.L l2010h 



or by mechanical exfoliation (|Hqng et al. 



2009; Zhang et al. 



To: 



2010; Shahil et al. 



201& iTeweldebrhan et Several theoretical works 



studie d thin films of the Bi2Se,3 family of topological insu - 
lators dLinder et ali 120091: iLiu et ali 1201 0t iLu et ali |2010|) . 
Interestingly, thin fil ms of proper thicknesses are pred icted to 
form a QSH insulator dLiu et ILu et which 

may constitute an approach for simpler realizations of the 2D 
QSH effect. 

Such a crossover from 3D to 2D topological insulators can 
be studied from two points of view, either from the bulk states 
of the 3D topological insulator or from the surface states. 
We first consider the bulk states. A thin film of 3D topo- 
logical insulator is described by restricting the bulk model 
( T3Tb to a QW with thickness d, outside which there is an in- 
finite barrier describing the vacuum. To establish the con- 
nection between the 2D BHZ model (O and the 3D topo- 
logical insulator model OTb . we start from the special case 
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A\ = and consider a finite Ai later on. For A\ = and Therefore, the complete effective Hamiltonian is given by 
^2 = k y = 0, the Hamiltonian PIT ) becomes diagonal and the 
Schrodinger equation for the infinite QW can be easily solved. 
The Hamiltonian eigenstates are simply given by \E n (H n )) = 

v /| s i n («|£ + ihl) |A), with |A) = |P1+, t (I)) for elec- 
tron subbands and | A) = |P2~,f (I)) for hole subbands. The 
corresponding energy spectrum is E e (n) = C + M + (£>i — 



H SUI f(k x , ky) — A.2 
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o \ 




—ik + 
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o / 



BO and £7 h (n) = C- M+ {D 1 + B 1 ){^f)\ respec- 

tively. We assume M < and Pi < so that the system 
stays in the inverted regime. The energy spectrum is shown in 
Fig. l2TTa). When the width d is small enough, electron sub- 
bands E n have a higher energy than the hole subbands H n 
due to quantum confinement effects. Because the bulk bands 
are inverted at the V point (M < 0), the energy of the elec- 
tron subbands will decrease with increasing d towards their 
bulk value M < 0, while the energy of the hole subbands 
will increase towards — M > 0. Therefore, there must exist a 
crossing point between the electron and hole subbands. 

When a finite Ai is turned on, the electron and hole bands 
are hybridized so that some of the crossings between the QW 
levels are avoided. However, as shown in Fig. l2TT b). some 
level crossings cannot be lifted, which is a consequence of in- 
version symmetry. When the band index n is increased, the 
parity of the wave functions alternates for both electron and 
hole subbands. Moreover, the atomic orbitals forming elec- 
tron and hole bands are |Pl+,t (I)) and |P2~,t (I)) re- 
spectively, which have opposite parity. Consequently, \E n ) 
and \H n ) with the same index n have opposite parity, so that 
their crossing cannot be avoided by the A\ term. When a 
finite k x , k y is considered, each level becomes a QW sub- 
band. The bottom of the lowest conduction band and the top 
of the highest valence band are indicated by Sf and S 2 in 
Fig. l2TTb). Since these two bands have opposite parity, each 
level crossing between them is a topological phase transition 
between trivial and QSH insulator phases (iBernevig et ail 
20061 iFu and KaneL [2007). Since the system must be triv- 
ial in the limit d — >• 0, we know that the first QSH insulator 
phase occurs between the first and second level crossing. In 
the A\ — >• limit, the crossing positions are given by the crit- 
ical well thicknesses d cn = nn^J-^jj. In principle, there is 

an infinite number of QSH phases between e? C! 2n-i and d C: 2n- 
However, as seen in Fig. l2TT b). the gap between \S^) and 
|£^~) decays quickly for large d. In the 3D limit d —> oo, the 
two states become degenerate and actually form the top and 
bottom surface states of the bulk crystal [Fig. l2TT c)l. 

This relation between the QW valence and conduction 
bands and the surface states in the d — > oo limit suggests an 
alternative way to understand the crossover from 3D to 2D, 
i.e. from the surface states. In the 3D limit the two surfaces 
are decoupled and are the only low-energy states. The top sur- 
face is described by the effective Hamiltonian (f34-b while the 
bottom surface is obtained from the top surface by inversion. 



When a slab of finite thickness is considered, the two surface 
states start overlapping, such that off-diagonal terms are intro- 
duced in the effective Hamiltonian. An effective Hamiltonian 
consistent with inversion and TR symmetry and incorporating 
inter-surface tunneling is given by 



H smf (k x ,k y ) = A 2 
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,(36) 



where M2D is a TR invariant mass term due to inter-surface 
tunneling, which generally depends on the in-plane momen- 
tum. Equation (l36l l is unitarily equivalent to the BHZ Hamil- 
tonian (O for HgTe QWs. Whether the Hamiltonian corre- 
sponds to a trivial or QSH insulator cannot be determined 
without studying the behavior of this model at large momenta. 
Indeed, we are missing a regularization term which would 
play the role of the quadratic term Bk 2 in the BHZ model 
(Sec. III. Al l. However, the transitions between trivial and non- 
trivial phases are accompanied by a sign change in M2D, in- 
dependently of the regularization scheme at high momenta. 
Upon variation of d, the sign of the inter-surface coupling 
M 2 d oscillates because the surface state wave functions os- 
cillate [Fig. l2TT b)l. Therefore, we reach the same conclusions 
as in the bulk approach. 

The results above obtained from calculations using an ef- 
fective model are also confirmed by first-principle calcula- 
tions. The parity eigenvalues of occupied bands have been 
calculated as a fu nction of the thi ckness of the 3D topologi- 
cal insulator film ( ILiu et al. ■ l2010h . from which the topologi- 
cal nature of the film can be inferred. The result is shown in 
Fig-Eld), which confirms the oscillations found in the effec- 
tive model. The first nontrivial phase appears at a thickness of 
three quintuple layers, i.e. about 3 nm for Bi2Se3. 



D. Electromagnetic properties 

In previous subsections, we have reviewed bulk and sur- 
face properties of 3D topological insulators, as well as their 
relation to 2D topological insulators (QSH insulators), based 
on a microscopic model. From the effective model of surface 
states, one can understand their robustness protected by TR 
symmetry. However, similarly to the quantized Hall response 
in QH systems, the topological structure in topological insu- 
lators should not only lead to robust gapless surface states, 
but also to unique, quantized electromagnetic response coef- 
ficients. The quantized electromagnetic response of 3D topo- 
logical insulators turns out to be a topological magnetoelectric 
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FIG. 21 Energy levels versus quantum well thickness for (a) Ai — 
eV-A, (b) Ai = 1.1 eV-A. Other parameters are taken from 
IZhang efa/.Ll2009l . Shaded regions indicate the QSH regime. The 
blue dashed line in (b) shows how the crossing between \Ei(Hi)) 
and l-f/2 (£^2)) evolves into an anti-crossing when Ai 7^ 0. (c) 
Probability density in the state ) (same for |S^~)) for Ai = 
1.1 eV-A and d = 20 nm. (d) Band gap and total parity from ab 
initio calculations on Bi2Se3, plotted a s a function of the number of 
quintuple layers. From lLiu et a/.Ll2blO] . 



effect (TME) dOi et all [20081 120091) . which occurs when TR 
symmetry is broken on the surface, but not in the bulk. The 
TME effect is a generic property of 3D topological insula- 
tors, which can be obtained theoretically fr om generic models 
and from an effective field theory ap proach ( lEssin et al. . 120091: 



Fu and Kanel l2007t lOi et ali [2008), independently of micro 



scopic details. However, in order to develop a physical intu- 
ition for the TME effect, in the section we review this effect 
and its physical consequences based on the simplest surface 
effective model, and postpone a discussion in the framework 
of a general effective theory to Sec. [IV] We shall also discuss 
various experimental manifestations of the TME effect. 



1 . Half quantum Hall effect on the surface 

We start by analyzing generic perturbations to the 
effective surface state Hamiltonian fl34l ). The only 
momentum-independent perturbation one can add is Hi = 
J2a=x y z m o ff °' an d tne perturbed Hamiltonian has the 

spectrum = ±\J (A 2 k y + m x ) 2 + (A 2 k x — m y ) 2 + ml. 
Thus, the only parameter that can open a gap and destabi- 
lize the surface states is m z , and we will only consider this 
perturbation in the following. The mass term m z a z is odd 
under TR, as expected from the topological stability of sur- 
face states protected by TR symmetry. By comparison, if the 
surface states consist of even number of Dirac cones, one can 
check that a TR invariant mass term is indeed possible. For 
example, if there are two identical Dirac cones, an imaginary 



coupling between the cones can be introduced, which leads to 
the gapped TR invariant Hamiltonian 



A 2 {<J x k y — a v k x ) —ima z 

ima z A 2 (a x k y — a y k x ) 



From such a difference between an even and an odd number 
of Dirac cones, one sees that the stability of the surface theory 
( l34l ) is protected by a Z 2 topological invariant. 

Although the surface state with a single Dirac cone does 
not remain gapless when a TR breaking mass term m z a z is 
added, an important physical property is induced by such a 
mass term: a half-integer quantized Hall conductance. As 
discussed in Sec. III. El the Hall conductance of a generic 
two-band Hamiltonian /i(k) = d a (k)o- a is determined by 
Eq. J291 l. which is the winding number of the unit vector 
d(k) = d(k)/|d(k)| on the Brillouin zone. The perturbed 
surface state Hamiltonian 



Hsmf (k) = A 2 a x k y - A 2 a y k x + m z a z 



(37) 



corresponds to a vector d(k) = (A 2 k y , — A 2 k x , m z ). At 
k = 0, the unit vector d(k) = (0, 0, m z /\m z \) points to- 
wards the north (south) pole of the unit sphere for m z > 
(m z < 0). For |k| \m z \/A 2 , the unit vector d(k) ~ 
A 2 {k y , —k x , 0)/|k| almost lies in the equatorial plane of the 
unit sphere. From such a "meron" configuration one sees that 
d(k) covers half of the unit sphere, which leads to a winding 
number ±1/2 and corresponds to a Hall conductance 



OH 



(38) 



From this formula, it can be seen that the Hall conductance 
remains finite even in the limit m z — > 0, and has a jump at 
m z = 0. As a property of the massive Dirac model, such a 
half Hall conductance has been studied a long time ago in high 
energy physic s. In that context, the effect is ter med the "par- 
ity anomaly" dRedlichl Il984t ISemenoffl Il984l) . because the 
massless theory preserves parity (and TR) but an infinitesimal 
mass term necessarily breaks these symmetries. 

The analysis above only applies if the continuum effec- 
tive model ( [34-b applies, i.e. if the characteristic momentum 
\ m z\/A 2 is much smaller than the size of the Brillouin zone 
2n/a with a the lattice constant. Since deviations from this 
Dirac-type effective model at large momenta is not included in 
the above calcul ation of the Hall conductance dFu and Kane , 
20071: ILeelE009l) . one cannot unambiguously predict the Hall 
conductance of the surface. In fact, if the effective theory de- 
scribes a 2D system rather than the surface of a 3D system, 
additional contributions from large-momentum corrections to 
the effective model are necessary, since the Hall condutance 
of any gapped 2D band insula t or mu st be quantized in integer 
units of e 2 /h dThouless et al. . 1982 ). For example, the QAH 
insulator [Eq. d28l l with mass term AI — > is also described 
by the same effective the ory as ([37|. bu t has H all conductance 
or 1 rather than ±1/2 dFradkin et a/.[|l986l) . 
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Interestingly, the surface of a 3D topological insulator is 
different from all 2D insulators, in the sense that such contri- 
butions from la rge momenta va nish due to the requirement of 
TR symmetry JOi et all [2008). This fact is discussed more 
rigorously in Sec.|IV]based on the general effective field the- 
ory. Here we present an argument based on the bulk to surface 
relationship. To understand this, consider the jump in Hall 
conductance at m z = 0. Although deviations from the Dirac 
effective model at large momenta may lead to corrections to 
the Hall conductance for a given m z , the change in Hall con- 
ductance Act h = <r H {m z -> + )-<T H (m z ->■ 0") = 
is independent of the large-momentum contributions. Indeed, 
the effect of the mass term m z a z on the large-momentum sec- 
tor of the theory is negligible as long as m z — > 0. Therefore, 
any contributions to ct# from large momenta should be con- 
tinuous functions of m z , and thus cannot affect the value of 
the discontinuity Act# . On the other hand, since the surface 
theory with m z = is TR invariant, TR transforms the system 
with mass m z to that with mass — m z . Consequently, from TR 
symmetry we have 

<?H(m z -> 0+) = -<r H (m z -> CT). 

2 

Together with the condition Actw = r^T-, we see that the 

& " \ni z | h ' 

half Hall conductance given by Eq. <f38b is robust, and the con- 
tribution from large-momentum corrections must vanish. By 
comparison, in a 2D QAH Hamiltonian discussed in Sec. III.EI 
i.e. the upper 2x2 block of Eq. (0, the Hamiltonian with 
mass M is not the TR conjugate of that with mass —M, and 
the above argument does not apply. Therefore, the half Hall 
conductance is a unique property of the surface states of 3D 
topological insulators which is determined by the bulk topol- 
ogy. This property distinguishes the surface states of 3D topo- 
logical insulators from all pure 2D systems, or topologically 
trivial surface states. 

The analysis above has only considered translationally in- 
variant perturbations to the surface states, but the conclu- 
sions remain robust when disorder is considered. A 2D 
metal without SOC belongs to the orthogonal or unitary 
symmetry classes of random Hamiltonians, the eigenfunc- 
tions of which are always localized when random disorder 
is in troduced. This effect is known as Anderson localiza- 
tion ( Abrahams et al. , 19791) . Anderson localization is a quan- 
tum interference effect induced by constructive interference 
between different backscattering paths. By comparison, a sys- 
tem with TR invariance and SOC belongs to the symplec- 
tic class, where the constructive interference becomes de- 
structive. In that case, the system has a metallic phase at 
weak disorder, which turns into an insulator phase by go- 
ing through a metal-insulator t ransition at a certai n disor - 



der strength fevers and Mirlini 120081: iHikami et all 1 1980b 



Naively, one would expect the surface state with nonmagnetic 
disorder to be in the symplectic class. However, Nomura, 
Koshino, and Ryu showed that the su rface state is metalli c 



On the contrary, with TR symmetry breaking disorder, 
the system belongs to the unitary class, which exhibits lo- 
calization for arbitrarily weak disorder strength. While 
the longitudinal resistivity flows to infinity due to local- 
ization, t he Hall conductivity flows to the quantized value 
±e 2 /2h dNomura et al. , 2008). Therefore, the system enters 
a half QH phase once an infinitesimal TR symmetry break- 
ing perturbation is introduced, independently of the detailed 
form of the TR breaking perturbation. Physically, TR break- 
ing disorder is induced by magnetic impurities, the spin of 
which not only contributes a random TR breaking field, but 
also has its own dynamics. For example, the simplest ex- 
change interaction between impurity spin and surface state 
can be written as H{ at = J^i Ji^i ■ ip'<rip(Ri) with S, the 
impurity spin, ip<(Tip(Ri) the spin density of surface elec- 
trons at the impurity position Rj, and Ji the exchange cou- 
pling. To understand the physical properties of the topologi- 
cal insulator surface in the presence of magnetic impurities, it 
is instructive to study the interacti on between impu rity spins 



mediated by the surface electrons dLiu et ail 120091) . As in a 
usual Fermi liquid, if the surface state has a finite Fermi wave 
vector fop, a Ruderman-Kittel-Kasuya-Yosida (RKKY) inter- 
action between the impurity spins is introd uced, the sign of 
which oscillate s with wave length oc l/2fcF (lLiu et al. . 20091: 
Ye et al. , 2010h . If the Fermi level is close to the Dirac point, 
i.e. kp 0, the sign of the RKKY interaction does not 
oscillate but is uniform. The sign of the resulting uniform 
spin-spin interaction is determined by the coupling to the sur- 
face electrons, which turns out to be ferromagnetic. Phys- 
ically, the interaction is ferromagnetic rather than antiferro- 
magnetic, because a uniform spin polarization can maximize 
the gap opened on the surface, which is energetically favor- 
able. Due to this ferromagnetic spin-spin interaction, the sys- 
tem can order ferromag netically when th e chemical potential 



is near the Dirac point dLiu et all 120091) . This mechanism is 



of great practical importance, because it provides a way to 
generate a surface TR symmetry breaking field by coating the 



tential near the Dirac point dCha et al. 
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even for an arbitrary impurity strength dNomura et all 120070 
which is consistent with the topological robustness of the sur- 
face state. 



2. Topological magnetoelectric effect 

As discussed above, the surface half QH effect is a unique 
property of a TR symmetry breaking surface, and is deter- 
mined by the bulk topology, independently of details of the 
surface TR symmetry breaking perturbation. A key differ- 
ence between the surface half QH effect and the usual inte- 
ger QH effect is that the former cannot be measured by a dc 
transport experiment. An integer QH system has chiral edge 
states which contribute to the quantized Hall current while be- 
ing connected to leads. However, as one can easily convince 
oneself, it is a simple mathematical fact that the surface of a 
finite sample of 3D topological insulator is always a closed 
manifold without an edge. If the whole surface of a topolog- 



ical insulator sample is gapped by magnetic impurities, there 
are no edge states to carry a dc transport current. If the mag- 
netic impurities form a ferromagnetic phase and there is a do- 
main wall in the magnetic moment, the Hall conductance has 
a jump at the domain wall due to the formula (l38l l. In this 
case, the jump of Hall conductance is e 2 /h across the domain 
wall, so that a chiral gapl ess edge state pr opagates along the 
domain wall [Fig. f22V a~)l ( Oi etai . 2008b . This mechanism 
provides another route towards the QAH effect without any 
external magnetic field and the associated LLs. This is very 
much alike the boundary between two ordinary QH states with 
Hall conductance ne 2 /h and (n + l)e 2 /h. The wave func- 
tion for an edge state along a straight domain wall can also be 
solved for analytically following the same procedure as that 
used in Sec. III.BI Interestingly, if one attaches voltage and 
current leads to the domain wall in the same way as for an 
ordinary Hall bar, one should observe a Hall conductance of 
e 2 /h rather than e 2 /2h, since the domain wall chiral state be- 
haves in the same way as the edge state of a ajj — e 2 /h QH 
system. Thus again we see that from dc transport measure- 
ments, one cannot observe the half Hall conductance. 

Such a difference between integer QH effect and surface 
half QH effect indicates that the surface half QH effect is 
actually a new topological phenomenon which, in terms of 
its observable consequences, is qualitatively different from 
the usual integer QH effect. Alternatively, the proper de- 
tection of this new topological phenomenon actually probes 
a uniq ue electromagnetic response prop erty of the bulk, the 
TME dEssin et all l2009t lOi et al.i 120081) . A magnetoelectric 
effect is defined as a magnetization induced by an electric 
field, or alternatively, a charge polarization induced by a mag- 
netic field. To understand the relation between surface half 
QH effect and magnetoelectric effect, consider the configu- 
ration shown in Fig. 1221 b). where the side surface of a 3D 
topological insulator is covered by magnetic impurities with 
ferromagnetic order, so that the surface is gapped and exhibits 
a half quantized Hall conductance. When an electric field E 
is applied parallel to the surface, a Hall current j is induced 
[Eq. (|39lll, which circulates along the surface. This surface 
current perpendicular to E will then induce a magnetic field 
parallel to E, so that the system exhibits a magnetoelectric 
response. The Hall response equation is written as 



(a) 



E M, 



m e . _ 
- — - — n x E, 

to 2h 



(39) 



with n a unit vector normal to the surface, and the sign of 
the mass to/|to| is determined by the direction of the surface 
magnetization. Such a Hall response is equivalent to a mag- 
netization proportional to the electric field: 

m e 2 ^ 
M, = - — E. 

\m\ 2hc 

This magnetization is a topological response to the electric 
field, and is independent of the details of the system. Simi- 
larly, a topological contribution to the charge polarization can 
be induced by a magnetic field. The complete electromagnetic 
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FIG. 22 (a) Ferromagnetic layer on the surface of topological insu- 
lator with a magnetic domain wall, along which a chiral edge state 
propagates, (b) Relation between surface half QH effect and bulk 
topological magnetoelectric effect. A magnetizatio n is induced b 
an electric field due to the surface Hall current. From lOi et fl/.Ll200i 



response of the system is described by the following modified 
constituent equations, 



H = B 4ttM + 2P 3 aE, 
D = E + 4ttP - 2P 3 aB, 



(40) 



with a = e 2 /He the fine structure constant, and P 3 = 
to/2|to| = ±1/2 the quantum of Hall conductance. A de- 
tailed explanation of the coefficient P 3 and the effective field 
theory description of the TME effect is discussed in Sec. [TV] 
In this section we focus on the physical consequences of the 
TME effect. We simply note the fact that more generally, for 
topological insulators P 3 can take the value n + 1/2 with arbi- 
trary integer n, since the number of Dirac cones on the surface 
can be any odd integer. 



3. Image magnetic monopole effect 

One of the most direct consequenc es of the TME e ffect is 

~ 2009h . 



the image magnetic monopole effect (lOi et ail 120091) . Con- 
sider bringing an electric charge to the proximity of an ordi- 
nary 3D insulator. The electric charge will polarize the dielec- 
tric, which can be described by the appearance of an image 
electric charge inside the insulator. If the same thing is done 
with a topological insulator, in addition to the image electric 
charge an image magnetic monopole will also appear inside 
the insulator. 

This image magnetic monopole effect can be studied 
straightforwardly by solving Maxwell's equations with the 
modified constituent equations d40l i. in the same way as the 
image charge problem in an ordinary insulator. Consider the 
geometry shown in Fig. 1237 a). The lower half-space z < 
is occupied by a topological insulator with dielectric constant 
£2 and magnetic permeability [i<x, while the upper half-space 
z > is occupied by a conventional insulator with dielec- 
tric constant t\ and magnetic permeability An electric 
point charge q is located at (0, 0, d) with d > 0. We assume 
that the surface states are gapped by some local TR symmetry 
breaking field to, so that the surface half QH effect and TME 
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exist. The boundary of the topological insulator acts as a do- 
main wall where P 3 jumps from 1/2 to 0. In this semi-infinite 
geometry, an image point magnetic monopole with flux g 2 is 
located at the mirror position (0, 0, — d), together with an im- 
age electric point charge q 2 . Physically, the magnetic field of 
such an image monopole configuration is induced by circu- 
lating Hall currents on the surface, which are induced by the 
electric field of the external charg e. A similar effect has bee n 
studied in the integer QH effect (IHaldane and Chenl 119831) . 
Conversely, the electromagnetic field strength inside the topo- 
logical insulator is described by an image magnetic monopole 
gi and electric charge q% in the upper half-space, at the same 
point as the external charge. The image magnetic monopole 
flux and image electric charge (qi,gi) and (92,52) are given 
by 



0.1 = <?2 = 



'92 



1 (6 1 -6 2 )(l//ii + l/^2)-4a 2 P 3 2 
ei (ei + e 2 )(l/Mi + I/M2) + 4a 2 P 2 9 ' 
4aP 3 

= ~(e 1 +e 2 )(l/^i + l/M2) + 4a 2 P 2 



9.(41) 



Interestingly, by making use of electric-magnetic duali ty these 
expre ssions can be simplified to more compact forms ( IKarch , 
20091) . 

Moreover, interesting phenomena appear when we consider 
the dynamics of the external charge. For example, consider a 
2D electron gas at a distance d above the surface of the 3D 
topological insulator. If the motion of the electron is slow 
enough (with respect to the time scale h/m corresponding to 
the TR symmetry breaking gap m), the image monopole will 
follow the electron adiabatically, such that the electron form s 
an electron-monopole composite, i.e. a dyon (lWittenl . il 9791) . 
When two electrons wind around each other, each electron 
perceives the magnetic flux of the image monopole attached 
to the other electron, which leads to statistical transmutation. 
The statistical angle is determined by the electron charge and 
image monopole flux as 



gig 
2hc 



2a 2 P 3 



(ei+e2)(l/Mi + l//i2) + 4a 2 P 2 ' 



(42) 



The image monopole can be detected directly by local 
probes sensitive to small magnetic fields, such as scan- 
ning superconducting quantum interference devices (scan- 
ning SQUID) and scanning mag netic force microscopy 



2009) 



The current due to 



(scanning MFM) dOi et al. 

the image monopole can also be detected in princi- 



ple ( Zang and Nagaosa , 2010l) 



4. Topological Kerr and Faraday rotation 

Another way to detect the TME effect is through the 
transmission and reflection of polarized light. When lin- 
early polarized light propagates through a medium which 
breaks TR symmetry, the plane of polarization of the trans- 
mitted ligh t may be rotated, which is known as the Fara- 
day effect (ILandau and Lifshitzl 119841) . A similar rotation 
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FIG. 23 (a) Image electric charge and image magnetic monopole 
due to an external electric point charge. The lower half-space is 
occupied by a topological insulator (TI) with dielectric constant ei 
and magnetic permeability [Li . The upper half-space is occupied by 
a topologically trivial insulator (e.g. vacuum) with dielectric con- 
stant ei and magnetic permeability An electric point charge q 
is located at (0,0, d). Seen from the lower half-space, the image 
electric charge q\ and magnetic monopole 171 are at (0, 0, d). Seen 
from the upper half-space, the image electric charge qi and mag- 
netic monopole (72 are at (0, 0, — d). The red (blue) solid lines repre- 
sent the electric (magnetic) field lines. The inset is a top-down view 
showing the in-plane component of the electric field on the surface 
(red arrows) and the circulating surface current (black circles), (b) Il- 
lustration of the fractional statistics induced by the image monopole 
effect. Each electron forms a "dyon" with its image monopole. When 
two electrons are exchanged, a Aharonov-Bohm phase factor is ob- 
tained, which is deter- mined by half of the image monopole flux, 
independently of the exchang e path, leading to the phenomenon of 
statistical transmutation. From lOi et a/.Ll2009l . 



may occur for light reflected by a TR symmetry breaking 
surfa ce, which is known as th e magneto-optical Kerr ef- 
fect ( ILandau and Lifshitzl 1 1 9840 . Since the bulk of the topo- 
logical insulator is TR invariant, no Faraday rotation will oc- 
cur in the bulk. However, if TR symmetry is broken on the 
surface, the TME effect occurs and a unique kind of Kerr and 
Faraday rotation is induced on the surface. Physically, the 
plane of polarization of the transmitted and reflected light is 
rotated because the electric field Eo(r, t) of linearly polarized 
light generates a magnetic field B(r, t) in the same direction, 
due to the TME effect. Similarly as for the image monopole 
effect, the Faraday and Kerr rotation angles can be calculated 
by solving Maxwell's equations with the modified constituent 
equations <l40b - In the simplest case of a single surface be- 
tween a trivial insulator and a semi-infinite topological insu- 
lator [Fig. l24la)l. the rotation angle for light incident from 
the tr i vial insulator is given by ( Karchl 20091: M acieiko et al. 



20 1 & lOi et all 120081; fTse and MacDonaldl |20 1 d) 



tan 6 k 



tan Of = 



4aP 3V / e 1 //i 1 



e 2 //U2 - ei/Mi +4a 2 P 3 2 ' 
2«P 3 



(43) 
(44) 



where ei, \x\ are the dielectric constant and magnetic perme- 
ability of the trivial insulator, and e 2 , fi 2 are those of the topo- 
logical insulator. 
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so that the monopole carries an electric charge dWittenll 1979b 
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FIG. 24 (a) Illustration of the Faraday rotation on one surface of 
topological insulator, (b) A more complicated geometry with Kerr 
and Faraday rotation on two surfaces of topological insulator. In this 
geometry, the effect of non-universal properties of the material can be 
eliminated and the quant ized magnetoele ctri c coefficient otP^ can be 
directly measured. From lOi et a/.ll2008l and lMacieiko et q/.Ll2O10l . 



Although the surface Faraday and Kerr rotations are in- 
duced by the topological property of the bulk, and are deter- 
mined by the magnetoelectric response with quantized coeffi- 
cient aP 3 , the rotation angle is not universal and depends on 
the material parameters e and [i. The TME response always 
coexists with the ordinary electromagnetic response, which 
makes it difficult to observe the topological quantization phe- 
nomenon. However, recently new proposals have been made 
to a void the dependence on non-universal material pa rame- 
ters dMaciejko et all l201fj iTse and MacDonaldi |2010b . The 



key idea is to consider a slab of topological insulator of fi- 
nite thickness with two surfaces, with vacuum on one side 
and a substrate on the other [Fig. [24l" b)l. The combination 
of Kerr and Faraday angles measured at reflectivity minima 
provides en ough information to de termine the quantized coef- 
ficient aP 3 dMacieiko et all 1201(1 



cot Op + cot 9k 
1 + cot 2 9 F 



= 2aP, 



(45) 



provided that both top and bottom surfaces have the same sur- 
face Hall conductance an = P 3 e 2 /h. In Eq. d45l ). the quan- 
tized coefficient aP 3 is expressed solely in terms of the mea- 
surable Kerr and Faraday angles. This enables a direct exper- 
imental measurement of P3 without a separate measurement 
of the non-universal optical constants e, /x of the topological 
insulator film and the substrate. If the two surfaces have dif- 
ferent surface Hall conductances, it is still possible to deter- 
mine them separately through a meas urement of the Kerr an d 
Faraday angles at reflectivity maxima dMaciejko et qfl l2010h . 



5. Related effects 

The TME has other interesting consequences. The TME ef- 
fect corresponds to a term 8~E ■ B in the action (see Sec. lIVI for 
more details), which mediate s the tran s mutat ion between elec- 
tric field and magnetic field ( Oi et a/.ll2008h . In the presence 
of a magnetic monopole, such an electric-magnetic transmu- 
tation induces an electric field around the magnetic monopole, 



q = ej^-^ with 4>o = hc/e the flux quanta. Such a com- 
posite particle carries both magnetic flux and electric charge, 
and is called a dyon. In principle, the topological insula- 
tor provides a physical system which can detect a magneti c 
monopole through this effect ( Rosenberg and Franz , 2010l) . 
For a topological insulator, we have 8 = tt which corre- 
sponds to a half charge q = e/2 for a monopole with unit 
flux. Such a half charge corresponds to a zero energy bound 
state induced by the monopole. The charge of the monopole 
is e/2 when the bound state is occupied, and —e/2 when it 
is unoccupied. Such a half charge and z ero mode is sim ilar 



to charge fractionalization in ID systems dSu et a/.l 119791) . If 
the monopole passes through a hole in the topological insula- 
tor, the charge will follow it, which co rresponds to a charge 



pumping effect ( Rosenberg et al. , 2010l) 



All effects discussed up to this point are consequences of 
the TME in a topological insulator with surface TR symmetry 
breaking. No effects of electron correlation have been taken 
into account. When the electron-electron interaction is con- 
sidered, interesting new effects can occur. For example, if 
a topological insulator is realized by transition metal com- 
pounds with strong electron correlation effects, antiferromag- 
netic (AFM) long-range order may develop in this material. 
Since the AFM order breaks TR symmetry and inversion sym- 
metry, the magnetoelectric coefficient P 3 defined in Eq. d40b 
deviates from its quantized value n+1/2. Denoting by n(r, t) 
the AFM Neel vector, we have P 3 (n) = P 3 (n = 0) + SP 3 (n) 
where P 3 (n = 0) is the quantized value of P 3 in the ab- 
sence of AFM order. This change in the magnetoelectric co- 
efficient has interesting consequences when spin-wave excita- 
tions are considered. Fluctuations of the Neel vector <5n(r, t) 
induce in general fluctuations of r5i-3, leading to a coupling 
betwe en spin- waves and the electromagnetic field ( Li et al.[ 
2010). In high-energy physics, su ch a particle coupled to 
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the E ■ B term is called an "axion" dPeccei and Ouinnl 1 1977 ; 
WilczekL 120091) . Physically, in a background magnetic field 
such an "axionic" spin-wave is coupled to the electric field 
with a coupling constant tunable by the magnetic field. Con- 
sequently, a polariton can be formed by the hybridization of 
the spin-wave a nd photon, similar to the p olariton formed by 
optical phonons ( Mills and Bursteinl 1974 ). The polariton gap 
is controlled by the magnetic field, which may realize a tun- 
able optical modulator. 

Another interesting effect emerges from electron correla- 
tions when a thin film of topological insulator is considered. 
When the film is thick enough so that there is no direct tun- 
neling between the surface states on the top and bottom sur- 
faces, but not too thick so that the long-range Coulomb inter- 
action between the two surfaces are still important, a inter- 
surface p^ticjej^iojeexcitation, i.e. an exciton, can be in- 



duced dSeradieh et al. 120091) . Denoting the fermion annihi- 
lation operator on the two surfaces by ipi, ip2, the exciton 
creation operator is ip\tp2- In particular, when the two sur- 
faces have opposite Fermi energy with respect to the Dirac 
point, there is nesting between the two Fermi surfaces, which 
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leads to an instability towards exciton condensation. In the 
exciton condensate phase, the exciton creation operator ac- 
quires a nonzero expectation value (ip\ip2) 7^ 0, which cor- 
responds to an effective inter-surface tunneling. Interest- 
ingly, one can consider a vortex in this exciton condensate. 
Such a vortex corresponds to a complex spatially-dependent 
inter-surface tunneling amplitude, and is equivalent to a mag- 
netic monopole. According to the Witten effect mentioned 
above ( Rosenberg and Franz . 2010h . such a vortex of the ex- 
citon condensate carries charge ±e/2, which provides a way 
to test the Witten effect in the absence of a real magnetic 
monopole. 

Besides the effects discussed above, there are many other 
physical effects related to the TME effect, or surface half QH 
effect. When a magnetic layer is deposited on top of the topo- 
logical insulator surface, the surface states can be gapped and 
a half QH effect is induced. In other words, the magnetic 
moment of the magnetic layer determines the Hall response 
of the surface states, which can be considered as a coupling 
betw een the magnet ic moment and the surface electric cur- 
rent dOi et all |2008). Such a coupling leads to the inverse 
of the half QH effect, which means that a charge current on 
the s urface can flip the magn etic moment of the magnetic 
layer ( Garate and Franz . 2010h . Similar to such a coupling be- 
tween charge current and magnetic moment, a charge density 
is cou pled to magnetic textures suc h as domain walls and vor 



u p 

tices dNomura and Nagaosa , 2010h . This effect can be used 



to drive magnetic textures by electric fields. These effects on 
a topological insulator surface coupled with magnetic layers 
are relevant to potential applications of topological insulators 
in designing new spintronics devices. 



E. Experimental results 

1 . Material growth 

There have been many interesting theoretical proposals for 
novel effects in topological insulators, but perhaps the most 
exciting aspect of the field is the rapid increase in experimen- 
tal efforts focussed on topological insulators. High-quality 
materials are being produced in several groups around the 
world, and of all different types. Bulk materials were first 
grown for experiments on topological insulators in the Cava 
grou p at Princeton Un iversity including the Bii-^Sb.,, al- 
loy ( Hsieh et ali [2008b and Bi 2 Se 3 , Bi 2 Te3, Sb 2 Te 3 crys- 



tals (IHsieh et all 12009b IXia et all I2009I) . Crystalline sam- 
ples of Bi2Te3 have also been grown at St anford Univer- 

In ad dition 
2010t 



2009h. 



sity in the Fisher group (IChen et al 
to bulk samples, Bi2Sea nanoribbo ns dHong et al. 



Kong et ali \20ld IPeng et ali |2010|) have been fabricated 
in the Cui group at Stanford University, and thin films of 
Bi2Se3 and Bi 2 Te3 have been grown by MBE by the Xue 



2009 



group at Isinghua University (ILi et al. 
2009I) , as well as other groups dli et al. 
2009I) . Thin films can also be obtained by exfoliation 



2010; 



Zhang et al 



Zhang et al. 



from bulk samples dHong et all l2010t IShahil et ali l2010t 



Teweldebrhan et ali 120101) . The stoichiometric compounds 



Bi2Se3, Bi2Te3, Sb2Te3 are not extremely difficult to grow, 
which should allow more experimental groups to have access 



to high-quality topologic al insulator samples dButch et al 
20101 Zhang et all 20091) . due to intrinsic dopi ng from va 



cancy and anti-site defects, Bi 2 Se,3 and B i2Te3 dChen et al 



2009; Hsieh et al., 2009; Xia et al. 



20091) are shown to con 



( Hsieh etali 120091) is 



tain n-type carriers while Sb2Te3 
type. Consequently, controllable extrinsic doping is re- 
quired to tune the Fermi energy to the Dirac point of 
the surfa ce states. For e xample, Bi2Te3 can be doped 
wit h Sn (Chen et ali 20091) an d Bi2Se3 can be doped with 
Sb dAnalvtis et aZ.U2010l) ~or Ca dHor et aZ.lEoiotlWang eTal. . 
2010h . Furthermore, it is found that dopin g Bi2S e,3 with 



Cu can induce superconductivity ( Hor et al. , 2010|). while 



Fe and Mn dopants may yield ferromagnetism dCha et al. 
20 lot IChen et ali \20ld Oor et ali l201oT 
Xia et q/.U2008l 



Wrav et al., 201C; 



2. Angle-resolved photoemission spectroscopy 

ARPES experiments are uniquely positioned to detect the 
topological surface states. The first experiments on topo- 
logical insulators we re ARPES experim ents carried out on 



the Bii_ x Sba; alloy dHsieh et ail 120081) . The observation 
of five branches of surface states, together with the re- 
spective spin polarizations determined later by spin-resolved 
ARPES dHsieh et q/.l 120091) . confirms the nontrivial topologi- 
cal nature of the surface states of Bi i-^Sbj;. 

2009h 



ARP ES work on Bi 9 Sec, AXiaetal. 



and 

Bi 2 Te 3 dChen et ali 120091; IHsieh et ali l2009h ~ soon followed. 
Unlike the multiple branches of surface states observed for 
Bii-^Sbj;, these experiments report a remarkably simple 
surface state spectrum with a single Dirac cone located at the 
r point and a large bulk band gap, in accordance with the 
theoretical predictions. For Bi 2 Se3, a single Dirac cone with 
linear dispersion is clearly shown at the T point within the 
band gap in Fig. 1251a) and (b). Figurel25Td) shows the y com- 
ponent of the spin polarization along the k x (T — M) dir ection 



measured by spin-resolved ARPES dHsieh et al. , 20091) . The 
opposite spin polarization in the y direction for opposite k 
indicates the helical nature of the spin polarization for surface 
states. As discussed above, Bi2Se3 has a finite density of 
n-type carriers due to intrinsic doping. Therefore, the above 
ARPES data [Fig. l25Ta).(b)l shows that the Fermi energy is 
above the conduction band bottom and the sample is, in fact, 
a metal rather than an insulator in the bulk. To obtain a true 
topological insulating state with the Fermi energy tuned into 
the bulk gap, careful control of external doping is required. 
Such control was first reported by Chen e t al. |Eig .|26] for a 
sample of Bi 2 Te 3 with 0.67% Sn doping dChen et ali l2009h . 



Some recent work on Sb2Te 3 (Hsi eh et al. 



20091) supports the 



theoretica l prediction that thi s material is also a topological 
insulator dZhang et al. , 120091) . This family of materials is 
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FIG. 25 ARPES data for the dispersion of the surface states of 
Bi2Se3, along directions (a) T — M and (b) T — K in the surface 
Brillioun zone. Spin-resolved ARPES data is shown along T — M for 
a fixed energy in (d), from which t he spin polarizati on i n momentum 
space (c) can be extracted. From Ixia et al and iHsieh et ail 

I2009I . 
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FIG. 26 ARPES measurement of (a) shape of the Fermi surface and 
(b) band dispersion along the K — Y — K direction, f or Bi2Te3 nom- 
inally doped with 0.67% Sn. From lChen eTafl I2009I . 
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FIG. 27 ARPES data for Bi 2 Se 3 thin films of thickness (a) 1QL (b) 
2QL (c) 3QL (d) 5QL (e) 6QL, measured at room t emperature (QL 
stands for quintuple layer). From Zhang et a/.Ll2009l . 



moving to the forefront of research on topological insulators 
due to the large bulk gap and the simplicity of the surface 
state spectrum. 

Although the simple model dlOU captures most of the 
surface state physics of these systems, experiments report 
a hexagonal surface state Fermi surface [Fig. [26), while 
Eq. ( IIQU only describes a circular Fermi surface sufficiently 
close to the Dirac point. However, such a hexagonal warping 
effect can be easily taken into account by including an addi- 
tional term in the surface Hamiltonian which is cubic in k dFul 
l2009h . The surface Hamiltonian for Bi2Te.3 can be written 



H(k) = Eo(k) + v k (k x <jy-k y <7 x ) + -(kl + k^)<j z , (46) 



where £o(k) = k 2 /(2m*) breaks the particle-hole symmetry, 
the Dirac velocity Vu = v(l + ak 2 ) acquires a quadratic de- 
pendence on k, and A parameterizes the amount of hexagonal 
warping (Fu, 2009I) . 

In addition to its usefulness for studying bulk crystalline 
samples, ARPES has also b een used to characterize the thin 
films of Bi2Se3 and P^T ea dLi et al. J200S Sakamoto et ali 



20 lfj IZhang et all I2009I ). The thin films we re grown to initi- 
ate a study of the crossover ( Liu et al. . 20ld) from a 3D topo- 
logical insulator to a 2D QSH state (Sec. Illl.Cb . In Fig. l27l 
ARPES spectra are shown for several thicknesses of a Bi 2Se3 
thin film, which show the evolution of the surface states. 



3. Scanning tunneling microscopy 

In addition to the ARPES characterization of 3D topolog- 
ical insulators, scanning tunneling microscopy (STM) and 
scanning tunneling spectroscopy (STS) provide another kind 
of surface-sensitive technique to probe the topological sur- 
face states. A set of materials hav e been investigated i n 
STM/ST S experiments: Bi! 



(Roushan et al. 



Bia Te.s ( Alpichshev et ali l2010l: IZhang et al. 



2009) 



2009), and 



Sb ( Gomes et all I2OO9I) . (Although Sb is topologically non- 



trivial, it is a semi-metal instead of an insulator.) The compar- 
ison bet ween STM/STS and AR PES was first performed for 
Bi 2 Te 3 dAlpichshev et al. l20ld) . where it was found that the 
integrated density of states obtained from ARPES [Fig.[28ja)] 
agrees well with the differential conductance dl/dV obtained 
from STS measurements [Fig. [28jb)] . From such a compar- 
ison, different characteristic energies {Ep, Ea, Eb, Ec and 
Ed in Fig. [28} can be easily and unambiguously identified. 

Besides the linear Dirac dispersion which has already been 
well established by ARPES experiments, STM/STS can pro- 
vide further information about the topological nature of the 
surface states, such as the interferenc e patterns of impuri- 
ties or edges ( Alpichshev et ali 20ldi Gomes et al. 2009: 



(A 
. 2 



Roushan et ali I2009t IZhang et all I2009I) . When there are 



impurities on the surface of a topological insulator, the sur- 
face states will be scattered and form an interference pat- 
tern around the impurities. Fourier transforming the inter- 
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FIG. 28 Good agreement is found between (a) the integrated den- 
sity of states from ARPES and (b) a typical scanning tunneling spec- 
troscopy spectrum. Ef is the Fermi level, Ea the bottom of the 
bulk conduction band, Eb the point where the surface states become 
warped, E g the top of the bulk va lence band, and Ed the Dirac 
point. From lAlpichshev et a/.Ll20l"ol . 



ference pattern into momentum space, one can quantitatively 
extract the scattering intensity for a fixed energy and scatter- 
ing wave vector. With such information one can determine 
what types of scattering events are suppressed. Figure l29Ta) 
and (c) sho ws the interference pat tern in momen tum space for 
Bi r Sb i_ T (IRoushan et HI l2009h and Bi 2 Te3 (IZhang et all 
2009b . respectivel y . In order to analyze the interference 
pattern dLee et all |2009), we take Bi2Te3 as an example 
[Fig. l29l c).(d)l. The surface Fermi surface of Bi2Te3 is shown 



in Fig. l29f d), for which the possible scattering events are dom- 
inated by the wave vectors qi along the K direction, q 2 along 
the M direction and q3 between the K and M directions. 
However, from Fig. |29l c) we see that there is a peak along 
the r — M direction, while scattering along the f — K di- 
rection is suppressed. This observation coincides with the 
theoretical prediction that backscattering between k and — k 
is forbidden due to TR symmetry, which supports the topo- 
logical nature of the surface sta tes. Other related theoreti- 
cal analysis are also consistent dBiswas and Balatskv , 2010t 
Guo and Franz, 1201 fj Zho u et al 1 12009b . A similar analysis 
can be applied to the surface of B^Sbi-^, and the obtained 
pattern [Fig. l297b)l also agrees well with the experimental 



data [Fig.^a)] dRoushan et aZ.[|2009l) . More recently, STM 



experiments have further demonstrated that the topological 
surface states ca n penetrate barr iers while maintaining their 



extended nature ( Sep et al. , 2010h 
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FIG. 29 (a) Measured interference pattern in momentum space for 
impurities on the surface of Bi^Sbi-a;. (b) Pattern calculated from 
ARPES data on Bi^Sbi-a;, which agrees well with the interference 
pattern in (a), (c) Similar interference pa ttern and (d) po s sible scat- 
tering wave vecto rs for Bi2Te3. From IRoushan et a/1 120091 and 
IZhang et aZ.ll2009l . 



is th e observation of s urface state LLs in a magnetic 
field ( Cheng et al. , 2010t Hanaguri et al. , 2010l) . As shown 
in Fig.[30ia) and (b), discrete LLs appear as a series of peaks 
in the differential conductance spectrum (dl/dV), which sup- 
ports the 2D nature of the surface states. Further analysis on 
the dependence of the LLs on the magnetic field B shows that 
the energy of the LLs is proportional to \JnB where n is the 
Landau level index, instead of the usual linear-in-B depen- 
dence. This unusual dependence provides additional evidence 
for the existence of surface states consisting of massless Dirac 
fermions. Furthermore, the narrow peaks in the spectrum also 
indicate the good quality of the sample surface. 



4. Transport 



In addition to the above surface-sensitive techniques, 
large effort has been devoted to transport measure- 



a 

ments 



inclu ding dc transpor t (Analytis et al. 



Butch et al. 



Steinberg et al. 



2010^ 

2010T 



Chen et al. . £010; 



Tims et all 120 1C ) 



ments in the microwave dA nalvtis et al. 



frared regime s (|Butch et al. 



Eto et al. 



Sushkov et all 120 lOh . which are necessary steps towards 
the direct measurement of topological effects such as the 
TME, and for future device applications. However, transport 
experiments on topological insulators turn out to be much 
more difficult than surface-sensitive measurements such 
as ARPES and STM. The main difficulty arises from the 
existence of a finite residual bulk carrier density. Materials 
such as Bii-^Sbz or Bi2Se3 are predicted to be topological 



20101: 



2010; 



2010; 



and 
20101) 



LaForge et al. 



measure- 
an d in- 
20101 



Another important result of STM/STS measurements 
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FIG. 30 Tunneling spectra for the surface of Bi2Se;s in a magnetic 
field, showing a serie s of peaks attributed to t he occurrence of surfa ce 
Landau levels. From lCheng et a/ll2010l and lHanaguri et fl/.Ll20ld . 
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FIG. 3 1 Angular dependence of the cyclotron resonance field for a 
71 GHz microwave. The inset shows examples of transmission data 
offset for clarity (top to botto m: 90° to 0° in steps of 10°). From 
lAvala-Valenzuela et flZll2010l . 



insulators if they are perfectly crystalline. However, real 
materials always have impurities and defects such as anti-sites 
and vacancies. Therefore, as-grown materials are not truly 
insulating but have a finite bulk carrier density. As discussed 
in Sec. IIII.E.21 such a residual bulk carrier de nsity is also 



observed in ARPES for Bi2Se3 and Bi2Te3 dChen et al. 



20091: iHsieh et all 12009b . From the ARPES results it seems 
that the residual carrier dens i ty can be compensated for by 
chemical doping ( Chen et all l2009t) . Nevertheless, in trans- 
port experiments the compensation of bulk carriers appears to 
be much more difficult. Even samples which appear as bulk 
insulators in ARPES experiments still exhibit s ome finite bulk 



carrie r density in transport measurements (lAnalvtis et al. 



20101) . which suggests the existence of an offset between 
bulk and surface Fermi levels. Another difficulty in transport 
measurements is that a cleaved surface rapidly becomes 
heavily n-doped when exposed to air. This leads to further 
discrepancies between the surface condition observed in 
transport and surface-sensitive measurements. 

In spite of the complexity described above, the signature 
of 2D surface states in transport experiments has been re 



cently reported (lAnalvtis et lAvala- Valenzue la et al. 

2010h . For example, Fig. [3~T| shows results obtained 



by microwave spectroscopy ( Ayala- Valenzuela et al. . 2010h 



on Bi2Se3, where it is found that the cyclotron res- 
onance frequency only scales with perpendicular mag- 
netic field B±, suggesting the 2D nature of the reso- 
nance. Similarly, the dependence of Shubnikov-de Haas 
oscillations on the angle of the magnetic field can help 
to distinguish the 2D surface states from the 3D bulk 
states, bot h for Bi2Se3 ( Ayala- Valenzuela et al. . 2010b and 
Bi^Sbi-^ ( Taskin and Andol " 20091) . Signatures of the topo- 
logical surface states have also been searc hed for in the 



tempe r ature dependence of the resist ance (lAnalvtis et al. 



20101; ICheckelskv et all l20l61 l2009h . the magnetoresis- 
tance ( Tang et al. . 2010l) . and weak antilocalization ef- 



fects dCheckelskv et all \2Q\& IChen et all |2010|) . However, 
besides the experiments mentioned above with positive ev- 
idence for the existence of topological surface states, some 
experiments show that t he transport data can be entirely ex - 
plained by bulk carriers ( Butch et al. . 2Q\(k Eto et al. . 2010h . 
The resolution of this controversy requires further improve- 
ments in experiments and sample quality. 

To reach the intrinsic topological insulator state without 
bulk carriers, various efforts have been made to reduce the 
bulk carrier density. One approach consists in compensat- 
ing the b ulk carriers by chemi cal d oping, e.g. dopin g Bi2Se3 
with Sb (lAnalvtis et all l2010l). Ca dHor et all l2009h . or dop- 



ing Bi2Te3 with Sn (IChen et 

oE l2009h . Although chem 



ical doping is an efficient way to reduce the bulk carrier 
density, the mobility will be usually reduced due to foreign 
dopants. However, we note that the substitution of the isova- 
lent Bi with Sb can reduce the carrier density but still keep 



high mobilities ( Analvtis et al. . 2010b . Also, it is difficult 



to achieve accurate tuning of the carrier density by chem- 
ical doping, because each different chemical doping level 
needs to be reached by growing a new sample. The sec- 
ond method consists in suppressing the contribution of bulk 
carriers to transport by reducing the sample si ze down to 



the nanoscale, such as quasi- ID nanoribbons dKong et al. 



20 lOt IPeng et all l2010t ISteinberg etall l2010l), or quasi-2D 
thin fi lm dCheckelskv et a/.ll2010l:IChen et a/.Ll2010l:lLi etal. 



2009). In Fig. [32] the magnetoresistance of a nanorib- 



bon exhibits a primary hc/e oscillation, which corresponds 
to Aharonov-Bohm oscillations of the surfa c e state around 
the surface of the nanor ibbon dPeng et all l2010h . This 
oscillation also indicates that the bulk carrier density has 
been reduced greatly so that the contribution of the surface 
states can be observed. The Aharon ov-Bohm oscillation has 
also been investigated theore tically dBardarson et 
Zhang and Vishwanath , 2010l) . An important advantage of 




FIG. 32 Magnetoresistance for fields up to ±9 T. Left inset: 
magnetic fields at which well-developed resistance minima are ob- 
served. Right inset: fast Fourier transform of the resistance deriva- 
tive dR/dB, whe re peaks correspon d to he/ e and hc/2e oscillations 
are labeled. From lPeng et a/.LboiOl 



a sample of mesoscopic size is the possibility of tuning the 
carrier density by an external gate voltage. Gate control of 
the carrier density is rather important because it can tune the 
bulk carrier density continuously while preserving the quality 
of the sample. Gate control of the carrier density has been 



indeed observed in nanoribbons of Bi2Se.3 (Steinberg et al 



2010), mechanically exfoliated thin fil ms (IChecke 
2010 ). or epitaxially grown thin films ( Chen et oIT 



sky et al 



2010). In 



particular, tuning of the carrier polarity fro m n-type to p-type 
has been reported ( Checkelskv et al. , 2010l) . where the change 
in polarity corresponds to a sign change of the Hall resistance 
R yx in a magnetic field [Fig. l33l b)l. 



F. Other topological insulator materials 

The topological materials HgTe, Bi2Se3, Bi2Te3 and 
Sb2Te3 not only provide us with a prototype material for 2D 
and 3D topological insulators, but also give us a rule of thumb 
to search for new topological insulator materials. The non- 
trivial topological property of topological insulators originates 
from the inverted band structure induced by SOC. Therefore, 
it is more likely to find topological insulators in materials 
which consist of covalent compounds with narrow band gaps 
and heavy atoms with strong SOC. Following such a guiding 
principle, a large number of topological insulator materials 
have been proposed recently, which can be roughly classified 
into several different groups. 

The first group is similar to the tetradymite semiconduc- 
tors, where the atomic p-orbitals of Bi or Sb play an essential 
role. Thallium-based III-V-VI2 ternary chalcogenides, includ- 
ing T lBiQ 2 and TlSbQ 2 with Q = Te, S e and S, belong to this 
class dLin et ^Zl l2010t Yan et 0/1 120101) . These materials have 
the same rhombohedral crystal structure (space group -D| d ) as 
the tetradymite semiconductors, but are genuinely 3D, in con- 
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FIG. 33 (a) Gate voltage dependence of the resistance in zero mag- 
netic field for different temperatures, (b) Gate voltage dependence 
of the Hall resistan ce R vx in a magnetic fi eld of 5 T for different 
temperatures. From lCheckelskv et q/.Ll2O10l . 



trast to the layered tetradymite compounds. These materials 
have recen tly been experimentally observed to be topological 



have recen tly been ex perimentally obser ved to 
insulators dChen et artbOlOtlSato et q/.ll2010l) 



A typical material of the second group is distorted bulk 
HgTe. In contrast to conventional zincblende semiconductors, 
HgTe has an inverted bulk band structure with the r§ band 
being higher in energy than the Tq band. However, HgTe by 
itself is a semi-metal with the Fermi energy at the touching 
point between the light-hole and heavy-hole r§ bands. Con- 
sequently, in order to get a topological insulator, the crys- 
tal structure of HgTe should be distorted along the [111] di- 
rection to_o£ena gapbetween the heavy-hole and light-hole 
bands dDai et ail I2008I) . A si milar band structure also exist s 
in tern ary Heusler compounds ( Chadov et al. , 2010t Lin et al. 



20 id) , and around fifty of them are found to exhibit band in 
version. These materials become 3D topological insulators 
upon distortion, or they can be grown in quantum well form 
similar to HgTe/CdTe to realize the 2D or the QSH insula- 
tors. Due to the diversity of Heusler materials, multifunc- 
tional topological insulators can be realized with additional 
properties ranging from superconductivity to magnetism and 
heavy-fermion behavior. 

Besides the above two large groups of materials, there are 
also some other theoretical proposals of new topological insu- 
lator materials with electron correlation effects. An example 
is the case of Ir-ba sed materials. The QS H effect has been pro- 
posed in Na2lr03 dShitade et a/ll2009l) . and topological Mott 
insulator phases have been proposed in Ir-based pyrochlore 
oxides Ln2lr207 with L n = Nd, Pr 



Guo andFranzl I2009 : 



2010I: lYangandKim . 



Pesin and BalentsL I2010t Wan et al. 

2010h . Furthermore, a topological structure has also been 
considered in Kond o insulators, with a possible realization in 
SmB 6 and CeNiSn dDzero et a/.Ll2Q10b . 



IV. GENERAL THEORY OF TOPOLOGICAL INSULATORS 



The TFT dOi et a/.[|2008l) and t h e TBT dFu and Kanell2007 



2007 



Fu et ali \2QQH iKane and MeleL I2005 



Moore and Balents, 



RovU2009li " o different general theories of the 
topological insulators. The TBT is valid for the non- 
interacting system without disorder. The TBT has given sim- 
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pie and important criteria to evaluate which band insulators 
are topologically non-trivial. The TFT is generally valid for 
interacting systems including disorder, and the it identifies the 
physical response associated with the topological order. Re- 
markably, the TFT reduces exactly to the TBT in the non- 
interacting limit. In this section, we review both general theo- 
ries, and also discuss their connections. 



A. Topological field theory 

We are generally interested in the long-wavelength and low- 
energy properties of a condensed matter system. In this case, 
the details of the microscopic Hamiltonian are not important, 
and we would like to capture essential physical properties in 
terms of a low-energy effective field theory. For conventional 
broken-symmetry states, the low-energy effective field theory 
is fully deter mined by the orde r parameter, symmetry and di- 
mensionality dAndersonl 1 1 997b . Topological states of quan- 
tum matter are similarly described by a low-energy effective 
field theory. In this case, the effective field theory generally 
involve topological terms which describe the universal topo- 
logical properties of the state. The coefficient of the topolog- 
ical term can be generally identified as the topological order 
parameter of the system. A su ccessful example is the TFT 
of the QH effect (IZhangL 119921) . which captures the univer- 
sal topological properties such as the quantization of the Hall 
conductance, the fractional charge and statistics of the quasi- 
particles and the ground state degeneracy on a topologically 
nontrivial spatial manifold. In this section, we shall describe 
the TFT of the TR invariant topological insulators. 



1 . Chern-Simons insulator in 2 + 1 dimensions 

We start from the previously men tioned QH sy stem in (2 
1)D, the TFT for which is given as ( lZhanall992l) 

Ci 

4.7T 



S e s 



d 2 x 



dtA„e^ T d v A T , 



(47) 



where the coefficient C\ is generally given by (IWang et al. 
2010h 



can be calcul ated explicitly from a single Feynman diagram 
[Fig [34Ta)1 ( Golterman et ali 1 1 9931: iNiemi and Semenoff , 
19831: loief all l2008t IVolovikl 120021) . and one obtains d as 
given in Eq. ( 1481 . but with G replaced by the noninteract- 
ing Green's function Go- Carrying out the frequency integral 
explicitly, one obtains the TK NN invariant expresse d as an 
integral of the Berry curvature dThouless et al.U 1982b . 



with 



/xj/(k) : 
a l (k) = 



1 

2V 



dk x I dkyf X y(k.) G Z, 



(50) 



da y (k) da x (k) 



dk x 



dk v 



a£ occ 



\ak) 



cy- 



linder TR, we have Aq — > Aq, A — > —A, from which we 
see that the (2 + 1)D CS field theory in Eq. d47]i breaks TR 
symmetry. All the low-energy response of the QH system can 
be derived from this TFT. For instance, from the effective La- 
grangian in Eq. d47b . taking a functional derivative with re- 
spective to A^, we obtain the current 



3p. 



27T 



e^ T d v A T . 



The spatial component of this current is given by 



3i 



2tt 



(51) 



(52) 



(53) 



while the temporal component is given by 

This is exactly the QH response with Hall conductance an = 
Ci/(2ir), implying that an electric field induces a transverse 
current, and a magnetic field induces charge accumulation. 
The Maxwell term contains more derivatives than the CS term 
and is therefore less relevant at low energies in the renormal- 
ization group sense. Therefore, all the topological response 
of the QH state is exactly contained in the low-energy TFT of 
Eq. d47]i. 



d = | J ^^TT^PGd^Gd^GdpG- 1 ]^) 2 - Chern-Simons insulator in 4 + 1 dimensions 



and G(k) = G(k, ui) is the imaginary-time single-particle 
Green's function of a fully interacting insulator, and {1,1/, p = 
0,1,2 = t, x, y. For a general interacting system, assuming 
that G is nonsingular, we have a map from the three dimen- 
sional momentum space to the space of nonsingular Green's 
functions, belonging to the group GL( n, C), whose third ho- 
motopy group is labeled by an integer dWang et a/.Ll2010h : 



7r 3 (GL(n,C)) £* 



(49) 



The TFT of the QH effect does not only capture the univer- 
sal low-energy physics, it also points out a way to generalize 
the TR symmetry breaking QH state to TR invariant topolog- 
ical states. The CS field theory can be ge neralized to all odd 
dimensional spacetimes ( Nakahara , 1990b . This observation 
lead Zhang a nd Hu to discover a g eneralization of the QH in- 
sulator state dZhang and Hu , 2001 ) which is TR invariant, and 
defined in (4 + 1)D. It is the fundamental TR invariant in- 
sulator state from which all the lower-dimensional cases are 



The winding number for this homotopy class is exactly mea- 
sured by G\ defined in Eq. d48l . Here n > 3 is the num- 
ber of bands. In the noninteracting limit, C\ in Eq.d47li 



derived, and is described by the TFT dBernevig et al. 

C 2 



2002) 



5, 



eft 



24tt 2 



<?xdbe! lvpar A ll d u A p d a A T . 



(54) 
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(a) 



(b) 



FIG. 34 Fermion loop diagrams leading to the Chern-Simons term, 
(a) The (2 + 1)D Chern-Simons term is calculated from a loop dia- 
gram with two external photon lines, (b) The (4+ 1)D Chern-Simons 
term is calculated from a loop diagram with three external photon 
lines. 



Under TR, we have Ao — > Aq, A — > —A, and this term is 
explicitly TR invariant. Generally, the coefficient C 2 is ex- 
pressed injermsrf theGreen's function of an interacting sys- 
tem as ( Wang et al. , 201fj) 



Co 



15 



d 5 k 



Tr[e> lv i" TT Gd ll G- 



l Gd v G- 



Gd p G- 



(2tt) 5 

xGdvG^GdrG- 1 }, (55) 
which labels the homotopy group dWang et q/.[|2010h 

7r 5 (GL(n,C)) Z, (56) 

similarly to the case of the (2 + 1)D CS term. For a nonin- 
teracting system, Ci can be calculated from a single Feynman 
diagram [Fig. [34l" b)l and one obtains C2 as given in Eq. (f55b . 
with G replaced by the noninteracting Green's function Go- 
Explicit integration over the frequency gives the second Chern 



number (lOi et ail 12008) 



Co 



1 



32tt 2 



(57) 



with 



a/3 



a/3 



(k) = -i (a,k| — |/3, k 



d_ 



+ i \a%, a. 



,aj3 



where k,£= 1, 2, 3,4 = x, y, z, w). 

Unlike the (2 + 1)D case, the CS term is less relevant than 
the non-topological Maxwell term in (4 + 1)D, but is still of 
primary importance when understanding topological phenom- 
ena such as the chiral anomaly in a (3+l)D system, which can 
be reg arded as the boundary of a (4 + 1)D system dOi et al. , 
2008b . Similar to the (2 + 1)D QH case, the physical response 
of (4 + 1)D CS insulators is given by 

C 2 



~-° jUi/pCTTQ A f) A 
Z7T 



(58) 



which is the nonlinear response to the external field A„. To 
understand th is response bet ter, we consider a special field 
configuration dOi et a/.ll2008l) : 



A :i 



0, A 



B z x, A z 



-E z t, A u 



A t = 0, (59) 



where x, y, z, w are spatial coordinates and t is time. The only 
nonvanishing components of the field strength are F xy = B z 
and F zt = —E z . According to Eq. d58l , this field configura- 
tion induces the current 

If we integrate the equation above over the x, y dimensions, 
with periodic boundary conditions and assuming that E z does 
not depend on x, y, we obtain 



dxdy j v 



C2_ 

4tt 2 



dxdyB z E. 



C 2 N X 
2tt 



: ^,(60) 



where N xy = J dxdyB z /2ir is the number of flux quanta 
through the xy plane, which is always quantized to be an in- 
teger. This is exa ctly the 4D generaliza tion of the QH effect 
mentioned earlier dZhang and HuluOOll) . Therefore, from this 
example we can understand the physical response associated 
with a nonvanishing second Chern number. In a (4 + 1)D in- 
sulator with second Chern number C2, a quantized Hall con- 
ductance C2N xy /2ir in the zw plane is induced by a magnetic 
field with flux 2irN xy in the perpendicular {xy) plane. 

We have discussed the CS insulators in (2 + 1)D and 
(4 + 1)D. Actually, these discussions can be straightforwardly 
generalized to higher dimensions. In doing so, it is worth 
noting that there is a even-odd alternation of the homotopy 
groups of GL(n, C): we have Tt2k+i (GL(n, C)) = Z, while 
7T2fc(GL(n, C)) = 0. This is the mathematical mechanism 
underlying the fact that CS insulators of integer class exist in 
even spatial dimensions, but do not exist in odd spatial di- 
mensions. Reduced to noninteracting insulators, this becomes 
the alternation of Chern numbers. As a number characteristic 
of complex fiber bundles, Chern numbers exist only in even 
spatial dimensions. This is an example of the relationship 
between homotopy theory and homology theory. We shall 
see another example of this relationship in the following sec- 
tion: both the Wess-Zumino-Witten (WZW) terms and the CS 
terms are well-defined only modulo an integer. 



3. Dimensional reduction to the three-dimensional Z 2 
topological insulator 

The 4D generalization of the QH effect gives the fundamen 
tal TR invariant topological insulator from which all lower- 
dimensional topological insulators can be derived s ystemati 



cally by a procedure called dimensional reduction (lOi et a 
20081) . Starting from the (4 + 1)D CS field theory in 



Eq. d54l i. we consider field configurations where A ll (x) 
A ll (xo,xi,X2,X3) is independent of the "extra dimension" 
X4 = w, for n = (xo, xi, X2, and A4 = A w depending 
on all coordinates (xq, X\, X2, £3, We consider the ge- 
ometry where the "extra dimension" X4 forms a small circle. 
In this case, the X4 integral in Eq. d54l ) can be carried out ex- 
plicitly. After restoring the unit of electron charge e and flux 
hc/e following the convention in electrodynamics, we obtain 
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an effective TFT in (3 + 1)D: 



32tt 2 



d 3 xdt 9(x, t)e» vpT F^F^x, t), (61) 



where a = e 2 /He ~ 1/137 is the fine structure constant and 

9(x,t) = C20(x,i) = C2 f dxi A±(x, t, £4), (62) 



which can be interpreted as the flux due to the gauge field 
Ai(x, t, X4) through the compact extra dimension [Fig. [331 . 
The fiel d 9(x, t) is calle d the axion field in the field theory lit- 
erature dWilczek . 1987 ). In order to preserve the spatial and 
temporal translation symmetry, 9(x, t) can be chosen as a con- 
stant parameter rather than a field. Furthermore, we already 
explained that the original (4 + 1)D CS TFT is TR invariant. 
Therefore, it is natural to ask how can TR symmetry be pre- 
served in the dimensional reduction. If we choose C2 = 1, 
then 9 = 4> is the magnetic flux threading the compactified 
circle, and the physics should be invariant under a shift of 9 
by 2tt. TR transforms 9 to —0. Therefore, there are two and 
only two values of 9 which are consistent with TR symmetry, 
namely 9 = and 9 = n. In the latter case, TR transforms 
9 = it to 9 = — it, which is equivalent to 9 = n mod 2ir. We 
therefore conclude that there are two different classes of TR 
invariant topological insulators in 3D, the topologically trivial 
class with = 0, and the topologically nontrivial class with 

9 = 7T. 

As just seen, it is most natural to view the 3D topological 
insulator as a dimensionally reduced version of the 4D topo- 
logical insulator. However, for most physical systems in 3D, 
we are generally given an interacting Hamiltonian, and would 
like to define a topological order parameter that can be eval- 
uated directly for any 3D model Hamiltonian. Since the 9 
angle can only take the two values and it in the presence of 
TR symmetry, it can be naturally defined as the topological 
order par ameter itself . For a generally interacting system, it is 
given by ( Wang et al . 2010l) : 



xGdpG^GdaG^GduG- 1 }, 



(63) 



where the momentum k = (ki, fe, k^) is integrated over the 
3D Brillouin zone and the frequency ko is integrated over 
(-00, +00). G(k,u = 0) = G(fc ,k,w = 0) = G(/c ,k) 
is the imaginary-time single-particle Green's function of the 
fully interacting many-body system, and G(k,u) for 11 ^ 
is a smooth extension of G(k,u = 0), with a fixed reference 
value G(k, u = 1) corresponding to the Green's function of a 
topologically trivial insulating state. G(k, u = 1) can be cho- 
sen as a diagonal matrix with G aa = (iko — A) -1 for empty 
bands a and Gpp = (iko + A) -1 for filled bands (3, where 
A > is independent of k. Even th ough P3 is a physical 
quantity in 3D, a WZW ( Witteniri983 ) type of extension pa- 



rameter u is introduced in its definition, which plays the role 
of kn in the formula (l55i l defining the 4D TI. The definition in 





FIG. 35 Dimensional reduction from (4 + 1)D to (3 + 1)D. The 
X4 direction is compactified into a small circle, with a finite flux (j> 
threading through the circle due to the gauge field A4,. To preserve 
TR symmetry, the total flux can be either or 7r, resulting in a Z2 
classification of 3D topological insulators. 



momentum-space is analogous to the real-space WZW term. 
Similar to the WZW term, P3 is only well-defined modulo an 
integer, and can only take the quantized values of or 1/2 



modu lo an integer for an TR invariant insulator ( IWang et al. 
2010h . 



Essential for t he definition of P3 in Eq. d63l is the TR in- 
variance identity ( IWang et a/.ll2010l) : 



G(fc ,-k) =TG(k ,k) T T- 



(64) 



which is crucial for the quantization of P3. Therefore, we 
see that unlike the integer-class CS insulators, the Z2 insula- 
tors are symmetry-protected topological insulators. This can 
be clearly seen from the viewpoint of the topological order 
parameter. In fact, the quantization of P3 is protected by TR 
symmetry. In other words, if TR symmetry in Eq. d64l is bro- 
ken, then P3 can be tuned continuously, and can be adiabati- 
cally connected from 1/2 to 0. This is fundamentally differ- 
ent from the CS insulators, for which the coefficient of the 
CS term given by Eq. ( l55l l is always quantized to be an inte- 
ger, regardless of the presence or absence of symmetries. This 
integer, if nonzero, cannot be smoothly connected to zero pro- 
vided that the energy gap remains open. There exists a more 
exhaustive classifica tion scheme for topological insulators in 
various dimensions ( Kitaev , 20091 : Oi et al. , 2008 ; Ryu et al. 



20101) which takes into account the constraints imposed by 
various symmetries. 

For a noninteracting system, the full Green's function G in 
the expression for P3 [Eq. d63l l is replaced by the noninteract- 
ing Green's function Go- Furthermore, the frequency integral 
can be carried out explicitly. After some manipulations, one 
finds a simple and beautiful formula 

P3= vhJ ' /iA '" iT, -rA/- k ? - ^(kK^KCKftD 

which expresses P3 as the integral of the CS form over the 3D 
momentum space. For explicit models o f topological i nsula- 



tors, such as the model by Zhang et al. (IZhang et "^Zl l2009h 



discussed in Sec. IIII.A1 one can evaluate this formula explic- 
itly to obtain 



P 3 = 9/2ir = 1/2 



(66) 
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in the topologically nontrivial state (lOi et all |2008). Es- 
sin, Moore and Vand erbilt also calculat ed P3 for a variety 



of interesting models dEssin et al. , 120091) . In a generic sys 



tem without time-reversal or inversion symmetry, there may 
be non-topological contributions to the effective action (l6"TT l 
which modifies the formula of P3 given in Eq. d65b . How- 
ever, such corrections vanishes when time-reversal symme- 
try or inversion symmetry is present, so that the quantized 
valu e J-3 = 1/2 (mod 1) in topological insulator remains ro- 
bust dEssin et q/.lkoiOtlMalashevich et q/.Ll2010l) . 

Similar to the case of (4+l)D CS insulators, there is also an 
important difference between the 9 term for (3+ 1)D topologi- 
cal insulators and the (2 +l)DCS term for QH s ystems, which 
we shall briefly discuss dMaciejko et all I2OI0I) . In (2 + 1)D, 
the topological CS term dominates over the non-topological 
Maxwell term at low energies in the renormalization group 
flow, as a simple result of dimensional analysis. However, in 
(3 + 1)D, the 9 term has the same scaling dimension as the 
Maxwell term, and is therefore equally important at low ener- 
gies. The full set of modified Maxwell's equations including 
the topological term d62l is given by 

^-d v F^ + d v V» v + ^e^d„(P 3 F aT ) = -j*,(67) 

47T 47T C 

which can be written in component form as 

V-D = 47rp + 2a(VP 3 -B), 

VxH-^ = ^j-2 Q ((VP 3 x E) + - c (d t P 3 ) b) , 
1<9B 

VxE+-— = 0, 

c ot 

V • B = 0, 



(68) 



where D = E + AttP and H = B - AnM only include the 
non-topological contributions. Alternatively, one can use the 
ordinary Maxwell's equations with modified constituent equa- 
tions d40b . These set of modified Maxwell equations are ca lled 
the axion electrodynamics in field theory dWilczeklll987l) . 

Even though the conventional Maxwell term and the topo- 
logical term are both present, there exist experimental designs 
which c an in principle extract the purely topological contri- 
butions ( Macieiko et ah , 20101) . Furthermore, the topological 
response is completely captured by the TFT, which we shall 
discuss. Starting from the TFT [Eq. [6TJ, we take a functional 
derivative with respect to A^, and obtain the current as 



f = ^e^d u P 3 d a A T , 



(69) 



which is the general topological response of (3 + 1)D insula- 
tors. It is worth noting that we do not assume TR invariance 
here, otherwise P 3 should be quantized to be integer or half- 
integer. In fact, here we assume a inhomogeneous P 3 (x,t). 
It is interesting to notice that this electromagnetic response 
looks very similar to the 4D response in Eq. [58] with the only 
difference that A4 is replaced by P3 in Eq. (1691 . This is a man- 
ifestation of dimensional reduction at the level of the electro- 
magnetic response. A more systematic treatment on this topic 



in phase space can be also be performed dOi et a/.l 120081) . The 
physical consequences Eq. (|69l can be understood by study- 
ing the following two cases. 

1. Half-QH effect on the surface of a 3D topological insula- 
tor. Consider a system in which P 3 = P 3 (z) only depends 
on z. For example, th is can be realize d by the lattice Dirac 
model dOi et all 120081) with 9 = 9{z) dFradkin et all \l9S& 
WilczekUl987i In this case, Eq. (l69l becomes 



r 



d z P 3 
2tt 



^ vp d u A p , n,v,p = t,x,y, 



which describes a QH effect in the xy plane with Hall con- 
ductivity a xy = d z P 3 /2ir. A uniform electric field E x along 
the x direction induces a current density along the y direc- 
tion j y = (d z P 3 /2Tr)E x , the integration of which along the z 
direction gives the Hall current 



t2D _ 

Jy ~ 



d,Zjy = 



1 

2^ 



dP 3 E x 



which corresponds to a 2D QH conductance 



a™ = / dP 3 /27T. 



(70) 



For a interface between a topologically nontrivial insulator 
with P 3 = 1/2 and a topologically trivial insulator with 
P 3 = 0, which can be taken as the vacuum, the Hall conduc- 
tance is an = AP3 = ±1/2. Aside from an integer ambi- 
guity, the QH conductance is exactly quantized, independent 
of the details of the interface. As discussed in section llll.D.ll 
the half quantum Hall effect on the surface is a reflection of 
the bulk topology with P3 = 1/2, and can not be determined 
purely from the low energy surface models. 
2. Topological magnetoelectric effect induced by a temporal 
gradient ofP 3 . Having considered a time-independent P3, we 
now consider the case when P 3 = P 3 (t) is spatially uniform, 
but time-dependent. Equation ( |69l now becomes 



3 = 



d t P 3 
2tt ' 



tjk djA k , i,j,k = x,y,z, 



which can be simply written as 

d t P 3 

3 = -^r 



B. 



(71) 



Because the charge polarization P satisfies j = <9 t P, we can 
integrate Eq. |7T| in a static, uniform magnetic field B to get 

dtP = -d t (P 3 B/2tt), so that 



B 

'2tt 



[P 3 + const.) 



(72) 



This equation describes the charge polarization induced by a 
magnetic field, which is a magnetoelectric effect. The promi- 
nent feature here is that it is exactly quantized to a half-integer 
for a TR invariant topological ins ulator, which is called the 



topological magnetoelectric effect dOi et all 120081) (TME) 



Anothe r related effect originating from the TFT is the Wit- 
ten effect dOi et al.i 120081: IWittenl [l979h . For this discussion, 
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we assume that there are magnetic monopoles. For a uniform 
P3, Eq. (EB leads to 



Vj 



d t P 3 
2ir 



V B. 



(73) 



Even if magnetic monopole does not exist as elementary par- 
ticles, for a lattice system, the monopole density p m = V • 
B/2tt can still be nonvanishing, and we obtain 



dtPe = (d t P 3 )pr, 



(74) 



Therefore, when P3 is adiabatically changed from zero to 
/ 2tt, the magnetic monopole will acquire an electric charge 
of 



e 

— < 

2tt 



where Q m is the magnetic charge. Such a relation was 
first derived by Witten in the co ntext of the t opological 
term in quantum chromodynamics (I Witten , 1 19791) . and later 



discussed in the context of topol ogical insulators dOi et al. 



2008 ; Rosenberg and Franz , 2010l) . This effect could also ap- 
pear under a different guis e in topological exciton condensa- 
tion ( Seradjeh et al. , 20091) . where a e/2 charge is induced by 
a vortex in the exciton condensate, which serves as the "mag- 
netic monopole". 



4. Further dimensional reduction to the two-dimensional Z 2 
topological insulator 

Now we turn our attention to (2 + 1)D TR invariant Z2 
topological insulators. Similar to the WZW-type topological 
order parameter P3 in (3 + 1)D, there is also a topological 
order parameter defined for (2 + 1)D TR invariant insulators. 
The main difference between (2 + 1)D and (3 + 1)D is that in 
(2 + 1)D we need two WZW extension parameters u and v, in 
contrast to a single parameter u in (3 + 1)D, which is a man- 
ifestation of the fact that both d escend from the fund amental 
(4 + 1)D topological insulator dZhang and Hu , 2001). For a 
general interacting insulat or, the (2 + 1)D t opological order 



parameter is expressed as dWang et aZl l2010l) 



P> 



1 



[if par 



du 



dv 



d 3 k 



120 7-i J-i J (2tt) 

■xGd^G^GdpG^GdaG^GdrG- 1 
or 1/2 (modZ), 



Tr [Gc^G" 1 



(76) 



where e^ yf " JT is the totally antisymmetric tensor in five di- 
mensions, taking value 1 when the variables are an even 
permutation of (fco, kx, fa, u, v). The cases P2 — and 
P-2 = 1/2 modulo an integer correspond to topologically triv- 
ial and nontrivial TR invariant insulators in (2 + 1)D, respec- 
tively. This topological order parameter is valid for interact- 
ing QS H systems in (2 + 1)D, including states in the Mott 
regime ( Raghu et al. . 20081) . P2 can be physically measured 



by the fractional charge at the edge of the QSH state (lOi et al 
20081) . 



5. General phase diagram of topological Mott insulator and 
topological Anderson insulator 

So far, we have introduced topological order param- 
eters for TR invariant topological insulators in 4D, 3D 
and 2D. These topological order parameters are defined in 
terms of the full single-particle Green's function. A cau- 
tion in order is that these topological order parameters are 
not appli cable to fractional states with ground state de- 
generacy (Bernevig and Zhangl l2006t Levin and Stern . l2009t 
Maciejko et a/.ll2010l:ISwingle et al.[\20l(jt) . for which a TFT 
approach is still possible, but simple topological order pa- 
rameters are harder to find. In 3D, fractional topological in- 
sulators are characterized by a topol ogical order parameter 



(75) P 3 that is a rational multiple of 1/2 dMaciejko et all l2010t 



Swingle et al. . 2010h . Such states are consistent with TR 



symmetry if fractionally charged excitations and ground state 
degeneracies on spatial manifolds of nontrivial topology are 
present. When TR symmetry is broken on the surface, a frac- 
tional T ME gives rise to half of a fractional QH eff ect on the 
surface dMaciejko et oll bOlOtlSwingle et a/.lEoiol) . 

Next we shall discuss the physical consequences implied 
by the topological order parameter such as P2 and P3. The 
discussion we shall present is very general and its applica- 
bility does not depend on the spatial dimensions. Further- 
more, since the topological order parameters are expressed 
in terms of the full Green's function of an interacting sys- 
tem, they can be useful to general interacting systems. Sup- 
pose we have a family of Hamiltonians labeled by several 
para meters. To be speci fic, we consider a typical phase dia- 
gram dWang et al. 1 120101) [Fig. d36ll] for an interacting Hamil- 
tonian H = -Ho (A) + Hi(g), where Hq is the noninter- 
acting part including terms such as tijc\cj, and Hi is the 
electron-electron interaction part including terms such as the 
Hubbard interaction gn^n^. These two parts are determined 
by single-particle parameters A = (Ai, A2, ■ • • ) and coupling 
constants g = (gx,92, ■ ■ ■ )■ When (A, g) are smoothly tuned, 
the ground state also evolves smoothly, as long as the energy 
gap remains open and the topological order parameters such as 
P2 and P3 remain unchanged. Only when the gap closes and 
the full Green's function G becomes singular, these topolog- 
ical order parameters have a jump, as indicated by the curve 
ab in Fig. d36b . The most important point in Fig.[36]can be il- 
lustrated by considering the vertical line BF. Starting from a 
noninteracting state F, one adiabatically tunes on interactions, 
and eventually there is a phase transition at E to the topolog- 
ical insulator (TI) state. The interacting state B, which is an 
interaction-induced topological insulator state, has a different 
topological order parameter from the corresponding noninter- 
acting normal insulator (NI) state F. The difference in the 
topological order parameter thus provides a criterion for dis- 
tinguishing topological classes of insulators in the presence of 
general interactions. 

There have already been several theoretical proposals of 
strongly interact ing topologi c al ins ulators, i.e. topological 
Mott insulators tiRaghu et all 120081) . 2D topologically non- 
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FIG. 36 Phase diagram in the (A, g) plane. The dark curve ab is 
the phase boundary separating normal insulators (NI) and topolog- 
ical insulators (TI). All phases are gapped, except on ab. The true 
parameter space is in fact infinite dimensional , but th is 2D diagram 
illustrates the main features. From lWang et q/.ll2010l . 



trivial insulating states have been obtained from the com- 
bination of a tri vial noninterac ti ng band structure a nd in- 
teraction terms dGuo and Franzl l2009t iRaghu et all 12008 : 
Weeks and Franzl 120101) . Such topological insulators can be 
regarded as topologically nontrivial s tates arising from dy- 
namic ally generated SOC dWu et all l2007t IWu and Zhangj 
20041) . The effect of in teractions on the the QS H state has also 
been recently studied dRachel and Hurl [20101) . In 3D, strong 
topologic al insulators with topological excitations ha ve been 
obtained ( Pesin and Balents , 2010t Zhang et al. . 120091) . Topo- 
logical insula tors have been sugg ested to exist in transition 
metal oxides (IShitade et all [200% . where the correlation ef- 
fect is strong. It was also proposed that one coul d achieve the 



topol ogical insulator state in Kondo insulators (IDzero et al. 



2010). All these topological Mott insulator states can be un- 
derstood in the framework of the topological ord er parameter 



expre ssed in terms of the full Green's function (IWang et al 



2010h . Interaction-induced topological insulator states such as 



the to pological Mott insulators proposed in Ref. dRaghu et al 



2008b correspond to regions represented by the point B in 
Fig. [36] which has a trivial noninteracting unperturbed Hamil- 
tonian Hq(B) but acquires a nontrivial topological order pa- 
rameter due to the interaction part Hi(B) of the Hamilto- 
nian. The previously discussed topological order parameters 
are useful for determining the phase diagrams of interacting 
systems. 

For disordered systems, the topological order parameters 
are still applicable, provided that we use the disorder-averaged 
Green's functions. In this case, Fig. [36] can be regarded as 
a simple phase diagram of disordered systems, with g inter- 
preted as the disorder strength. The representative point B is 



induced TI state has been studied recently ( Groth et al. 


20091; Guo 


. 20101; Guo et all 2010t 


Imura et al. , 


2009; 


Jiang et all 


20091 iLi et a/.Ll2009llLoring and HastingsL 


2010; 


Obuse et al. 


, 20081 Olshanetskv et a/.Ll2010tlOstrovskv et al. 


2009: Shind 


ou and Murakami! 120091). Therefore, the topolog- 



ical order parameters previously discussed have the ability to 
describe both interacting and disordered systems. 



B. Topological band theory 

We shall now give a brief introduction to TBT. Even though 
this theory is only valid for non-interacting systems, it has 
become an important tool in the discovery of new topo- 
logical materials. Unfortunately, evaluating the Z2 invari- 
ants for a generic band structure is in general a difficult 
problem. Several approaches have bee n explored in the lit- 
erature includ i ng sp i n Chern numbers ( Fukui and Hatsugal 
20071; IProdanl 120091: ISheng et all 120061) . topo logical invari 
ants constructed from Bloch wave functions ( Fu and Kane , 



2006 
2009) 



Kane and Melel 120051: iMoore and Balentsl 120071: iRov 



and discrete indices calcul ated from single-pa rticle 
states at TRIM in the Brillouin zone dFu and Kanell2007l). We 
will f ocus on the last method for its simplicity dFu and Kane . 
20071) . 

This basic quantity in this approach is the matrix element of 
the TR ope rator T between stat es with TR conjugate momenta 
k and -k dFu and Kanell2006l) . 



B a p (k) = (-k,a|T|k,/3). 



(77) 



Since B a p is defined as a matrix element between Bloch states 
at TR conjugate momenta, it is expected that this quantity 
contains some information about the band topology of TR 
invariant topological insulators. At the TRIM Ti, B(k = 
Ti) is antisymmetric , so th at the following quantity can be 
defined dFu and Kanel 120061) : 



pf[s(r l )] 



(78) 



in which Pf stands for the Pfaffian of an antisymmetric ma- 
trix. Since Pf[B(T 4 )] 2 = det[B(r i )], we have Si = ±1. It 
should be noticed that the wavefunctions |k, a) must be cho- 
sen continuously in BZ to avoid ambiguity in the definition 
of 5i. In ID, there are on ly two TRIM, and a "TR polariza- 
tion" dFu and Kane , 120061) can be defined as the product of 
Si, 



7T 



-1 



which is 

tion ( King-Smith and Vanderbilt , 



= <5l^2, 

2 analog to the 



(79) 

ch arge po lariza 



19931: iRestaL 119941 



Thoulessl I1983T IZakl 119891) . A further analogy between 
the charge polarization and the TR polarization suggests 
the form of the Z2 invariant for TR invariant topological 
insulators. If an angular parameter 8 is tuned from to 2-7T, 
the change in the charge polarization P after such a cycle is 
expressed as the first Chern number C\ in the (k, 8 ) space . In 



fact, t he same C\ would gives the TKNN invariant dThouless , 
19831) if 8 were regarded as a momentum. By analogy with 
P, the Z 2 invariant for (2 + 1)D topological insulators can be 
defined as 

(_l)f2D _ ^_i^Pe(k 2 =0)-P B (k 2 =Ti) ^ ^g0) 

where Pe(k,2) = S1S2, Si is defined at the TRIM k\ = 0or7r, 
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FIG. 37 (a) The 2D bulk Brillouin zone projected onto the ID edge 
Brillouin zone. The two edge TRIM Ai and A2 are projections of 
pairs of the four bulk TRIM r i=(afl) . (b) Projection of the TRIM 
of the 3D Brillouin zone onto a 2D surface Brillouin zone. From 
iFu and Kanell2007l . 



and &2 is regarded as a parameter. Expanding Eq. ( f80b gives 

4 

( _ 1 )^ D= JJ ( j. j (81) 

i=l 

where i = 1, 2,3,4 labels the four TRIM in the 2D Bril- 
louin zone. (— \y 2D = +1 implies a trivial insulator while 
(— \y 2D = —1 implies a topological insulator. Furthermore, 
as a TR polarization, ^2D also determines the way in which 
Kramers pairs of surface states are connected [Fig. d38bl. 
which suggests that bulk topology and edge physics are inti- 
mately related. This is another example of the "holographic 
principle" for topological phenomena in condensed matter 
physics. 

We now discuss 3D topological insulators. It is interesting 
to note that in TBT the natural route is "dimensional increase", 
in contrast to the "dimensional reduction" procedure of TFT. 
From this dimensional increase , the 3D (strong) topological 



invari ant is naturally defined as dFu and KaneL 12007 Fli et all 
2007) 



(-1)" 



(82) 



In addition to the strong invariant, it has been shown that 
the product of any four Si's for which the I\ lie in the same 
plane is also gauge invariant, and defin es topologica l invari - 
ants characterizing the band structure dFu and Kanei 120071) . 
This fact leads to the definition of three additional invariants in 



3D known as weak topological invariants dFu and Kane U2007t 
Fu et all l2007t iMoore and Balentsl l2007t iRovl l2009h . These 
I2 invariants can be arranged as a 3D vector with elements Vk 
given by 



n 



5, = 



i=(nin 2 n 3 ) j 



(83) 



where {v\ v^) depend on the choice of reciprocal lattice vec- 
tors and are only strictly well defined when a well-defined lat- 
tice is present. It is useful to view these invariants as compo- 
nents of a mod 2 reciprocal lattice vector, 
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FIG. 38 Schematic representation of the surface energy levels of a 
crystal in either 2D or 3D, as a function of surface crystal momentum 
on a path connecting TRIM A a and At. The shaded region shows the 
bulk continuum states, and the lines show discrete surface (or edge) 
bands localized near one of the surfaces. The Kramers degenerate 
surface states at A a and Aj, can be connected to each other in two 
possible ways, shown in (a) and (b), which reflect the change in TR 
polarization APe of the cylinder between those points. Case (a) oc- 
curs in topological insulators, and guarantee s the surface ba nds cross 
any Fermi energy inside the bulk gap. From lFu et a/.L[2007l . 



When vq = 0, states are classified according to G u , and are 



2007) when the 



called weak topological insulators dFu et al. 
weak indices Vk are odd. 

Heuristically these states can be interpreted as stacked QSH 
states. As an example, consider planes of QSH stacked in the 
z direction. When the coupling between the layers is zero, the 
band dispersion will be independent of k z . It follows that the 
four Si's associated with the plane k z = ir/a will have prod- 
uct — 1 and will be the same as the four associated with the 
plane k z = 0. The topological invariants will then be given 
by vq = and G u = (2ir/a)z. This structure will remain 
when hopping between the layers is introduced. More gener- 
ally, when QSH states are stacked in the G, direction the in- 
variant will be G„ = G mod 2. This implies that QSH states 
stacked along differen t directions Gi and G2 are equivalent 
if Gj = G 2 mod 2 dFu and Kanel l2007h . As for the sur- 
face states, when the coupling between the layers is zero, it is 
clear that the gap in the 2D system implies there will be no 
surface states on the top and bottom surfaces; only the side 
surfaces will have gapless states. We can also think about the 
stability of the surface states for the weak insulators. In fact, 
weak topological insulators are unstable with respect to dis- 
order. We can heuristically see that they are less stable than 
the strong insulators in the following way. If we stack an odd 
number of QSH layers, there would at least be one delocalized 
surface branch. However, the surface states for an even num- 
ber of layers can be completely localized by disorde r or pertur- 



bation s. Despite this instability, it has been shown dRan et al. 



G„ = ^bi + i/ 2 b 2 + ^b 3 



(84) 



20091) that the weak topological invariants guarantee the exis- 
tence of gapless modes on certain crystal defects. For a dislo- 
cation with Burgers vector b it was shown that there will be 
gapless modes on the dislocation if G„ • b = (2?i + l)7r for 
integer n. 

Similar to the 2D topological insulator, there are connec- 
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FIG. 39 Diagrams depicting four different phases indexed by 
va; ivwivz)- The top panel depicts the signs of Si at the points I\ 
on the vertices of a cube. The bottom panel characterizes the band 
structure of a 001 surface for each phase. The solid and open circles 
depict the TR polarization n a at the surface momenta A a , which are 
projections of pairs of T, which differ only in their z component. 
The thick lines indicate p ossible Fermi arcs which enclose specific 
A . From lFuef flZll2007l . 



tions between the bulk invariants of 3D topological insulator 
and the corresponding 2D surface state spectrum. As a sam- 
ple, Fig. d39t shows four different topological classes for 3D 
band structures labeled with the corresponding (vq; v\ViVz). 
The eight I\ are represented as the vertices of a cube in mo- 
mentum space, with the corresponding Si shown as ± signs. 
The lower panel shows a characteristic surface Brillouin zone 
for a 001 surface with the four A a labeled by either filled or 
solid circles, depending on the value of 7r a = 5j=( a i)5j=( a 2V 
Genetically it is expected that the surface band structure will 
resemble Fig. l38l b) on paths connecting two filled circles or 
two empty circles, and will resemble Fig. l38la) on paths con - 
necting a filled circle to an empty circle dFu and Kanel 120071) . 
This consideration determines the 2D surface states qualita- 
tively. 

If an insulator has inversion symme try, there is a simpl e 
algorithm to calculate the Z2 invariant dFu and Kanel 120071) : 
indeed, the replacement in Eq. ( |8TT i and Eq. (l82l of Si by 



N 



Si = Y\ &m(Fi), 



(85) 



m— 1 



gives the correct Z 2 invariants. Here 62m = ±1 is the 
parity eigenvalue of the 2mth occupied energy band at r, 
[Fig. [37), which shares the same e igenvalue £ 2 m = 
with its degenerate Kramers partner ( Fu and Kanel 2007 ). The 
product is only over half of the occupied bands. Since the 
definition of the Si relies on parity eigenvalues, the Si are 
only we ll-defined in t his ca se when inversion symmetry is 
present dFu and Kanel 120071) . However, for insulators with- 
out inversion symmetry, this algorithm is very useful. In fact, 
if we can deform a given insulator to an inversion-symmetric 
insulator and keep the energy gap open along the way, the 
resultant Z2 invariants are the same as the initial ones due 
to t opological invarian ce, but can be calculated from par- 
ity dFu and Kanel 12007b . 



C. Reduction from topological field theory to topological 
band theory 

We now briefly discuss the relation between the TFT and 
TBT. On the one hand, the TFT approach is very powerful to 
reveal various aspects of the low-energy physics, and it also 
provide a deep understanding of the universality among dif- 
ferent systems. Furthermore, in contrast to TBT, TFT is valid 
for interacting systems. On the other hand, from a practical 
viewpoint, we also need fast algorithms to calculate topolog- 
ical invariants, which is the goal of TBT. An intuitive under- 
standing of the TBT of Z2 topological insulators is as follows. 
For integer-class CS topological insulators, the topological in- 
variant C n is expressed as the integral of Green's functions 
(or Berry curvature, in the noninteracting limit). Therefore, 
the knowledge of Bloch states over the whole Brillouin zone 
is needed to calculate C n . For TR invariant Z2 topological 
insulators, the TR symmetry constraint enables us to deter- 
mine the topological class of a given insulator with less in- 
formation: we do not need the information over the entire 
Brillouin zone. For insulators with inversion symmetry, the 
parity at several high -symmetry po i nts co mpletely determines 
the topological class dFu and Kanel 120071) . which explains the 
success of TBT. As naturally expected, the TBT approach is 
related to the TFT approa ch. In fact, it has been recently 
proved ( Wang et al. , 2010t) that the TFT description can be 
exactly reduced to the TBT in the nonint eracting limit. We 
now outline this proof dWang et aZl l2010b . Starting from the 
expression for P3 in Eq. d63l and Eq. d65l l. one can show that 

2P 3 (mod2) = / d 3 ke ijk Tr[(Bd i B r ){Bd j B r ) 



24tt 2 j 
•x{Bd k B^)] (mod 2). 



(86) 



By some topological argument, this expression for 
P3 is shown to g iv e the degree dcg/ of certain 
map dDubrovin et ail 1985) from the Brillouin zone 
three-torus T 3 to the SU(2) group manifold. There are two 
seemingly different expressions for deg/, one of which is 
of integral form as given by Eq. d86l ), while the other is 
of discrete form and given simply by the number of points 
mapped to a arbitrarily chosen image in SU(2). Due to TR 
symmetry, if we choose the image point as one of the two 
antisymmetric matrices in SU(2) (e.g. icr y ), we have an 
interesting " pair annihilation" o f those points other than the 
eight TRIM dWang et al.[ |2010|) . The final result is exactly 
the Z 2 invar iant from TBT. The explicit relation between TFT 
and TBT is dWang et al.i l2010b 



(-1) 



2P 3 



(-1)' 



(87) 



V. TOPOLOGICAL SUPERCONDUCTORS AND 
SUPERFLUIDS 

Soon after their discovery, the study of TR invariant topo 
logical insulators was generalized to TR invariant topologi 
cal superconductors and superfluids dKitaevl 120091 lOi et a 
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20091: IrovI 120081: ISchnvder et all I2008T) . There is a direct 
analogy between superconductors and insulators because the 
Bogoliubov-de Gennes (BdG) Hamiltonian for the quasipar- 
ticles of a superconductor is analogous to the Hamiltonian of 
a band insulator, with the superconducting gap corresponding 
to the band gap of the insulator. 

3 He-B is an example of such a topological superfluid state. 
This TR invariant state has a full pairing gap in the bulk, 
and g apless surface stat e s con s isting of a single Majorana 



cone dChung and ZhangL 120091: bi et all ExM iRoyl 12008 



Schnyder et al.i 2008h . In fact, the BdG Hamiltonian for 
3 He-B is identic al to the model Ha miltonian of a 3D topo- 
logical insulator dZhang et ali 120091) . and investigated exten- 



sively in Sec. IIII.AI In 2D, the classification of topological 
superconductors is very similar to that of topological insula- 
tors. TR_bj^aking_™pei^ndj^ by an inte- 
ger ( Read and Green , 2000 : Volovikt 1988 ). similar to quan- 
tum Hall insulators (Th ouless et aZ.[ll982l). while TR invariant 
supe r conductors are classifi e d (iKitaev , 120091: lOi et al.\ . 2009: 



RovL I2008t ISchnvder et all 120081) by a Z 2 i nvariant in ID 



and 2D, but by an in teger (Z) invariant in 3D (lKitaev , l2009t 



Schnvder et ai, 2008) 



Besides the TR invariant topological superconductors, 
the TR breaking topological superconductors have also at- 
tracted a lot of interest recently, because of their rela- 
tion with non-Abelian statistics and their potential applica- 
tion to topological quantum computation. The TR break- 
ing topological superconductors are described by an inte- 
ger AT The vortex of a topological superconductor with 
odd topological quantum nu mber Af carrie s an odd num- 
ber of Majorana zero modes (I Volovikl 1 19991). giving rise to 
non-Abelian statistics ( Ivanov , 200 it Read and Green , 120001) 
which could provide a plat form for topological quantum com- 



puting (|Nay^k_eLa/i|2008|). The simplest model for an Af = 1 
chiral topological superconductor is realized in the p x + ip y 
pairing state of spinless fermions ( Read and Green , 2000l) . A 



spinful version of the chi ral superconductor has been pre- 
dicted to exist in Sr2Ru04 ( Mackenzie and Maenoi 2003 ). but 
the experimental situation is far from definitive. Recently, 
several new proposals to realize Majorana fermion states in 
conventional supercond uctors have been investigated by mak- 



physical degrees of freedom, each chiral Majorana edge state 
has half the degrees of freedom of the chiral edge state of a 
QH system. Therefore, the chiral superconductor is the "min- 
imal" topological state in 2D. The analogy between a chiral 
superconductor and a QH state is illustrated in the upper pan- 
els of Fig. |40] Following the same analogy, one can consider 
the superconducting analog of QSH state — a "helical" su- 
perconductor in which fermions with up spins are paired in 
the p x + ip y state, and fermions with down spins are paired 
in the p x — ip y state. Such a TR invariant state has a full gap 
in the bulk, and counter-propagating helical Majorana states 
at the edge. In contrast, the edge states of the TR invariant 
topological insulator are helical Dirac fermions with twice 
the degrees of freedom. As is the case for the QSH state, a 
mass term for the edge states is forbidden by TR symmetry. 
Therefore, such a superconducting phase is topologically pro- 
tected in the presence of TR symmetr y, and can be described 
by a Z2 topological quantum number ( Kitaevi l2009t lOi et al 



2009llRo\(l2008tlSchnyder et a/.[l2008l) . The four types of 2D 
topological states of matter discussed here are summarized in 
Fig. |40] 

As a starting point, we first consider the Hamiltonian of the 
simplest nontr ivial TR breaking super conductor, the p + ip su- 
perconductor ( Read and Green . 2000l) for spinless fermions: 



H 




, (88) 



with e p = p 2 /2m — p and p± = p x ± ip y . In the weak 
pairing phase with p > 0, the p x + ip y chiral superconductor 
is known to have chiral Majorana edge states propagating on 
each boundary, described by the Hamiltonian 



-Hedge = ^2 V F kyip-k y ipk 



(89) 



ky>0 



ing use of stron g SOC dFu and Kanet l2008t bi et all l2010t H S Te d Bernevig et al. 



Sau et aZ.[|2oToh 



wh ere jj-k,, = if 1 ], is the quasiparticle creation opera- 
tor dR^adjindGreen , 120001) and the boundary is taken parallel 
to the y direction. The strong pairing phase /i < is triv- 
ial, and the two phases are separated by a topological phase 
transition at fx = 0. 

In the BHZ m odel for the QSH state in 
200l), if we ignore the coupling 



A. Effective models of time-reversal invariant 
superconductors 

The simplest way to understand TR invariant topological 
superconductors is through their analogy with topological in- 
sulators. The 2D chiral superconducting state is the supercon- 
ductor analog of the QH state. A QH state with Chern num- 
ber N has N chiral edge states, while a chiral superconductor 
with topological quantum number Af has Af chiral Majorana 
edge states. Since the positive and negative energy states of 
the BdG Hamiltonian of a superconductor describe the same 



terms between spin up and spin down electrons, the system is 
a direct product of two independent QH systems in which spin 
up and spin down electrons have opposite Hall conductance. 
In the same way, the simplest model for the topologically 
nontrivial TR invariant superconductor in 2D is given by the 
following Hamiltonian: 



H 



2 ^ 

v 



A*p_ 










with *(p) = 





( p 

T 






-A*p. 



\ 



*(90) 



From Eq. d90l ) we 



see that spin up (down) electrons form p x + ip y (p x — ip y ) 
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FIG. 40 (Top row) Schematic comparison of 2D chiral superconduc- 
tor and QH state. In both systems, TR symmetry is broken and the 
edge states carry a definite chirality. (Bottom row) Schematic com- 
parison of 2D TR invariant topological superconductor and QSH in- 
sulator. Both systems preserve TR symmetry and have a helical pair 
of edge states, where opposite spin states counter-propagate. The 
dashed lines show that the edge states of the superconductors are Ma- 
jorana fermions so that the E < part of the quasiparticle spectrum 
is redundant. In terms of the edge state degrees of freedom, we have 



symb olically QSH = (QHf 
From lOi et fl7U2009l 



(Helical SC) = (Chiral SC) 



Cooper pairs, respectively. Comparing this model Hamilto- 
nian (f90T > for the topological superconductor with the BHZ 
model of the HgTe topological insulator [Eq. [2], we first see 
that the term proportional to the identity matrix in the BHZ 
model is absent here, reflecting the generic particle-hole sym- 
metry of the BdG Hamiltonian for superconductors. On the 
other hand, the terms proportional to the Pauli matrices a a are 
identical in both cases. Therefore, a topological superconduc- 
tor can be viewed as a topological insulator with particle-hole 
symmetry. The topological superconductor Hamiltonian also 
has half as many degrees of freedom as the topological insu- 
lator. The model Hamiltonian d90l l is expressed in terms of 
the Nambu spinor 5*(p) which artificially doubles the degrees 
of freedom as compared to the topological insulator Hamilto- 
nian. Bearing these differences in mind, in analogy with the 
QSH system, we know that the edge states of the TR invariant 
system described by the Hamiltonian d90l l consist of spin up 
and spin down quasiparticles with opposite chiralities: 



tfcdgC = VFk y (^-VI^M - tP-kyi^kyi) ■ (91) 

fe H >0 

The quasiparticle operators tpk t> ^k 4. can be expressed in 
terms of the eigenstates Uk (%),Vk ( x ) of the BdG Hamilto- 
nian as 



ipkyf = J d 2 x (u ky (x)c t (x) + v ky (x)c\{x)^j , 
V>fc B i = / d 2 x (u*_ ky (x)c±(x) + v*_ ky (x)cl(x)j , 



from which the TR transformation of the quasiparti- 
cle operators can be determined to be Ti/jk^T^ 1 = 

ip-kyU Ti^k^T^ 1 = -V>-ft v t" In other words, 



(V'fc fj i) transforms under TR as a Kramers doublet, 
which forbids a gap in the edge state spectrum when TR 
is preserved by preventing the mixing of spin up and spin 
down modes. To see this explicitly, notice that the only fey- 
independent term that can be added to the edge Hamiltonian 
( |9T1 i is imJ2k ^-k fV'fc i> w i tn TO rea l- However, such a 
term is odd under TR, which implies that any backscattering 
between quasiparticles is forbidden by TR symmetry. The dis- 
cussion above is exactly parallel to the Z2 topological charac- 
terization of QSH system. In fact, the Hamiltonian (|9Qt has 
exactly the same form as the four-band effect ive Hamiltonian 



of the QSH effect in HgTe quantum wells dBernevig et al. 
20061) . The edge states of the QSH insulator consist of an 
odd number of Kramers pairs, which re main gapless un 



der any small TR in variant perturbation dWu et al. , l2006t 



Xu and Moore , 120061) . Such a "helical liquid" with an odd 
number of Kramers pairs at the Fermi energy cannot be real- 
ized in any bulk ID system, and can only appear holograph- 



ic ally as the edge theory of a 2D QSH insulator dWu et al. 



20061) . Similarly, the edge state theory Eq. d9TT > can be called 
a "helical Majorana liquid", and can only exist on the bound- 
ary of a Z2 topological superconductor. Once such a topolog- 
ical phase is established, it is robust under any TR invariant 
perturbations such as Rashba-type SOC and s-wave pairing, 
even if spin rotation symmetry is broken. The edge helical 
Majorana liquid can be detected by electric transport through 
a quant um point contact be tween two topological supercon- 



ductors dAsano et aZl l2010l) 



The 2D Hamiltonian d90b describes a spin-triplet pairing, 
the spin polarization of which is correlated with the orbital 
angular momentum of the pair. Such a correlation can be nat- 
urally generalized to 3D where spin polarization and orbital 
angular momentum are both vectors. The Hamiltonian of such 
a 3D superconductor is given by 



e P hx2 ia 2 a a A^ Pj 

h.C. -£ p l2x2 



(92) 



where we use 



c tp> c Ipi c t- 



different basis ^(p) = 
A QJ is a 3 x 3 matrix with 



p 

a = 1,2,3 and j = x, y, z. Interestingly, an example of such 
a Hamiltonian is given by the well-known 3 He-B phase, for 
which the order parameter A"-* is dete rmined by an orthogo- 
nal matrix A aj = Au aj , u e SO (3) dVollhardt and Wolfle . 
19901) . Here and below we ignore the dipole-dipole inter- 
action term dLeggett , 1975). since it does not affect any 
essential topological properties. Performing a spin rotation, 
A" 3 can be diagonalized to A" 3 = A<5 Q -' , in which case the 
Hamiltonian (l92l can be expressed as: 



H 



d 2 x^ 







-A* Pz 





-A* Pz 



A P+ 
-A Pz 

- £ p 




-A Pz 




(SB) 



/ 



Compared with the model Hami ltonian Eq. ( l3"TT l 
plest 3D topological insulators dZhang et all 2( 



for the sim- 
2009b . we see 
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that the Hamiltonian d93l ) has the same form as that for Bi2Se3 
(up to a basis transformation), but with complex fermions 
replaced by Majorana fermions. The kinetic energy term 
p 2 /2m — /i corresponds to the momentum dependent mass 
term TVf (p) = M — Bip 2 , — B2P 2 of the topological insulator. 
The weak pairing phase /.t > corresponds to the nontriv- 
ial topological insulator phase, and the strong pairing phase 
/i < corresponds to the trivial insulator. From this analogy, 
we see that the superconductor Hamiltonian in the weak pair- 
ing phase describes a topological superconductor with gapless 
surface states protected by TR symmetry. Different from the 
topological insulator, the surface states of the topological su- 
perconductor are Majorana fermions described by 



ifsurf = ■= v F ip T k (k x cr y - k y a x ) ifa, 



(94) 



with the Majorana condition ?/>_k = 0" x ^ k ■ We see that 
this Hamiltonian for the surface Majorana fermions of a topo- 
logical superconductor takes the same form as the surface 
Dirac Hamiltonian of a topological insulator in this special 
basis. However, because of the generic particle-hole symme- 
try of the BdG Hamiltonian for superconductors, the possi- 
ble particle-hole symmetry breaking terms for surface Dirac 
fermions such as a finite chemical potential is absent for sur- 
face Majorana fermions. Because of the particle-hole sym- 
metry and TR symmetry, the spin lies strictly in the plane 
perpendicular to the surface normal, and the integer wind- 
ing number of the spin around the momentum is now a well- 
defined quantity. This integer winding number gives a Z clas- 
sification of the 3D top ological superconductor (IKitaev . 120091 
Schnyder et all 120081) . The surface state remains gapless un- 



der any small TR invariant perturbation, since the only avail- 
able mass term m ^ k V'^k^V'k is TR odd. The Majorana 
surface state is spin-polarized, and can thus be detected by its 
special contribution to the spin relaxation of an electron on 
the surface of 3 He-B, similar to the measurement of electron 
spin correl ation in a solid state sys tem by nuclear magnetic 
resonance ( Chung and Zhan j, 20091) . 



B. Topological invariants 

From the discussion above, we see that the model Hamil- 
tonian for the topological superconductor is the same as that 
for the topological insulator, but with the additional particle- 
hole symmetry. The simultaneous presence of both TR and 
particle-hole symmetry gives a different classification for the 
2D and 3D topological superconductors, in that the 3D TR in- 
variant topological superconductors are classified by integer 
(Z) classes, and the 2D TR invariant topological supercon- 
ductors are classified by the Z2 classes. To define a n integer- 



valued topological invariant (ISchnyder et all 120081) . we start 



from a generic mean-field BdG Hamiltonian for a 3D TR in- 
variant superconductor, which can be written in momentum 




FIG. 41 Setting for detecting the Majorana surface states of the He- 
B phase, which consist of a single Majorana cone. When electrons 
are injected into 3 He-B, they exist as "bubbles". If the injected elec- 
trons are spin-polarized, the spin will relax by interaction with the 
surfa ce Majorana modes, and this relaxation is strongly anisotropic. 
From Chung and Zhang, 200^. 



space as 



In a different basis we have H = 'I'^-ffk^k with 



2\hu- «TA 







(95) 



In general, ij\ is a vector with TV components, and /ik and Ak 
are TV x TV matrices. The matrix T is the TR matrix satisfying 
TUkT = h T k , T 2 = -I and T+T = I, with I the iden- 
tity matrix. We have chosen a special basis in which the BdG 
Hamiltonian has a special off-diagonal form. It should be 
noted that such a choice is only possible when the system has 
both TR symmetry and particle-hole symmetry. These two 
symmetries also require TA fc to be Hermitian, which makes 
the matrix Tik + ^TA^ generically non-Hermitian. The matrix 
ft k + iT'A^ can be decomposed by a singular value decompo- 
sition as Tik + zTA^ = [Tj[Z?kVk with C/k, Vk unitary matri- 
ces and £>k a diagonal matrix with nonnegative elements. One 
can see that the diagonal elements of are actually the pos- 
itive eigenvalues of H^. For a fully gapped superconductor, 
£>k is positive definite, and we can adiabatically deform it to 
the identity matrix I without closing the superconducting gap. 
During this deformation, the matrix hu + iTA^ is deformed 
to a unitary matrix Qk = C^Vk € U(N). The integer-valued 
topological invariant characterizing topologi cal superconduc 
tors is defined as the winding number of Qk ( ISchnyder et al. 
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FIG. 42 (a,b) Superconducting pairing on two Fermi surfaces of a 
3D superconductor, (c) An example of 2D TR invariant topologi- 
cal supercond uctor, (d) ID TR invariant topological superconductor. 
Adapted from lOi et a/1.120101 . 



2008): 



N w 



24tt 2 



d 3 k e yfc Tr 



Q{diQ^QidjQ^Qid k Qy$6) 



We note that the topological invariant d96l ) is expressed as 
an integral over the entire Brillouin zone, similar to its coun- 
terpart for topological insulators. However, there is a key 
difference. Whereas the insulating gap is well defined over 
the entire Brillouin zone, the superconducting pairing gap in 
the BdG equation is only well defined close to the Fermi sur- 
face. Indeed, superconductivity arises from a Fermi surface 
instability, at least in the BCS limit. Therefore, one would 
like to define topological invariants for a topological super- 
conductor strictly in terms of Fermi surface quantities. The 
desired topological invariant can be obtained by reducing the 
wind ing number in E q. d96*l ) to a integral over the Fermi sur- 



face dOiefa/.ll2010l) 



w 



(97) 



where s is summed over all disconnected Fermi surfaces and 
sgn(<5 s ) denotes the sign of the pairing amplitude on the sth 
Fermi surface. C\ s is the first Chern number of the sth Fermi 
surface (denoted by FS S ): 



C\ s — — 

2tt 



dfl lj (dia S j(k) - dja S i(k)) 



(98) 



with a S i = —i (sk| d/dki |sk) the adiabatic connection de- 
fined for the band |sk) which crosses the Fermi surface, and 
dfl tJ the surface element 2-form of the Fermi surface. 

As an example, we consider a two-band model with nonin- 
teracting Hamiltonian = k 2 /2m — fi + ak ■ a, for which 
there are two Fermi surfaces with opposite Chern number 



C± = ±1 [Fig.@2ta),(b)]. If we choose A k = iA a v , which 
has the same 6 S for both Fermi surfaces, we obtain Nw = 
[Fig. 1427 b)!. If we instead choose Ak = iAo<j y a ■ k, we ob- 
tain Nw = 1 FFig. l42l a)l. In the latter case, if we take the 
a — > limit, we arrive at the result Nw = 1 f° r the 3 He-B 
pha se, which indi cates that 3 He-B is topologically nontriv- 



ial dVolovikl. 120031) . 

For 2D TR invariant superconductors, a procedure of di- 
mensional reduction leads to the following simple Fermi sur- 
face topological invariant: 



^2D =IJ(sgn(*.)) r 



(99) 



The criterion d99b is quite simple: a 2D TR invariant super- 
conductor is topologically nontrivial (trivial) if there is an odd 
(even) number of Fermi surfaces, each of which encloses one 
TR invariant point in the Brillouin zone and has negative pair- 
ing. As an example, see Fig. 1427c). where Fermi surfaces 2 
and 3 have negative pairing. Fermi surfaces 3 and 4 enclose 
an even number of TR invariant momenta, which do not affect 
the I2 topological invariant. There is only one Fermi surface, 
surface 2, which encloses on odd number of TR invariant mo- 
menta and has negative pairing. As a result, the Z2 topological 
invariant is ( — l) 1 = — 1. 

For ID TR invariant superconductors, a further dimensional 
reduction can be carried out to give 



(100) 



where s is summed over all the Fermi points between and 
7r. In geometrical terms, a ID TR invariant superconductor is 
nontrivial (trivial) if there is an odd number of Fermi points 
between and 7r with negative pairing. We illustrate this for- 
mula in Fig.l427d). where the sign of pairing on the red (blue) 
Fermi point is —1 (+1), so that the number of Fermi points 
with negative pairing is 1 if the chemical potential /i = \x\ or 
fjb = fj,2, and if yU = /j.^. The superconducting states with 
/x = Hi and fj, = fj,2 can be adiabatically deformed to each 
other without closing the gap. However, the superconductor 
with fi = /X3 can only be obtained from that with ^2 through 
a topological phase transition, where the pairing order param- 
eter changes sign on one of the Fermi points. It is easy to see 
from this example that there are two classes of ID TR invari- 
ant superconductors. 



C. Majorana zero modes in topological superconductors 

1 . Majorana zero modes in p + ip superconductors 

Besides the new TR invariant topological superconductors, 
the TR breaking topological superconductors have attracted a 
lot of interest because of their relevance to non-Abelian statis- 
tics and topological quantum computation. In a p + ip super- 
conductor described by Eq. (l88l , it can be shown that the core 
of a superconducting vortex contains a localized quasiparticle 
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with exactly zero energy dVolovik , Il999h . The correspond- 
ing quasiparticle operator 7 is a Majorana fermion obeying 
[7, H] = and 7^ = 7. When two vortices wind around 
each other, the two Majorana fermions 71, 72 in their cores 
transform nontrivially. Because the phase of the charge-2e 
order parameter winds by 2tt around each vortex, an elec- 
tron acquires a Berry phase of it when winding once around 
a vortex. Since the Majorana fermion operator is a superpo- 
sition of electron creation and annihilation operators, it also 
acquires a tt phase shift, i.e. a minus sign when winding 
around another vortex. Consequently, when two vortices are 
exchanged, the Majorana operators 71, 72 must transform as 
7i — > 72, 72 — > ~7i- The additional minus sign may be as- 
sociated with 71 or 72, but not both, so that after a full wind- 
ing we have 7! ( 2 ) — ^ —71(2)- Since two Majorana fermions 
71 and 72 define one complex fermion operator 71 + 272, 
the two vortices actually share two internal states labeled by 
*7i72 = ±1- When there are 2N vortices in the system, the 
core states span a 2 JV -dimensional Hilbert space. The braid- 
ing of vortices leads to non-Abelian unitary transformations 
in this Hilbert space, implyin g that the vortices in this system 
obey non-Abelian statistics ( llvanovi 1200 It iRead and Green , 



2000). Because the internal states of the vortices are not local- 



ized on each vortex but shared in a nonlocal fashion between 
the vortices, the coupling of the internal state to the environ- 
ment is exponentially small. As a result, the superposition 
of different internal states is immune to decoherence, which 
is ideal for the purpose of quantum computation. Quantum 
computation with topologically protected q-bits is generally 
known as topological qua ntum computation, and is currently 



an active field of research (N avak et ail 2008). 



Several experimental candidates for p-wave superconduc- 
tivity have been proposed, among which Sr2Ru04, which 
is considered as the most promising candidate for 2D chi- 
ral superconductivity ( Mackenzie and Maenol 2003 ). How- 
ever, many properties of this system remain unclear, such as 
whether this superconducting phase is gapped and whether 
there are gapless edge states. 



Fortunately, there is an alternate route towards topological 
supercond uctivity without p-wave pairing. In 1981, Jackiw 
and Rossi (IJackiw and Rossilll98ll) showed that adding a Ma- 
jorana mass term to a single flavor of massless Dirac fermions 
in (2 + 1)D would lead to a Majorana zero mode in the vortex 
core. Such a Majorana mass term can be naturally interpreted 
as the pairing field due to the proximity coupling to a conven- 
tional s-wave superconductor. There are now three different 
proposals to realize this route towards topological supercon- 
ductivity: the superconducting proximity effect on the 2D sur 
face state of the 3D topological insulator ( Fu and Kane), 



on the 2D TR breaking topological insulator dOi et al 



and o n semiconductors with strong Rashba SOC ( ISau et al. 



2008) 



2010) 



20101) . We shall review these three proposals in the following. 



2. Majorana fermions in surface states of the topological 
insulator 



Fu and Kane dFu and KaneL 12008) proposed a way to real- 
ize the Majorana zero mode in a superconducting vortex core 
by making use of the surface states of 3D topological insula- 
tors. Consider a topological insulator such as Bi2Se3, which 
has a single Dirac cone on the surface with Hamiltonian from 
Eq.El 



H = V' t [v(<r x p) • Z - fj] ip, 



(101) 



where ip = (V't' i'i) an( l we nave taken into account a fi- 
nite chemical potential [i. Consider now the superconduct- 
ing proximity effect of a conventional ,s-wave superconduc- 
tor on the 2D surface states, which leads to the pairing term 
Ha = AV'jV'l + h.c. The BdG Hamiltonian is given by 

#BdG = |E P &H P V, where #t = ( ^ ^ 



and 



H 



I v(cr x p) • z 

= \ -ia v A* 



— v(cr X p) • z + /! 



The vortex core of such a superconductor has been shown to 
have a sing le Majorana zero mode, similar to a p + ip su per- 
conductor dFu and Kanel l2008t Ijackiw and Rossi 119811) . To 
understand this phenomenon, one can consider the case of fi- 
nite fi, and introduce a TR breaking mass term ma z in the 
surface state Hamiltonian ( 1 1 li t. As discussed in Sec. IIII.Bl 
this opens a gap of magnitude \m\ on the surface. Consider- 
ing the case /i > m > 0, // — m <C m, the Fermi level in the 
normal state lies near the bottom of the parabolic dispersion, 
and we can consider a "nonrelativistic approximation" to the 
massive Dirac Hamiltonian, 

H = ipi [v(er x p) • z + ma z — //] ip 



d 2 xip\ [ ^—+m-n) ip+, 
* 2m 



(102) 



where ip + is the positive energy branch of the surface states. 
In momentum space, ip+ p = u p^Pt + v p^Pl w i m u p = 

. Consider- 



and v p = S± 



2^p2 +m 2 "P - Ipl Y 2 2 ^~i 

ing the projection of the pairing term Ha onto the ip + band, 
we obtain 

H A ^Y1 ^+, P V-'+,-p Au p' i; p + h - c - 



2m 



(103) 



We see that in this limit, the surface Hamiltonian is the same 
as that of a spinless p + ip superconductor [Eq. d88ll. When 
the mass m is turned on from zero to a finite value, it can be 
shown that as long as m < \i, the superconducting gap near 
the Fermi surface remains finite, so that the Majorana zero 
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mode we obtained in the limit < fi — m <C m must re- 
main at zero energy for the original m = system. Once we 
have shown the existence of a Majorana zero mode at finite [i, 
taking the fj, — > limit for a finite A also leaves the super- 
conducting gap open, so that the Majorana zero mode is still 
present at fj, = 0. 

From the analogy with the p + ip superconductor shown 
above, we also see that the non-Abelian statistics of vortices 
with Majorana zero modes apply to this new system as well. A 
key difference between this system and a chiral p + ip super- 
conductor is that the latter necessarily breaks TR symmetry 
while the former can be TR invariant. Only a conventional 
s-wave superconductor is required to generate the Majorana 
zero modes in this proposal and in the other proposals dis- 
cussed in the following subsection. This is an important ad- 
vantage compared to previous proposals requiring an uncon- 
ventional p + ip pairing mechanism. 

There is also a lower dimensional analog of this nontrivial 
surface state superconductivity. When the edge states of 2D 
QSH insulator are in proximity with an s-wave superconduc- 
tor and a ferromagnetic insulator, one Majorana fermion ap- 
pears at each domain wall between ferromagn etic region and 
superconducting region f (|Fu and Kanel l2009h ). The Majo- 
rana fermion in this system can only move along the ID QSH 
edge, so that non-Abelian statistics is not well-defined. Be- 
cause an electron cannot be backscattered on the QSH edge, 
the scattering of the edge electron by a superconducting re- 
gion indu ced by proximity effect is always perfect Andreev 
reflection ( Adroguer et ah . 20ld : Guigou and Cayssol 2010t 
Sato etaimOldf) . 



(a) 




(b) 



QAH 
m 



TSC,JV=1 
TSC, N=2 

Metal 

Metal 




NSC, N= 



NI 



FIG. 43 (a) Phase diagram of the QAH-superconductor hybrid sys- 
tem for p — 0. m is the mass parameter, A is the magnitude of the 
superconducting gap, and M is the Chern number of the supercon- 
ductor, which is equal to the number of chiral Majorana edge modes, 
(b) Phase diagram for finite /i, shown only for A > 0. The QAH, 
normal insulator (NI) and metallic (M etal) phases are well-defined 
only for A = 0. From lOi et aZl^Old . 



3. Majorana fermions in semiconductors with Rashba 
spin-orbit coupling 

From the above analysis, we see that conventional s-pairing 
in the surface Hamiltonian d 1 1 b induces topologically non- 
trivial superconductivity with Majorana fermions. There is a 
2D system which is described by a Hamiltonian very similar 
to Eq. ( 11011 ), i.e. a 2D electron gas with Rashba SOC. The 

Hamiltonian is H = j d 2 x ip^ (j^i + a ( cr x P) ' ^ — ty ^' 
which differs from the surface state Hamiltonian only by 
the spin-independent term p 2 /2m. Consequently, when 
conventional s-wave pairing is introduced, each of the two 
spin-split Fermi surfaces forms a nontrivial superconductor. 
However, the Majorana fermions from these two Fermi sur- 
faces annihilate each other so that the s-wave superconduc- 
tor in t he Rashba system is trivial. It was pointed out re- 
cently ( Sau et all l2010h that a nontrivial superconducting 
phase can be obtained by introducing a TR breaking term 
Ma z into the Hamiltonian, which splits the degeneracy near 
k = 0. If the chemical potential is tuned to < \M\, the 
inner Fermi surface disappears. Therefore, superconductivity 
is only induced by pairing on the outer Fermi surface, and be- 
comes topologically nontrivial. Physically, one cannot induce 
a TR breaking mass term by applying a magnetic field in the 



perpendicular direction, because the magnetic field may de- 
stroy superconductivity. Two ways to realize a TR breaking 
mass term have been proposed: by applying an in-pla ne mag- 
netic field and making use of the Dresselhaus SOC dAlicea , 
20101) . or by exchange coupling to a ferromagnetic insulating 
layer ( Sau et al. , 2010l) . The latter proposal requires a het- 
erostructure consisting of a superconductor, a 2D electron gas 
with Rashba SOC, and a magnetic insulator. 

This mechanism can also be generalized to the ID semicon- 
ductor wires wi th Rashba SOC coupling in proximity w ith a 
superconductor dOreg et 0/1 120101: IWimmer et aZL l2010h . De- 
spite the ID nature of the wires, non-Ab elian statistics is stil l 
possible by making use of wire networks (lAlicea et 



4. Majorana fermions in quantum Hall and quantum 
anomalous Hall insulators 

More recently, a new approach to reali ze a top o logica l 



superconductor phase has been proposed ( Oi et al. , 2010h 



which is based on the proximity effect to a 2D QH or QAH 
insulator. Integer QH states are classified by an integer N 
corresponding to the first Chern number in momentum space 
and equal to the Hall conductance in units of e 2 /h. Con- 
sider a QH insulator with Hall conductance Ne 2 /h in close 
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proximity to a superconductor. Even if the pairing strength 
inducing by the superconducting proximity effect is infinites- 
imally small, the resulting state is topologically equivalent to 
a chiral topological superconductor with Z topological quan- 
tum number Af = 2N. An intuitive way to understand 
such a relation between QH and topological superconducting 
phases is through the evolution of the edge states. The edge 
state of a QH state with Chern number N = 1 is described 
by the effective ID Hamiltonian fledge = J2 P v Py r lpy r lpy> 
where 77+ , r\ Vy are creation/annihilation operators for a com- 
plex spinless fermion. We can decompose rj Py into its real 
and imaginary parts, r, Py = l/V2(-y Pyl + ij Py2 ) and rj^ = 
l/v / 2(7-p B i - ij-p y 2), where j Pya are Majorana fermion 
operators satisfying 7^ = 7 _ P!ja and {"f-p y a,l p ' yb } = 
S a bS„ „' • The edge Hamiltonian becomes 



edge mode vortex 



fledge 



E 

Py>0 



Py (7-p y i7p v i + l~ Py 2l Py 2) 



(104) 



up to a trivial shift of the energy. In comparison with the edge 
theory of the chiral topological superconducting state, the QH 
edge state can be considered as two identical copies of chiral 
Majorana fermions, so that the QH phase with Chern number 
N = 1 can be considered as a chiral topological supercon- 
ducting state with Chern number Af = 2, even for infinitesi- 
mal pairing amplitudes. 

An important consequence of such a relation between QH 
and topological superconducting phases is that the QH plateau 
transition from N = 1 to N = will generically split into 
two transitions when superconducting pairing is introduced. 
Between the two transitions, there will be a new topologi- 
cal superconducting phase with odd winding number Af = 1 
[Fig. @3). Compared to other approaches, the emergence of 
the topological superconducting phase at a QH plateau tran- 
sition is determined topologically, so that this approach does 
not depend on any fine tuning or details of the theory. 

A natural concern raised by this approach is that the strong 
magnetic field usually required for QH states can suppress su- 
perconductivity. The solution to this problem can be found 
in a special type of QH state — the QAH state, which 
is a TR breaking gapped state with nonzero Hall conduc- 
tance in the absence of an external orbital magnetic field 
(Sec. III. El l. There exist now two realistic proposals for real- 
izing the QAH state experimentally, both of which make use 
of the TR i nvariant topologi cal insulator materials Mn-doped 
HgTe OWs dLiu et a/.Ll2008h . and Cr- or Fe-doped Bi 2 Se 3 thin 
films (lYu et ali 120101) . The latter material is proposed to be 
ferromagnetic, and can thus exhibit a quantized Hall conduc- 
tance at zero magnetic field. The former material is known to 
be paramagnetic for low Mn concentrations, but only a small 
magnetic field is needed to polarize the Mn spins and drive 
the system into a QAH phase. This requirement is not so pro- 
hibitive, because a nonzero magnetic field is already neces- 
sary to generate superconducting vortices and the associated 
Majorana zero modes. 




FIG. 44 3D topological insulator in proximity to ferromagnets 
with opposite polarization (Mf and M4.) and to a superconductor 
(S). The top panel shows a single chiral Majorana mode along the 
edge between superconductor and ferromagnet. This mode is elec- 
trically neutral, and therefore cannot be detected electrically. The 
Mach-Zehnder interferometer in the bottom panel converts a charged 
current along the domain wall into a neutral current along the su- 
perconductor (and vice versa). This allows electrical detection of 
the parity of the num ber of enclosed vortices/flux quanta. From 
lAkhmerovef aUl2009l . 



5. Detection of Majorana fermions 

The next obvious question is how to detect the Ma- 
jorana fermion if such a proposal is experimentally re- 
alized. There exist two similar theoretical proposals of 
electrical transport measurements to detect these M ajorana 
fermions ( Akhmerov et al. . l2009t iFu and KaneL 120091) . Con- 
sider the geometry shown in Fig. [44] This device is a com- 
bination of the inhomogeneous structures on the surface of 
a topological insulator discussed in the previous subsections. 
The input and output of the circuit consist of a chiral fermion 
coming from a domain wall between two ferromagnets. This 
chiral fermion is incident on a superconducting region where 
it splits into two chiral Majorana fermions. The chiral Ma- 
jorana fermions then recombine into an outgoing electron or 
hole after traveling around the superconducting island. More 
explicitly, an electron incident from the source can be trans- 
mitted to the drain as an electron, or converted to a hole by 
an Andreev process in which charge 2e is absorbed into the 
superconducting condensate. To illustrat e the idea we discus s 
the behavior for a E = quasiparticle dFu and Kane , l2009h . 
A chiral fermion incident at point a meets the superconduc- 



tor and evolves from an electron c\ into a fermion i\> built 
from the Majorana operators 71 and 72. The arbitrariness in 
the sign of 71.2 allows us to choose ip = 71 + 272. After 
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the quasiparticle winds around the superconducting region, ip 
recombines into a complex fermion at point d. This fermion 
must be either or cj, since a superposition of the two is not a 
fermion operator and is thus forbidden. To determine the cor- 
rect operator we can use adiabatic continuity. When the size 
of the superconductor shrinks continuously to zero, points a 
and d continuously tend to each other. Adiabatic continuity 
implies that an incident E = electron is transmitted as an 
electron, cj — > c\. However, if the ring encloses a quantized 
flux $ = nhc/2e, this adiabatic argument must be reconsid- 
ered. When n is an odd integer, the two Majorana fermions 
acquire an additional relative phase of ir, since each flux quan- 
tum hc/2e is a ir flux for an electron, and thus ir for a Majo- 
rana fermion. Up to an overall sign, one can take 71 — > —71 
and 72 —J- 72- Thus, when the ring encloses an odd number of 
flux quanta, c\ — > eg, and an incident E = electron is con- 
verted to a hole. The general consequences of this were cal - 
culated in detail dAkhmerov et al. l2009tlFu and Kanell2009l) . 
and it was shown that the output current (through arm d in the 
lower panel of Fig.l44li changes sign when the number of flux 
quanta in the ring jumps between odd and even. This unique 
behavior of the current provides a way to electrically detect 
Majorana fermions. 

Besides these two proposals reviewed above, several other 
theoretical proposals have also been made recently to ob- 
serve the Majorana fermion state, whic h make use of 
the Coulomb chargi ng energy © 120101) or a flux qubit 



( Hassleref al. , 2010h . More indirectly, Majorana fermions 



can also be detected through their contribution to Josephson 
coupling dFu and Kanel 120081 l2009t iLinder and Sudbol |20 1 ' 



Linder eVa/.H2010l:lLutchvn et a/.ll2010l:lTanaka et q/.ll200^ 
For a topological superconductor ring with Majorana fermions 
at both ends, the period of Josephson cu rrent is doubled, in- 
dependent from the physical realiz ation (IFu and Kane , l2009t 
Kiteevll200lllLutchvn et aZ.ll2010l) . 



VI. OUTLOOK 

The subject of topological insulators and topological super- 
conductors is now one of the most active fields of research in 
condensed matter physics, developing at a rapid pace. The- 
orists have systemati cally cla s sified topological states in all 



dimensions. Ref. (lOi et ail 120081) initiated the classifica 



tion program of all topological insulators according to discrete 
particle-hole symmetry and the TR symmetry, and noticed a 
periodic structure with period eight, which is known in math- 
ematics as the Bott periodicity. More extended and systematic 
classification of all topological insulator and superconductor 
states are obt ained according to TR, particle - hole and b i partite 
symmetries dKitaevl [20091: lOi et all l2008t iRvu eToZl koiot 



Schnvder era/., 2008 



Stone et al. , 2010). Such classification 



scheme gives a "periodic table" of topological states, which 
may play a similar role as the familiar period table of ele- 
ments. For future progress on the theoretical side, the most 
important outstanding problems include interaction and dis- 



order effects, realistic predictions for topological Mott insu- 
lator materials, a deeper understanding of fractional topolog- 
ical insulators and realistic predictions for materials realiza- 
tions of such states, the effective field theory description of 
the topological superconducting state, and realistic materials 
predictions for topological superconductors. On the experi- 
mental side, the most important task is to grow materials with 
sufficient purity so that the bulk insulating behavior can be 
reached, and to tune the Fermi level close to the Dirac point 
of the surface state. Hybrid structures between topological in- 
sulators and magnetic and superconducting states will be in- 
tensively investigated, with a focus on detecting exotic emer- 
gent particles such as the image magnetic monopole, the axion 
and the Majorana fermion. The theoretical prediction of the 
QAH state is sufficiently realistic and its experimental discov- 
ery appears to be imminent. The topological quantization of 
the TME effect in 3D and the spin-charge separation effect 
in 2D could experimentally determine the topological order 
parameter of this novel state of matter. 

Due to space limitations, we did not discuss in detail the 
potential for applications of topological insulators and super- 
conductors. It would be interesting to explore the possibility 
of electronic devices with low power consumption based on 
the dissipationless edge channels of the QSH state, spintron- 
ics devices based on the unique current-spin relationship in 
the topological surface states, infrared detectors, and thermo- 
electric applications. Topological quantum computers based 
on Majorana fermions remain a great inspiration in the field. 

Topological insulators and superconductors offer a platform 
to test many novel ideas in particle physics — a "baby uni- 
verse" where the mysterious 9 vacuum is realized, where ex- 
otic particles roam freely and where compactified extra di- 
mension can be tested experimentally. In the introduction to 
this article we drew an analogy between the search for new 
states of matter and the discovery of elementary particles. Up 
to now, the most important states of quantum matter were first 
discovered empirically and often serendipitously. On the other 
hand, the Einstein-Dirac approach has been most successful 
in searching for the fundamental laws of nature: pure logical 
reasoning and beautiful mathematical equations guided and 
predicted subsequent experimental discoveries. The success 
of theoretical predictions in the field of topological insulators 
shows that this powerful approach works equally well in con- 
densed matter physics, hopefully inspiring many more exam- 
ples to come. 
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